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1.  INTRODUCTION 

Station  02  -  the  main  streamflow  gauging  station  of  the 
Nahal  Yael  watershed  -  has  recorded  eight  flows  during 
the  14  runoff  events  that  occurred  in  the  catchment  sin¬ 
ce  the  inception  of  the  program  in  1965.  The  six  remai¬ 
ning  flows  were  not  large  enough  to  traverse  the  mid¬ 
basin  alluvial  channel  reach,  into  which  they  infiltrated 
The  last  of  these  eight  "throughf lows"  took  place  as  a  re 
suit  of  event  12,  on  20-21  February  1975.  Since  then, 
for  a  period  of  f i ve-and-a-hal f  years,  not  a  drop  of 
floodwater  has  wetted  the  weir  of  station  02. 

This  example  is  one  of  many  that  could  be  presented  for 
the  extreme  vari abi 1 i ty  i n  time  of  phenomena  which  con¬ 
trol  the  mechanics  of  arid  zone  geosystems.  We  have  for 
Nahal  Yael  a  data  bank  which,  however,  valuable  and 
unique,  parallels  -  from  several  points  of  view  -  a 
record  of  one  or  two  years'  duration  only  in  a  non-arid 
watershed.  This  has  been  expected,  of  course,  at  the 
time  the  program  originated  15  years  ago,  but  the  extre¬ 
mity  of  the  behavior  could  not  have  been  fully  appre¬ 
ciated.  It  is  clear  now  that  a  period  of  a  few  decades 
is  the  minimum  period  of  observation  for  supplying  a 
reasonable  measurement  record  of  cumulatively  dominant 
terrain  processes  in  the  extreme  desert. 

Given  this  situation,  the  other  approach  possible  for 
understanding  quantitatively  the  mechanics  of  hydro- 
geomorphic  watershed  response  on  an  event-to-event  sca¬ 
le  is  a  strictly  deterministic  process  approach.  This 
can  be  accomplished  only  on  the  basis  of  responses  of 
small  unit  areas  during  short  time  intervals,  with  re¬ 
sults,  subsequently  appropriately  integrated  over  time 
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and  space,  meaningfully  compared  with  data  measured 
on  the  watershed  scale. 

The  report  herewith  presented  reflects  the  transition 
of  the  emphasis  of  the  Nahal  Yael  Research  Project 
from  the  direct  measurement  approach  of  the  late  six¬ 
ties  and  early  seventies,  to  the  deterministic  model¬ 
ling  approach  of  recent  years,  which  we  hope  to  de¬ 
velop  in  the  early  eighties.  As  in  the  earlier  major 
report  (Schick  &  Sharon,  1974),  we  elected  also  in  the 
present  report  to  include  certain  studies  in  detail, 
while  others  are  treated  only  in  summary  form.  In  this 
sense  the  central  parts  of  this  report  are  chapter  8, 
which  presents  an  areal  framework  for  the  development 
of  a  geosystem-oriented,  distributed,  deterministic 
model,  and  chapter  9,  which  summarizes  an  innovative 
series  of  field  experiments  designed  to  functionally 
correlate  infiltration  parameters  with  terrain  proper¬ 
ties. 

The  "instrumental"  part  of  the  work  is  summarized  in 
chapter  2  and  3.  As  shown  in  chapter  2  the  last  few 
years  were  lean  ones  in  terms  of  major  flood  events. 
This  fact,  however,  gave  us  the  opportunity  to  complete 
the  design,  construction,  and  instrumentation  of  the 
Nahal  Yael  dam,  over  a  period  of  more  than  a  year,  with 
out  the  loss  of  either  data  or  construction  equipment 
andmaterial.  The  dam,  described  in  chapter  3,  is  expec¬ 
ted  to  impound  total  flow  volumes  and  total  sediment 
yields  for  all  events  up  to  a  near-catastrophic  magni¬ 
tude,  and  thus  will  considerably  enhance  the  credibili¬ 
ty  of  the  water  and  sediment  data  recorded  since  1965. 

Work  more  oriented  toward  modelling  the  watershed  pro¬ 
cesses  is  summarized  in  chapter  4,  which  presents  a  re¬ 
lationship  between  rainfall  energy  and  watershed  sedi¬ 
ment  yield,  and  in  chapter  6,  which  discusses  the  bed- 
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load  transport  problem  in  arid  streams.  The  results  of 
an  attempt  at  conventional,  simplified  lumped  modelling 
for  runoff  is  summarized  in  chapter  10.  Chapter  7  is  a 
brief  report  on  a  geophysical  survey  conducted  along 
the  main  alluvial  reach  of  the  Naha!  Yael  watershed  in 
order  to  obtain  depth  to  bedrock,  and  through  it,  an 
estimate  of  its  water  storage  capacity  and  flood  ab¬ 
straction  potential.  Finally,  chapter  5  presents  the 
initial  results  of  our  search  for  a  method  to  quantita¬ 
tively  represent  such  small-unit  terrain  characteristics 
which  affect  the  parameters  of  the  Nahal  Yael  arid  zone 
geosystem  model . 

Contrary  to  many  other  modelling  studies,  which  depend 
heavily  on  sophisticated  mathemati cal -hydraul ic  exer¬ 
cises  with  but  little  organic  input  of  real-world  con¬ 
ditions  and  processes,  the  study  reported  here  is  essen¬ 
tially  field-based.  Although  an  effort  is  continuously 
made  to  work  with  general  principles  and  universally 
applicable  techniques  (such  as,  for  example,  the  areal 
framework  and  unit  terrain  characterization),  the  op¬ 
portunity  of  constant  intercourse  with  a  real-life,  ope¬ 
rating,  and  process-oriented  instrumented  watershed 
while  developing  the  model  is  regarded  by  us  as  a  defi¬ 
nite  advantage.  We  hope  to  be  able  to  continue  this  dual 
task  in  future. 


2.  BASIC  DATA 


The  Nahal  Yael  watershed  instrumental  network  continued 
to  function  along  the  lines  indicated  in  the  previous 
report  (Schick  &  Sharon,  1974),  with  the  following 
changes : 

Rainfall:  Four  recorders  are  now  regarded  as  satisfacto¬ 
rily  representative  of  the  catchment  rainfall,  as  conclu- 
dedon  the  basis  of  the  elaborate  network  maintained  for 
a  number  of  years  in  the  past.  Rainfall  values  applica¬ 
ble  to  the  small  upstream  Nahal  Yael  catchments  04  and 
05  are  adequately  represented  by  averaging  one  hilltop 
and  one  valley  station. 

Instrument  30  had  to  be  moved  a  few  meters  uphill,  as 
it  stood  originally  too  close  to  the  rock-cut  central 
part  of  the  reservoir  overflow  channel  (see  chapter  3 
below) . 

Considerable  use  was  made  of  information  available  from 
radar  screens  with  regard  to  alerts  for  probable  events 
and  decisions  as  to  organizing  flood  survey  teams. 

Runoff:  Station  01  at  the  foot  of  the  Nahal  Yael  allu¬ 
vial  fan  near  Nahal  Roded  confluence  was  discontinued 
with  the  beginning  of  the  construction  of  the  dam.  Its 
instrument  was  moved  to  the  upstream  side  of  the  dam 
and  installed  as  a  reservoir  water  level  gauge  (see 
chapter  3  below). 

The  hydraulic  controls  associated  with  the  gauging  sta¬ 
tions  of  the  Nahal  Yael  watershed  have  gradually  silted 
up  to  an  unacceptable  degree  during  their  decade  of 
operation.  On  several  occasions  in  the  past,  sediment 
from  the  immediate  vicinity  of  the  connecting  intake 
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pipes  was  cleared,  but  the  continued  piling  up  of  sediment 
in  the  larger  area  behind  the  weirs  posed  a  severe  pro¬ 
blem.  Plans  for  the  removal  of  this  sediment  and  the 
restoration  of  the  original  stilling  basin  properties 
were  complicated  by  the  inaccessibility  of  all  four 
sites  for  vehicles  or  for  any  other  mechanical  equip¬ 
ment  too  heavy  to  carry  up  the  dam  and  the  rock  water¬ 
fall  by  hand. 

In  July-August  1978  the  sediment  clearing  was  accom¬ 
plished  by  a  group  of  four  student-assistants.  About 
50  cubic  meters  of  sediment  were  excavated  by  shovel¬ 
ling  according  to  predetermined  specifications,  hauled 
by  wheelbarrow  to  the  weir  and  dumped  down  the  down¬ 
stream  side  of  the  control.  Thereby  the  gauging  sta¬ 
tions  have  been  restored  to  their  original  level  of 
accuracy  at  least  for  the  few  coming  events. 

Sediment;  Automatic  suspended  samplers  71  and  72,  pre¬ 
viously  located  on  or  downstream  of  the  dam  site,  have 
been  discontinued.  The  middle  and  lower  bedload  traps 
have  also  been  dismantled. 

In  line  with  our  increased  interest  in  the  small  up¬ 
stream  watersheds  due  to  their  importance  for  evalua¬ 
ting  watershed  processes,  two  new  automatic  suspended 
sediment  samplers  (Nos.  78  and  79)  were  installed  in 
1979  in  tributaries  of  watershed  04. 

N§nC§ti on_of _event§ ;  Since  event  12,  which  occurred  in 
February  1975,  a  period  of  "non-flows"  set  in.  This  was 
nearly  broken  early  in  1977,  when  on  one  occasion  run¬ 
off  threshold  was  closely  approached,  but  not  actually 
reached . 

On  7  February  1978  simulatneous  reports  were  received 
from  several  voluntary  informers  and  through  the  radio 
of  a  storm  which  hit  Eilat  and  yielded  14  mm  of  rainfall. 
On  the  assumption  that  such  a  considerable  amount  of 
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rainfall  (nearly  one  half  of  the  mean  annual  amount)  is 
bound  to  have  some  regional  influence,  and  therefore 
also  affect  Nahal  Yael  3  km  away,  a  flood  survey  team 
of  six  was  quickly  organized.  The  team  left  Jerusalem 
within  a  few  hours  of  the  event,  but  upon  reaching  the 
watershed  was  disappointed  to  find  that,  due  to  a 
particularly  sharp  rainfall  gradient  which  existed  bet¬ 
ween  Eilat  and  Nahal  Yael,  only  3  mm  fell  in  the  latter. 
This  amount  was  insufficient  to  generate  any  flow  ex¬ 
cept  on  a  few  rocky  slopes. 

Event  13  finally  arrived  on  11  December  1978.  It  began 
with  a  few  intermittent  rains  of  low  intensity  (0.05  mm/ 
min  and  less).  Gradually  the  rains  became  stronger,  but 
the  peak  intensities  for  the  3-hour  rain  duration  rea¬ 
ched  only  medium-low  values:  0.33  mm/min  for  station 
24;  0.27  mm/min  for  station  30.  These  values  obtain  for 
periods  of  2-3  minutes  only;  hence  the  total  storm  rain¬ 
fall  of  11.0  mm  (a  mean  value  between  one  hilltop  gauge 
and  one  valley  gauge)  fell  in  an  exceptionally  tempe¬ 
rate  and  inerosive  manner. 

The  response  of  the  watershed  to  this  kind  of  rainfall 
is  of  particular  interest.  Despite  the  large  rainfall 
amount,  which  in  all  earlier  events  did  not  fail  to 
yield  at  least  some  flow  at  the  downstream  stations 
(03,  02),  only  the  upstream  stations  recorded  a  flow. 

The  main  data  are  listed  in  Table  2.1. 

The  peak  discharges  are  very  low.  All  of  the  conside¬ 
rable  volume  of  water  percolated  into  the  alluvial  rea¬ 
ches  downstream  of  station  04,  as  well  as  in  the  short 
alluvial  reach  upstream  of  station  03.  The  capacity  of 
the  alluvial  fill  to  absorb  the  floodwaters  must  have 
been  close  to  the  absolute  maximum,  as  nearly  four  full 
years  have  elapsed  since  the  previous  flow.  Thus  this 
event  also  provides  important  information  as  to  the  sto¬ 
rage  volume  of  the  alluvial  fill. 
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Fortunately  event  13  occurred  during  the  presence  in 
the  watershed  of  the  team  doing  rai nfal 1 -i nf i 1 tration- 
runoff  experiments  (see  chapter  9  below).  This  team 
managed  to  organize  itself  and  perform  several  useful 
observations  and  measurements  during  the  event,  as  well 
as  to  make  sure  that  all  instruments  functioned  proper¬ 
ly.  Only  the  vacuum-type  liquid  sampler  did  not  yield 
data,  because  its  intake  level  was  above  the  stream  water 
level  during  the  entire  flow.  Originally  a  very  low  and 
protracted  flood  wave  such  as  in  event  13  was  not  anti¬ 
cipated,  but  steps  have  since  been  taken  to  lower  the 
intake. 

From  manual  water  samples  collected  by  the  team  the 
following  suspended  sediment  concentration  values  were 
obtained : 


Slope  runoff  in  the  vicinity  of  stn  02  1500  ppm 

Channel  flow,  peak  stage  1500 

"  "  ,  immed.  post-peak  600 

"  "  ,  40  min.  after  peak  100 


A  "near-event"  occurred  on  7  February  1979,  and  was 
caused  by  a  low-intensity  rainfall  of  5.4  mm  (station  33) 
spread  over  a  period  of  nearly  3  hours;  4.7  mm  fell 
evenly  during  40  minutes.  As  on  7  February  1978  -  exact¬ 
ly  one  year  earlier  -  only  some  rocky  slopes  generated 
surface  runoff,  and  channel  alluvium  quickly  absorbed 
all  incoming  flow. 

Event  14  came  2  days  later,  on  Friday,  9  February.  A 
small  team  was  working  in  the  watershed  on  routine  main¬ 
tenance  assignments.  Arriving  just  in  time  for  the  effec¬ 
tive  part  of  the  rainstorm,  useful  work  during  the  flow 
could  be  accomplished. 

The  mean  total  rainfall  over  the  watershed  for  event  14 
was  6.9  mm.  Intensities  were  low,  but  the  rainfall  amount 
was  sufficient  to  generate  a  flow  in  the  upper  part  of 
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the  watershed.  Undoubtedly  the  partial  saturation  of 
the  colluvial  and  alluvial  deposits  by  the  rainfall  two 
days  earlier  contributed  to  the  flow.  At  stations  05 
and  04  the  flow  was  very  low,  but  it  lasted  for  about 
1  ^  hours;  its  farthest  downstream  extension  was  only 
200  meters  below  04.  Evidently  the  storage  of  the  main 
alluvial  reach  was  still  in  a  state  of  nearly  total  de¬ 
pletion,  due  to  the  lack  of  any  large  scale  replenish¬ 
ment  during  the  preceding  period  of  four  years  (since 
event  12  -  February  1975). 

The  hydrologic  data  for  event  14  are  listed  in  Table  2.1 

The  presence  of  a  team  in  the  watershed  at  the  time  of 
the  flow  enabled  the  detailed  sampling  of  the  floodwa- 
ters  for  sediment  content.  As  in  a  previous  event,  sta¬ 
tion  05  was  selected  for  the  purpose.  Manual  sampling 
of  the  total  load  was  made  at  2  ^  -minute  intervals, 
starting  with  the  peak  flow.  In  addition,  data  are  avai¬ 
lable  from  the  vacuum-operated  sampler,  whose  level  had 
been  appropriately  lowered  following  event  13,  as  well 
as  one  sample  from  Hayim-type  instrument  77. 

Preliminary  results  indicate  that,  for  total  load,  con¬ 
centration  jumped  from  310  ppm  at  the  time  of  peak  flow 
to  1750  ppm  3  minutes  thereafter.  Concentration  remai¬ 
ned  over  1000  ppm  for  another  10  minutes,  and  fell  gra¬ 
dually  to  110  ppm  in  the  flow  recession.  The  vacuum- 
type  samples  (fine  suspended  material  only)  had  concen¬ 
trations  between  730  and  260  ppm. 

A  detailed  grain-size  analysis  of  all  samples  is  current 
ly  underway. 

Late  May  and  the  first  half  of  June  1979,  had  numerous 
cloudy  days  in  the  project  area.  Situations  of  near¬ 
rain  developed  several  times.  Accordingly,  it  was  deci¬ 
ded  to  prolong  the  recording  season  through  the  end  of 
July.  However,  by  mid-late  June,  weather  returned  to  the 
normal  summer  heat  and  dryness,  and  no  rain  fell  until 
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Table  2.1:  Events  13  and  14,  Nahal  Yael  Watershed 


Hydrologic  Data 


Event  13  -  11  Dec  1978 

Peak  dischrage 
Rai nf al 1 
Rainfall  volume 
Runoff  volume 
Runoff-rainfall  ratio 

Event  14-9  Feb  1979 

Peak  discharge 
Rai nfal 1 
Rainfall  volume 
Runoff  volume 
Runoff-rainfall  ratio 


Stn  05  Stn  04 


0.038 

0.02 

m3/s 

11.0 

11.0 

% 

55  0 

1240 

m3 

53 

no 

m3 

9.6 

9 

% 

0.029 

0.029 

m3/s 

6.9 

6.9 

mm 

345 

770 

m3 

31 

30 

m3 

9 

4 

% 
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the  closing  date  of  this  report. 

Contrary  to  1978,  when  50%  of  the  instruments  were 

kept  in  operation  in  the  watershed  (in  rotation,  to  faci¬ 
litate  summer  maintenance),  it  was  decided  in  1979  to 
dismantle  the  entire  system  for  the  three  summer  months. 

A  few  years  earlier,  it  was  decided  to  continue  one 
half  of  the  basic  instrumentation  throughout  the  summer 
on  a  rotating  basis,  in  order  to  take  care  of  a  possible, 
though  quite  unlikely,  occurrence  of  a  summer  cloudburst 
of  the  1926  type.  This  resulted  in  some  mechanical  pro¬ 
blems  in  certain  instruments,  probably  due  to  excessive 
heat,  which  had  generated  ambient  temperature  levels 
inside  the  instrument  housings  which  were  beyond  the 
tolerance  of  certain  parts.  This  caused  considerable 
problems,  which  will  be  hopefully  prevented  in  future. 

Contrary  to  the  very  rainy  character  of  the  1979/80 
water  year  in  the  entire  country,  including  the  northern 
and  even  central  Negev,  no  runoff  event  occurred  in  Na- 
hal  Yael  during  the  entire  year.  Even  the  destructive 
floods  of  the  southern  Sinai,  which  caused  so  much  da¬ 
mage  in  and  around  Ophira  (Sharm-esh-Shei kh)  left  the 
Nahal  Yael  stream  channels  dry. 
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3.  THE  NAHAL  YAEL  DAM 

2 

A  tentative  water  and  sediment  budget  for  the  0.5  km 
Mahal  Yael  watershed  02  has  been  prepared  on  the  basis 
of  detailed  observations  of  rainfall,  streamflow,  sedi¬ 
ment  transport  and  deposition,  and  of  quatitative  chan¬ 
ges  in  geomorphic  features  (Schick,  1977).  For  the  ten 
year  period  considered ,  (1965/66  to  1974/75),  the  mean 
annual  rainfall  was  31.6  mm.  99  per  cent  of  the  geomor¬ 
phic  work  was  accomplished  during  five  days  by  seven 
discrete  events. 

Sediment  in  suspension,  computed  from  data  obtained  from 
automatic  sampling  program,  accounted  for  a  mean  annual 
yield  of  127  tons.  Bed  load  yield  was  estimated  on  the 
basis  of  (1)  distance  of  transport  determined  from  tra¬ 
ceable  particles,  (2)  determination  of  the  area  and 
depth  of  the  scour  layer  for  the  in-channel  bed  and  for 
the  gravel  bars,  and  (3)  a  comparison  of  grain  sizes 
on  the  bed  and  bars  with  sediment  in  transport  as  sam¬ 
pled  by  liquid  samplers  and  bed  load  traps.  A  mean  annu¬ 
al  yield  of  66  tons  was  found.  The  dissolved  load  was 

only  about  one  per  cent  of  the  total  load.  The  resul- 

2 

ting  mean  annual  sediment  yield  of  388  tons/km  consi¬ 
derably  exceeded  the  accepted  norm  for  arid  environments. 

It  also  exceededby  a  factor  of  3  the  estimated  sediment 

2 

yield,  corrected  for  drainage  area,  for  the  3,100  km 
watershed  of  Wadi  Watir,  located  in  a  similar  environ¬ 
ment  in  Eastern  Sinai. 

While  the  internal  sediment  delivery  ratios  of  the  Na- 
hal  Yael  drainage  system  as  found  were  reasonably  con¬ 
sistent,  the  aggradation  rates,  as  measured  directly  on 
the  alluvial  fan  over  a  period  of  six  years,  were  only 
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one  fifth  of  what  they  should  be  on  the  basis  of  the  se¬ 
diment  transport  computations.  The  main  reasons  for 
this  large  difference  seem  to  be  an  insufficient  under¬ 
standing  of  the  transport  mechanism  during  violent  de¬ 
sert  floods,  exchanges  between  suspended  and  bed  load 
transport  modes,  and  the  importance  of  transient  allu¬ 
vial  storage  and  its  intense  localization.  But  inade¬ 
quacies  of  the  sampling  and  measurement  network  may  also 
be  responsible  for  the  discrepancy,  at  least  in  part. 

The  idea  to  construct  a  100%  trap  efficiency  structure 
for  water  and  sediment  within  the  Nahal  Yael  project  is 
not  new.  Actually  watershed  03  has  been  kept  "free"  for 
this  purpose  since  the  detailed  instrumentation  network 
was  installed  at  the  beginning  of  the  1967/68  water  year. 

The  need  for  calibrating  the  measurements  during  trans¬ 
port  of  water  and  sediment  in  relation  to  uncontroversi - 
al  data  of  total  flow  volume  and  total  sediment  load  be¬ 
came  much  more  pressing  in  the  light  of  the  results  of 
the  ten  year  budget. 

But  watershed  03  did  not  seem  any  longer  to  be  a  suitable  can¬ 
didate  for  calibration  by  damming:it  had  only  one 
gauging  station  with  little  variety  in  its  predominant¬ 
ly  bedrock  channel  (the  alluvial  reach  is  limited  to  the 
lowest  70  meters);  its  suspended  sediment  record  is  dull, 
no  bedload  data  (traps  or  traced  particles)  were  availa¬ 
ble,  and  the  network  of  channel  cross  sections  was  prac¬ 
tically  non-exi stant.  Above  all,  watershed  03  was  too 

2 

small  a  tributary  (drainage  area  0.13  km  )  for  reasonable 
generalizations  of  the  anticipated  results. 

Therefore  the  choice  turned  toward  damming  the  entire  Na¬ 
hal  Yael  watershed,  approximately  at  the  apex  of  the  al¬ 
luvial  fan.  The  decision  to  do  so  was  not  taken  lightly. 

The  construction  meant  the  termination  of  gauging  sta¬ 
tion  01  at  the  foot  of  the  alluvial  fan  near  confluence 
to  Nahal  Roded,  as  well  as  of  the  automatic  sediment  samp- 
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lers  71  and  72  and  of  the  lower  and  middle  bedload  traps. 
Moreover,  it  meant  the  discontinuation  of  the  topogra¬ 
phic  sections  which  monitored  the  changes  on  the  allu¬ 
vial  fan  from  flood  to  flood. 

However,  the  prospect  of  obtaining  an  independent  check 
on  the  runoff  and  sediment  data  obtained  from  the  various 
station  in  the  middle  and  upper  part  of  the  watershed 
outweighed  the  reservations  mentioned  above.  Monitoring 
the  water  flow  -  total  volume  in  the  reservoir  as  well 
as  the  inflow  (and  overflow)  discharge  -  would  conside¬ 
rable  improve  the  rating  curves  of  the  gauging  stations. 
Accurate  volumetric  determinations  of  the  sediment  vo¬ 
lume  which  accumulates  in  the  reservoir  during  discrete 
events  would  lend  much  confidence  to  the  sediment  dis¬ 
charge  figures  which  were  up  to  now  based  mostly  on 
floodwater  samples  abstracted  during  the  rise  of  the  flow. 

Additional  benefits  of  such  a  venture  would  be:  I)  ob¬ 
servation  on  and  accumulation  of  experience  with  regard 
to  the  construction  operational  aspects  of  small,  simple 
and  inexpensive  earth  dams  in  an  extremely  arid  envi¬ 
ronment;  and  II)  evaluation  of  longer-term  environmen¬ 
tal  effects,  caused  by  the  expected  increased  percola¬ 
tion  in  the  reservoir  bed,  with  subsequent  vegetation 
response  (mainly  acacia  trees)  and  the  concomitant  in¬ 
crease  of  faunal  activity  (mainly  gazelles)  -  all  fac¬ 
tors  which  could  result  in  changing  the  spot  into  a  plea¬ 
sant  rest  site  for  hikers  and  motorists. 

Detailed  plans  were  accordingly  drawn  in  1977  for  the 
construction  of  a  7  meter  high  earth  dam,  with  an  im- 

3 

pounding  capacity  of  about  8,500  m  .  A  flood  volume  of 
such  a  magnitude  is  caused  by  an  event  frequency  of  one 
in  15  to  20  years.  The  amount  of  gravel  and  sand  invol- 
ved  in  such  a  structure  was  estimated  at  6,000m  ,  to  be 
derived  locally  from  the  present  and  Pleistocene  fans 
immediately  downstream  the  dam  s i te .  T he  pi ans  also  called 
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for  the  reduction  of  the  dam  permeability,  conditioned 
by  the  coarse-grained  character  of  the  material  on  si¬ 
te,  by  suitable  methods  such  as  a  plastic  coating. 

A  concrete  spillway  with  an  appropriate  hydraulic  con¬ 
trol  were  designed  for  the  right-hand  valley  side  near 
rainfall  recorder  30;  in  such  a  way  the  potentially 
dangerous  overflow  would  be  kept  away  from  the  earth 
embankment . 

Because  of  financial  constraints,  the  project  as  des¬ 
cribed  above  had  to  be  subsequently  redesigned  on  a 

more  modest  scale.  The  structure  as  actually  built 

3 

used  only  3,000  m  of  earth.  It  has  an  overall  height 
of  5.5  meters,  a  maximum  depth  of  impoundment  of  5  m, 

3 

and  a  water/sediment  capacity  of  somewhat  over  3,000  m  . 
The  material  was  local  sand  and  gravel,  with  large 
boulders  forming  the  base  of  the  downstream  embankment 
slope.  Construction  started  in  December  1977  and  was 
completed  in  summer  1973.  (Photographs  1-3).  Ancillary 
work  completed  a  few  months  thereafter  included: 

I)  A  concrete-lined  overflow  channel  with  a  prefabri¬ 
cated  metal  crump  weir  were  constructed  and  installed 
on  the  right-hand  valley  side  over  a  bedrock  sill  se¬ 
veral  meters  away  from  the  end  of  the  embankment. 

II)  Two  independent  water  level  recorders,  to  monitor 
the  impounded  water  level  and  the  depth  of  overflow. 

One  is  a  conventional  float-operated  Ott  recorder  and 
the  other  is  a  pneumatic  type  Rimco  (Australia)  recor¬ 
der  . 

III)  The  reservoir  bottom  was  completely  graded,  in  order 
to  facilitate  post-flood  accurate  surveys  of  water  and 
sediment.  Several  marker  plates,  made  of  sheet  metal, 
were  installed.  These  should  facilitate  augering  and 
sampling  after  the  expected  silting.  The  plates  were 
coated  with  a  layer  of  sediment  particles,  to  keep  their 
hydraulic  roughness  similar  to  that  of  the  "natural" 
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bottom.  Approximately  25  benchmark  referenced  cross- 
sections  were  accurately  surveyed  across  the  expected 
lake  area. 

IV)  The  presence  of  the  earthmoving  equipment  on  the 
site  was  utilized  for  an  exploratory  excavation,  2.5m 
deep,  at  the  apex  of  the  active  alluvial  fan.  The  ex¬ 
posed  section  of  alluvium  is  being  studied  in  relation 
to  a  similar  section  in  the  Pleistocene  fan,  also  fresh¬ 
ly  exposed  during  the  construction.  The  purpose  of  this 
study  is  to  evaluate  the  spectrum  of  hydrologic  events 
in  each  of  the  fans,  with  a  bearing  on  extreme  events 
recorded  and  on  climatically  effected  changes. 

V)  The  steep  downstream  side  of  the  overflow  channel 
has  been  protected  with  combined  stone-concrete  masonry 
work,  to  prevent  erosion  of  the  toe  of  the  embankment 
by  a  major  overflow.  The  dam  was  designed  in  such  a  way 
as  to  lose  water  by  controlled  seepage  through  the  em¬ 
bankment.  Although  the  rate  of  seepage  is  as  yet  unknown, 
it  is  expected  to  be  slow  enough  so  that  the  interpreta¬ 
tion  of  the  inflow  rate  into  the  lake  and  the  eventual 
outflow  rate  through  the  overflow  channel  will  not  be 
seriously  jeopardized.  For  future  flows  which  will  fill 
the  reservoir,  the  seepage  losses  are  expected  to  gradu¬ 
ally  decl i ne . 

During  the  first  two  water  years  after  the  completion  of 
the  dam  (1978/79  and  1979/80)  not  a  single  drop  of  water 
passed  through  Nahal  Yael  station  02  to  enter  the  reser¬ 
voir.  During  the  thirteen  previous  years  a  flood  which 
would  have  flowed  into  the  reservoir  occurred  eight  times. 
In  fact,  no  flow  was  recorded  at  station  02  (at  the  top  of 
the  dry  waterfall  immediately  upstream  of  the  reservoir) 
since  event  12  (20/21  February  1975),  i.e.  during  the 
last  five  and  a  half  years. 


3  A  reservoir  sedimentation  reference  plate 
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It  is  hoped  that  the  dam  will  yield  important  informa 
tion  for  research  as  well  as  of  an  applied  nature. 
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4.  THE  RELATIONSHIP  BETWEEN 
RAINFALL  AND  EROSION  IN  NAHAL  YAEL 


by  Judith  Lekach  and  Asher  P. Schick 


In  order  to  better  understand  some  of  the  factors 
which  underlie  the  fluvial  erosion  process  of  the  ex¬ 
treme  desert,  as  represented  by  the  by  now  numerous 
data  obtained  from  the  network  of  multiple  stage  auto¬ 
matic  suspended  sediment  samplers  (Schick  &  Sharon, 

1974,  and  subsequent  data),  several  preliminary  analy¬ 
ses  of  the  Nahal  Yael  record  were  undertaken. 

A  fairly  large  number  of  two-variable  and  error  regres¬ 
sions  between  various  assumed  causative  factors  (most¬ 
ly  rainfall  derivatives)  and  the  erosion  data  -  repre¬ 
sented  by  the  mean  vertical  suspended  sediment  concen¬ 
tration  -  were  conducted.  For  small  non-alluvial  water¬ 
sheds  such  as  03,  04,  and  05  at  Nahal  Yael,  the  best 
fit  was  consistently  obtained  between  a  rainstorm  deri¬ 
vative  defined  as  effective  rainfall  energy  calculated 
for  3-minute  time  intervals,  and  the  mean  vertical  sus¬ 
pended  sediment  concentration.  The  linear  regressions 
yielded  correlation  coefficients  between  0.94  and  0.83. 
Then  the  goodness  of  fit  was  further  improved  by  varying 
I)  the  initial  rainfall  loss  prior  to  the  initiation  of 
runoff,  and  II)  the  infiltration  rate  during  the  period 
of  effective  rainfall  (assumed  to  be  constant).  By  best- 
fit  procedures,  the  initial  loss  was  found  to  be  a 
function  of  drainage  area,  as  follows: 


Fig.  4.1:  Fluvial  suspended  sediment  content  and  rainfall  ener¬ 
gy  for  floods  in  extremely  arid  small,  non-alluvial 
watersheds,  based  on  Nahal  Yael  watersheds  03,  04,  and 
05.  The  lines  inside  the  diagram  represent  drainage 
area  in  square  kilometers. 


watershed  05  -  4  mm 
watersheds  04  &  03  -  5.5  mm 

The  procedure  further  indicated  a  preferred  constant 
infiltration  rate  of  0.5  mm/3  min.  Both  findings  corro¬ 
borate,  independently,  the  general  range  of  the  conclu¬ 
sions  concerning  runoff  generation  summarized  in  Schick 
(1970),  and  are  also  in  general  agreement  with  results 
of  infiltration  tests  summarized  in  this  report  (see 
chapter  9  below) . 

The  resulting  erosion  rate  model  is  summarized  in  fi¬ 
gure  4.1,  which  correlates  the  effective  rainfall  in¬ 
tensity  of  the  event  with  the  mean  vertical  suspended 
sediment  concentration,  for  various  watershed  sizes.  A 
similar,  though  not  identical,  relationship  was  found 

satisfactory  for  erosion  values  derived  from  small 
2 

(2m  )  plots  (Pearce  1976),  indicating  the  validity  of 

the  Hortonian  overland  flow  model  in  preference  to  the 

partial  area  contribution  model.  The  fact  that  physio- 

graphically  fully  developed  watersheds  in  the  0.05-0.5 
2 

km  size  range  behave  identi cal ly, further  supports  the 
Hortonian  overland  flow  interpretation  over  the  satura¬ 
tion  overland  flow  concept  in  arid  terrain. 

The  relationship  shown  on  fig.  4.1  breaks  down  when  da¬ 
ta  from  watersheds  which  include  a  sizeable  portion  of 
continuous  channel  alluvium  are  used.  Under  these  con¬ 
ditions,  which  are  practically  universal  for  desert  wa- 
tersheds  whose  area  exceeds  0.3  km  ,  an  additional  fac¬ 
tor  must  be  introduced. 


5.  PHYSIOGRAPHIC  CHARACTERISTICS 
OF  AREAL  UNITS 

by  Oded  Salmon,  Eyal  Taubenhaus  and  Asher  P. Schick 

5.1  Selected  elements 

In  order  to  obtain  a  "feel"  for  data  collection,  hand¬ 
ling,  and  relevance  for  the  purposes  of  developing  a 
deterministic  and  distributed  hydrologic  and  geomor- 
phic  model  for  Nahal  Yael ,  six  areal  elements  in  water¬ 
shed  05  were  selected  for  detailed  field  study.  This 
work  was  carried  out  at  a  time  when  our  idea  was  to  de¬ 
velop  the  model  on  the  basis  of  a  triangular  grid.  In 
such  a  framework,  each  areal  element  is  an  equilateral 

triangle,  side  21.65  m,  height  18.75  m,  and  unit  area 

?  2 
very  close  to  200  m  (actually  202.975  m  ). 

This  triangular  grid  was  subsequently  abandoned  as  the 
basis  for  the  model  in  favor  of  a  rectangular  grid.  The 
information  assembled,  however,  continues  to  be  of  im¬ 
portance  as  a  basis  for  correlations  between  field 
truths  and  information  derived  manually  or  automatically 
by  means  of  remote  sensing  methods. 

Each  of  the  six  sites  was  first  located  on  an  aerial  pho¬ 
tograph  by  using  criteria  of  terrain  representativity. 
Then  each  triangle  was  sited  in  the  field.  The  corners 
were  staked  and  the  sides  marked  and  bounded  by  colored 
plastic  tape.  Each  triangle  was  photographed  in  totali¬ 
ty  from  behind  of  each  of  its  three  apices.  Several  ter¬ 
rain  closeups  were  also  taken. 
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A  level  survey  determined  the  relative  elevation  of  the 
three  apices  and  of  the  centermass  point  of  the  tri¬ 
angle.  This  procedure  was  aimed  specifically  at  getting 
a  measure  of  "internal  relief"  -  a  factor  of  potential 
relevance  to  the  approximation  of  several  physiographic 
controls  for  the  model,  as  well  as  one  possible  objec¬ 
tive  yardstick  for  a  further  subsequent  sub-division  of 
certain  areal  elements  until  a  desired  degree  of  homo¬ 
geneity  is  achieved  (see  chapter  8). 

A  field  sketch  of  each  triangle  was  made,  showing  the 
distribution  of  colluvium  and  of  rock  (defining  diffe¬ 
rent  types),  and  topographic  irregularities  such  as 
scarplets,  joints,  rills,  hollows,  etc.  Information  on 
colluvium  character  and  particle  size  was  also  noted. 

In  addition,  the  following  measurements  or  semi-quan¬ 
titative  estimates  were  recorded  for  each  triangle: 

(1)  percent  bedrock  and  percent  colluvium; 

(2)  depth  of  colluvium; 

(3)  density  of  colluvium  (qualitative  estimate  in  three 
classes); 

(4)  roughness  of  colluvium  surface  (four  classes); 

(5)  median  size  of  colluvium  particles; 

(6)  concentration  potential  of  runoff  (five  classes) 
judged  mainly  by  connectivity  of  bedrock  patches; 

(7)  "internal  relief",  defined  as  the  difference  in 
elevation  between  the  actual  height  of  the  center- 
mass  point  and  the  mean  value  of  the  heights  of  the 
three  apices. 

In  most  parameters  significant  differences  among  tri¬ 
angles  was  obtained.  The  prospects  of  obtaining  most  of 
values  from  remote  sensing  procedures  seem  good.  As  a 
first  step,  a  relationship  will  be  established  between 


|o  „°°J  colluvium 
gg  bedrock  &  colkjvium 
titHi  bedrock 
HU  dyke 


Fig.  5.1:  Terrain  lithology,  Nahal  Yael  watershed  05.  Based  on 
aerial  photograph  interpretation  with  limited  ground 
control.  Map  has  not  been  corrected  for  photo  distortion 
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the  detailed  site  character i zati on  and  the  mapping  off 
aerial  photographs  of  bedrock/coll uvi urn  areas  (para¬ 
graph  5.2  below). 

5.2  Bedrock/  colluvium  ai rphotomap ,  watershed  05 

Figure  5.1.  shows  the  distribution,  in  Nahal  Yael  wa¬ 
tershed  05,  of  bedrock,  colluvium,  and  col luvium/bed- 
rock  combination.  Major  dyke  rocks  are  also  shown. 

The  "colluvium/bedrock  combination”  represents  an  inter¬ 
fingering  of  areas  of  the  two  lithologies  too  fine  to 
differentiate  at  the  scale  employed.  In  many  cases  it  re¬ 
presents  intermittent  and/or  thin  patches  of  colluvium 
which  overlie  bedrock  which  crops  out  at  fairly  close 
intervals  beneath  the  colluvium. 

The  distribution  was  mapped  on  the  basis  of  aerial  pho¬ 
tography  interpretation,  with  only  a  limited  amount  of 
ground  control . 

A  rough  approximation  of  the  areal  extent  of  the  three 
main  formations  is: 

bedrock  30  per  cent 

colluvium  10  "  " 

colluvium/bedrock  60  "  " 


The  colluvium/bedrock  combination  is  estimated  to  com¬ 
prise,  areally,  two-thirds  colluvium  and  one-third  bed¬ 
rock.  Substituting  these  figures,  we  find  that  roughly 
one  half  of  watershed  05  is  covered  by  a  colluvial  ve¬ 
neer;  the  other  half  is  bare  rock. 

The  areal  distribution  (fig.  5.1)  will  be  used  in  fu¬ 
ture,  together  with  the  geological  map,  the  detailed 
triangles  study,  and  additional  work,  to  define  numeri¬ 
cally  each  unit  area  for  physical  characteristics  such 
as  infiltration  and  roughness. 
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6.  BEDLOAD  TRANSPORT  AND  STREAM  POWER 
EXAMPLES  FROM  DESERT  FLOODS 


by  Judith  Lekach  and  Asher  P.  Schick 


6.1  Bedload  transport  and  stream  power  for  two  floods, 
Nahal  Yael  watershed  05. 

Events  12A  (20  February  1975)  and  14  (9  February  1979) 
occurred  in  Nahal  Yael,  fortuitously,  during  field  pe¬ 
riods  of  the  study  team,  whose  presence  on  site  facili¬ 
tated  the  collection  of  usually  unobtainable  manually 
sampled  and  time-discharge  controlled  sediment  data. 

Both  flows  were  really  of  a  relatively  small  magnitude, 

3 

peak  discharges  between  0.04  and  0.12  m  /s.  But  given 
the  extreme  paucity  of  such  data,  even  these  flows 
should  provide  an  important  insight  into  the  factors 
and  mechanisms  of  sediment  transport  in  small  extreme¬ 
ly  arid  watersheds. 

Event  12A  was  rather  untypical:  38.4  mm  of  low-inten¬ 
sity  rain  (mean  between  valley  station  30  and  ridge 
station  26)  fell  over  a  period  of  9  ^  hours,  and  ge¬ 
nerated  during  its  latter  half  an  extended  hydrograph 
with  several  low  peaks.  During  the  flow  APS  collected 
manually  25  samples  of  the  total  load  at  the  foot  of  a 
bedrock  riffle  immediately  upstream  of  gauging  station 
05.  Most  of  the  samples  were  taken  from  the  first  flood- 
wave  and  from  its  recession,  but  a  few  were  also  taken 
at  certain  intervals  from  the  subsequent  floodwaves. 


unit  bedload  transport  rate  ib  .  kg /ms 


1 


6.1 


unit  stream  power  « ,  kg/ms 


Fig.  6.1:  Bedload  transport  and  stream  power,  Nahal  Yael  watershed 
05.  Based  on  events  12A  (20  Feb  1975)  and  14  (9  Feb  1979). 
Figures  inside  graph  indicate  median  grain  size  of  the 
bedload  (>  0.063  mm)  only. 


V 


Event  14  was  a  small  flood,  very  typical  of  desert  con¬ 
ditions:  a  short  high- intons i ty  rain  yield  a  steep  flow 
rise  with  a  slower  declina,  uncomplicated  by  any  addi¬ 
tional  rai nfal 1 . 

Sediment  concentration  peaked  simultaneously  with  dis¬ 
charge  in  event  12A,  and  was  slightly  later  than  peak  dis¬ 
charge  in  event  14.  There  is  some  suggestion  of  a  secon¬ 
dary  peak  in  concentration  down  the  recession  of  event 
12A,  but  this  needs  further  corroboration. 

The  samples  of  both  floods  were  subjected  to  a  particle 
size  analysis.  For  the  following  analysis  only  the  bed¬ 
load  fraction,  defined  here  as  larger  in  size  than  0.063mm, 
is  considered.  The  median  size  of  the  bedload  fraction  is 
determined,  and  the  associated  unit  bedload  transport  ra¬ 
te  and  unit  stream  power  are  computed.  The  unit  power  nee¬ 
ded  for  initial  motion  proved  to  be  negligible,  and  is 
disregarded.  Table  6.1  shows  the  detailed  data  for  event  14. 

The  data  for  both  events  align  well  on  the  standard  ib  ver¬ 
sus  (o  graph,  and  seem  to  parallel,  more  or  less,  data  for 
other  streams  (see  6.2  below).  A  comprehensive  treatment 
of  the  bedload  transport  rate  in  Nahal  Yael  must  incorporate 
data  on  larger  particles  such  as  available  from  bedload 
traps  and  tracing  experiments. 

6.2.  Bedload  transport  and  stream  power  for  an  extreme 
event,  Wadi  Mikeimin,  Sinai 

§i?il__Ab§tract 

A  reconstruction  of  a  low-frequency  high-magnitude  flood 
event,  which  occurred  in  1971  in  Wadi  Mikeimin  in  the  arid 
SE  Sinai  Mts.,  enabled  the  evaluation  of  bedload  transport 
rates  in  relation  to  flow  characteristics.  Unit  transport 
rates  between  20  and  100  kg/m.s  (at  least  one  order  of  mag¬ 
nitude  higher  than  nearly  all  data  previously  reported), 
accord  with  data  on  other  rivers.  The  proportionality  of 
bedload  transport  rate  to  unit  stream  power  in  excess  of 
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that  necessary  for  initial  moti on, rai sed  to  the  power  of 
3/2,  has  been  validated  for  these  very  high  transport  ra¬ 
tes  which  approach  debris  flows  in  character. 


6i2^2__Introductign 

Direct  measurement  of  the  transport  of  coarse  bedload 
in  heavily  laden  desert  floods  is  practically  impossible 
to  accomplish.  In  certain  fortuitous  cases,  in  which  the 
total  load  is  found  as  a  deposit  and  measurable,  and 
enough  data  are  available  to  reconstruct  the  stream- 
flow  and  evaluate  its  characteristics,  indirect  estima¬ 
tes  of  bedload  transport  are  possible.  The  value  of  the¬ 
se  estimates  is  two-fold:  They  may  offer  additional  ex¬ 
perimental  validation  of  the  proportionality  between 
bedload  transport  rate  and  the  3/2  power  of  excess  unit 
stream  power  (1,2);  and  they  extend  the  body  of  data  be¬ 
yond  the  range  presently  available  from  direct  measure¬ 
ments  (3,  4).  Such  estimates  for  an  exceptional  event 
in  the  Jordan  and  Heshushim  rivers,  northern  Israel, 
yielded  compatible  data  of  bedload  transport  rate  and 
stream  power  one  order  of  magnitude  higher  than  previous¬ 
ly  known  (5).  The  exceptional  event  of  January  1971  in 
Wadi  Mikeimin,  southeastern  Sinai,  which  caused  large 
quantities  of  bedload  to  be  transported  and  deposited 
them  in  entirety  as  a  new  alluvial  fan  at  the  watershed 
outlet,  enables  another  controlled  estimate  of  the  same 
type. 

Wadi  Mikeimin  is  an  uninhabited  tributary  of  Wadi  Watir, 

2 

which  drains  3100  km  of  extremely  arid  mountainous 

country  towards  the  large  alluvial  fan  of  Neviot  (Nou- 

2 

eiba)  on  the  shore  of  the  Gulf  of  Aqaba.  The  12.9  km 
watershed  is  entirely  composed  of  precambrian  basement 
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rocks  with  tertiary  volcanic  intrusions,  crossed  by  nu¬ 
merous  NNW-SSE  runoff  faults.  The  mean  annual  rainfall  is 
estimated  at  30-40  mm  and  is  highly  erratic.  Due  to  a 
relief  difference  of  600  m  in  the  5  km  long  watershed, 
slopes  are  rocky  and  steep  but  covered  in  part  by  coar¬ 
se  grained  colluvium.  Average  stream  channel  gradients 
of  Wadi  Mikeimin  and  of  its  main  tributaries  are  bet¬ 
ween  5  and  10  per  cent.  Most  of  the  valley  network  is 
alluvium-covered,  braided  between  rocky  banks,  and  its 
low  bars  are  composed  of  particles  mean  size  between 
29  and  53  mm.  Large  boulders  of  diameter  up  to  1500  mm 
are  also  present. 

6^2^4 _ Mikeirnln_Flgod 

In  January  1971  the  sudden  appearance  of  an  alluvial 
barrier  at  the  mouth  of  Wadi  Mikeimin  was  reported. 
Being  located  at  the  western  margins  of  the  Ein  Fortagha 
springs  area,  the  barrier  caused  the  emergence  of  an  up 
to  400  m  long  "finger"  lake  on  its  upstream  side. 

The  barrier  proved  to  be  an  alluvial  fan  deposited  by 
a  highly  localized  and  short-lived  flood  event  in  Wadi 
Mikeimin.  Detailed  field  work  both  before  and  after  the 
complete  erosion  of  the  fan  by  the  November  1972  flood 
in  Wadi  Watir  (7)  enabled  a  reconstruction  of  the  water 
and  sediment  characteristics  of  the  Mikeimin  event  (8). 

6^2^5 _ Methods 

The  analysis  involved  in  the  reconstruction  of  the  event 
proceeded  along  the  following  stages:  (a)  high  water 
marks  at  8  judiciously  selected  stations  in  the  water¬ 
shed  enabled  the  evaluation  of  peak  discharges  for  the 
event;  (b)  comparison  with  other  flood  events  which 
occurred  over  the  last  15  years  in  Nahal  Yael  and  in 
Eastern  Sinai  (9,  10),  evaluated  in  the  light  of  vari¬ 
ous  studies  on  rain  characteristics  (11)  enabled  the 
construction  of  probable  hypothetical  hydrographs. 


Table  6.2.  Summary  reconstructed  data  of  river  hydrau¬ 
lics  and  bedload  transport.  Wadi  Mikeimin, 
Sinai 


Station 

M, 

main  channel 

tributary  D 

tributary  ! 

Discharge  range 

5  -  68 

5-80 

2-47 

Q,  m3/s 

Width 

w,  m 

18.8 

28.5 

40 

Slope 

s,  m/m 

0.087 

0.085 

0.075 

Depth 

d,  m 

0.4 

0.36 

0.20 

a 

Unit  stream  power 
u,  kg/m.s 

169 

127 

46 

Bedload  size*5 
°50,  mm 

31C 

16d 

8d 

Bedload  yield 
tons® 

10,000f 

8,000g 

3.4009 

Duration^1 

min 

60 

50 

35 

k 

Unit  transport  rate 

1b,  kg/m.s 

92 

58 

25 

Unit  power1  at  initial 
motion 
w0»  kg/m.s 

0.37-3.6 

0.15-1.5 

0.05-0.52 

Excess  unit  power11 
o>-(i>o»  kg/m.s 

165 

125.5 

46 

Ratio  of  flow  depth  to 

median  bdld  size  Y/Dc_, 

ou 

m/m 

12.9 

22.5 

25 

Unit  transport  raten 

45 

62 

21 
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a  to  =  Qs/w;  mean  value  for  discharge  range. 

k  Estimated  median  diameter  of  particles  (b  axis). 

c  Based  on  measurements  along  the  truncated  face  of  the 
Mi keimi n  al 1 uv i al  fan. 

d  Based  on  Leopold  (13)  counts  at  the  station. 
e  Immersed  weight,  metric  tons. 

^  Estimated  from  deposit. 

9  Estimated  by  the  Einstein  (12)  procedure,  compared 
with  the  deposit. 

L. 

Estimated  from  reconstructed  hydrographs. 

i/ 

Based  on  this  research. 

^  Estimated  by  multiplying  the  Shields  critical  shear 
stress  by  the  reconstructed  mean  velocity;  range  re¬ 
flects  two  different  variants  of  the  computation. 

m  Using  highest  value  of  estimated  Uo . 

n  Computed  by  the  Bagnold  (2,  equation  4)  procedure. 
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(c)  the  Einstein  bedload  computation  procedure  (12)  was 
applied  to  three  of  the  cross  sections,  in  relation  to 
the  computed  hydrographs  and  bed  material  size.  The 
bulk  size  distribution  of  the  alluvial  fan  deposit  was 
also  taken  into  account,  (d)  Comparison  between  the 
total  bedload  yield  as  computed  with  the  actual  bed¬ 
load  volume  as  deposited  in  the  Mikeimin  alluvial  fan 
enabled  the  validation  of  the  components  of  the  proce¬ 
dure:  i.e.,  because  of  the  agreement  between  the  compu¬ 
ted  and  observed  bedload  volumes,  the  data  on  which  the 
bedload  computation  is  based  (bedload  transport  rate, 
hydrographs)  are  assumed  to  be  reasonably  correct. 

§i?i§__B§§ylts_and_Discussion 

Unit  transport  rates  in  the  range  of  20-100  kg/m.s  are 
indicated  (Table  6.2)  as  a  mean  value  for  a  30-60  min 
long  flow  event  in  a  small  arid  watershed  with  very 
steep  stream  channels.  These  rates  derive  from  the 
Bagnold  (2)  procedure.  The  data  suggest  that  the  propor¬ 
tionality  of  bedload  transport  rate  to  unit  stream  power 
in  excess  of  that  necessary  for  initial  motion,  raised 
to  the  power  of  3/2,  is  valid  also  for  short-lived  vio¬ 
lent  torrents  which  approach  debris  flows  in  character. 

The  data  for  Wadi  Mikeimin  agree  well  with  those  for 
Elbow  River  (14)  ( fi g .  6.2) ; both  rivers  transport  bed¬ 
load  of  similar  median  size. 

The  Mikeimin  data  also  agree  well,  despite  their  very 
special  characteristics,  with  data  from  other  rivers  re¬ 
presenting  widely  varying  conditions  (2,3,4,5,15).  Sub¬ 
ject  to  the  morphological  distribution  between  rivers 
whose  bed  material  grain  size  approximates  the  size  of  the 
transported  bedload  and  those  whose  effective  bed  surface 
roughness  size  is  significantly  larger  (2),  and  provided 
bedload  material  is  available  in  unlimited  supply,  these 
new  data  provide  further  encouragement  in  the  search  for 


w  -  w0, excess  unit  power,  kg/ms 


Bedload  transport  as  a  function  of  excess  stream  power 
in  relation  to  data  from  Wadi  Mikeimin,  Eastern  Sinai. 
Bedload  transport  rate  is  by  immersed  weight.  Figures 
inside  graph  indicate  representative  median  grain  size 
in  meters. 
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a  satisfactory  universal  method  for  predicting  bedload 
transport  rates  of  rivers  on  the  simple  basis  of  excess 
power,  size  and  depth  alone. 
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7 .  GEOPHYSICAL  SURVEY  OF  NAHAL  YAEL  CHANNEL  BED 
by  Asher  P.  Schick 

In  sequence  to  a  few  seismic  soundings  carried  out  by 
Bull  and  Schick  four  years  ago,  a  more  systematic  sur¬ 
vey,  consisting  of  about  forty  seismic'  soundings,  was 
carried  out  by  means  of  a  portable  refraction  instru¬ 
ment. 

The  objective  of  the  survey  was  to  obtain  a  reasonably 
accurate  quantitative  estimate  of  the  potential  water 
storage  in  the  alluvium  beneath  the  channel  bed,  so  that 
this  factor  can  be  evaluated,  together  with  other  water 
losses,  which  affect  the  generation  of  floods  within 
the  watershed. 

Since  the  main  alluvial  body  of  the  channel  system  of 
the  catchment  above  the  waterfall  (station  02)lies  in 
the  one  kilometer  long  main  channel  between  stations  02 
and  04,  the  sounding  activity  was  confined  to  this  reach. 

Excepting  the  central  part  of  the  alluvial  fan  beyond  the 
mountain  front,  where  bedrock  depths  in  excess  of  8  me¬ 
ters  were  indicated,  all  data  within  the  alluvial  chan¬ 
nel  range  between  0.5  and  2.0  meters.  In  general,  for  the 
main  alluvial  reach  between  the  beginning  of  continuous 
alluvial  fill  (station  02)  and  the  waterfall  (station  04), 
the  following  picture  emerges  (Table  7.1):  for  the  upper 
and  lower  sub-reaches,  a  fill  of  around  0.6  to  1.2  meters 
is  indicated,  while  for  the  middle  sub-reach  the  figure 
is  more  like  1.5  meters.  This  is  in  line  with  the  fact  that 
the  alluvial  terraces  in  the  middle  and  middle-upper  sub¬ 
reaches  are  also  considerably  higher  above  the  present 
channel  bed  than  those  at  the  two  ends  of  the  main  allu¬ 
vial  reach.  This  similarity  of  behavior  between  the 
highest  and  deepest  flow  profiles  is  remarkable  in  that 
it  suggests  a  unified  geomorphic  mechanism  for  periods 
of  cut  and  of  fill,  such  as  proposed,  for  example, by  the 
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superflood  theory. 

Three  years  ago  Bull  and  Schick  dug  a  few  pits  in  an  at¬ 
tempt  to  obtain  a  correlation  between  their  seismic 
soundings  and  the  actual  section  exposed.  In  one  of  tho¬ 
se  pits,  near  470,  bedrock  was  not  encountered  even  at 
the  depth  of  3.30  meters;  at  this  point  the  pit  had  to 
be  abandoned,  as  it  was  completely  silted  up  by  event 
12  which  occurred  on  the  next  morning.  Because  of  the 
fact  that  in  the  present  survey  no  sounding  showed  a 
depth  more  than  2  meters,  it  could  well  be  that  the  pit 
dug  by  Bull  and  Schick  was  at  a  chance  and  unrepresenta¬ 
tive  site.  But  it  could  also  mean  that  very  close  to 
concave  bends  undercutting  an  alluvial  terrace  -  as  was  the 
case  with  that  pit  -  deeper  fills  are  present.  Usually 
it  is  impossible  to  align  a  proper  sounding  line  at  such 
sites  because  of  topographic  difficulties. 

The  dominant  width  of  the  "probable"  flood  channel  varies 
systematically  from  about  3.5  m  at  station  400  (60  meters 
downstream  of  station  04),  to  about  24.5  m  at  station 
485  (70  meters  upstream  of  station  02)  (see  fig.  9,  cur¬ 
ve  2,  in  Schick  1977).  Using  these  figures  in  conjunc¬ 
tion  with  the  data  of  table  10.1,  we  arrive  at  the  fol¬ 
lowing  volumetric  estimate  of  the  sub-channel  alluvial 
fill  of  Nahal  Yael  above  the  waterfall: 


Length 

(m) 

Width 

(m) 

Depth 

(m) 

Volume 

(m3) 

Upper  part,  400-430 

300 

7.4 

1.12 

2,486 

Middle  part,  430-477 

470 

18.4 

1.21 

10,464 

Lower  part,  477-490 

130 

23 

1.04 

3,110 

Total,  400-490 

16,060 

Clearly  this  figure 

of  total 

,  sub 

-channel 

level  alluvi- 

al  fill  volume  of  the  Nahal 

Yael 

watershed 

is  only  some 

what  more  accurate  than  an  order  of  magnitude.  In  its 
computation  we  have  neglected  a  number  of  factors,  such 
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Table  7.1:  Summary  of  seismic  soundings,  Nahal  Yael 
alluvial  channel 

Distance  marker*  Depth  to  bedrock**  Remarks 

lO^n  m 


510-508 

13,  9,  8.0? 

on  alluvial  fan  in 

505-502 

1.76 

damned  area 

- — - 

-  Station  02 

488 

0.8  0.55 

486 

1.8  1.1 

480 

0.9  1.1 

-  Station  03 


474 

0.6  0.8 

472 

1.1 

470 

1.0 

468 

0.7 

464 

1.95 

462 

1.3 

458 

1.6? 

454 

1.0 

450 

1.05 

446 

1.0 

442 

0.8 

440 

1.3? 

438 

1.8 

436 

1.1 

434 

1.3 

432 

1.6 

—  Station  75 


428  1.4 

0.9?  on  terrace  (conglome- 


■V  W  IH-l.-jMi* 
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Table  7.1  (cont.) 

416  1.4 

414  1.2 

410  0.6  1.7  1.8  1.2  1.0 

398  upstream  end  of 

alluvial  channel; 
Station  04 


Distance  alongside  main  channel  of  Nahal  Yael ,  in  tens 
of  meters.  Note  that  upper  divide  is  not  zero.  System 
originates  at  a  point  located  between'  the  dam  and  sta¬ 
tion  02,  labelled  500  in  order  to  avoid  negative  values. 


Result  rounded  off  to  nearest  decimeter.  Due  to  the  si¬ 
milarity  of  sound  velocity  through  alluvium  and  through 
intensely  weathered  granite  and  schist,  the  depth  values 
may  include  a  bottom  layer  of  weathered  bedrock  1-2  de¬ 
cimeters  thick. 
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as  the  usual  curved  shape  of  a  fluvially  formed  cross 
section  instead  of  a  rectangular  form,  or  the  disregard 
of  lateral  effluence  during  floods  at  alluvial  banks. 
Incidentally,  these  two  factors  may  quantitatively  off¬ 
set  one  another. 

Assuming  that  the  volume  of  alluvial  channel  fill  of 
3 

16000  m  is  correct,  and  assuming  further  a  mean  poro¬ 
sity  of  25-30  per  cent  of  the  alluvial  material,  we  ar- 

3 

rive  at  a  total  potential  water  storage  of  4000-5000  m 
This  represents  equivalent  watershed  rainfall  of  8  - 
10  mm  -  a  reasonable  figure  in  agreement  with  the  avai¬ 
lable  runoff  data . 

The  question  whether  the  entire  sub-channel  alluvial  bo 
dy  must  indeed  be  saturated  before  runoff  occurs  at 
station  02  cannot  as  yet  be  conclusively  answered. 

Other  factors  as  yet  unknown  are  the  amount  of  shallow 
groundwater  stored  near  the  alluvium-bedrock  contact 
(preferably  in  pools,  perhaps  at  the  foot  of  sub-allu¬ 
vial  waterfalls?),  as  well  as  the  rate  of  depletion  of 
this  storage  with  increasing  time  interval  between  suc¬ 
cessive  runoff  events. 


8.  INFILTRATION  TESTS 


by  Oded  Salmon  and  Asher  P. Schick 


8.1  Introduction 

The  work  reported  here  deals  with  the  determination  of 
approximate  quantitative  values  of  infiltration  rates 
for  several  lithologies  in  an  extremely  arid  environ¬ 
ment.  These  values  were  arrived  at  by  means  of  infil¬ 
tration  trials  which  utilized  simulated  rainfall  on  very 
small  plots,  supplemented  by  comparisons  with  other  da¬ 
ta  whenever  available. 

The  impetus  for  the  work  was  two-fold.  First,  the  in¬ 
filtration  data  are  of  prime  importance  for  understan¬ 
ding  the  mechanism  of  surface  runoff  in  deserts.  In¬ 
asmuch  as  the  Nahal  Yael  watershed  is  composed  of  mag¬ 
matic  rocks  like  much  of  the  major  deserts  of  the  world, 
it  seemed  reasonable  to  make  the  field  study  within  this 
watershed,  in  which  runoff  data  available  for  the  past 
15  years  offer  a  means  of  verifying  the  results. 

Second,  actual  infiltration  data  are  needed  as  an  input 
to  the  Nahal  Yael  runoff  and  sediment  model  which  is 
curently  being  developed  (see  chapter  9).  This  model  is 
based  on  very  small  unit  areas,  and  it  is  necessary  to 
understand  and  quantify  the  infiltration  process  on  a 
similar  level  of  detail. 

The  physiography  of  the  Nahal  Yael  watershed  has  been 
described  in  detail  in  several  sources.  The  lithologic 
composition  is  of  immediate  relevance  to  the  infiltra- 
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tion  tests  reported  here,  and  will  be  briefly  summari¬ 
zed.  Most  of  the  catchment  is  underlain  by  metamorphic 
rocks,  mainly  schist,  amphibolite,  and  some  gneiss 
(Shimron,  1974).  The  lower  part  of  the  catchment  is  com¬ 
posed  of  granite  of  the  Elat  type.  All  rocks  are  cros¬ 
sed  by  two  fairly  dense  systems  of  dykes.  Occasional 
colluvium  on  the  slopes,  which  are  mostly  very  steep 
(up  to  roughly  40°),  merges  into  alluvium  on  terraces 
above  the  larger  channels.  The  first  and  second  order 
channel  beds  are  almost  exclusively  bedrock.  Only  in 
some  lower  reaches  does  the  channel  bed  widen  to  more 
than  10  meters  and  is  completely  underlain  by  alluvi¬ 
um.  Such  channels  as  well  as  the  few  "soil"  pockets 
carry  the  only  vegetation  present  in  the  region  -  low 
shrubs.  Acacias  grow  only  in  the  l.wer  part  of  the 
main  channel . 

The  infiltration  trials  were  conducted  over  a  time  span 
of  nearly  two  years.  In  total,  19  plots  were  construc¬ 
ted,  and  24  infiltration  trials  were  accomplished  over 
8  field  periods. 

The  first  stages  of  the  field  work  were  afflicted  by 
difficulties  in  making  the  plot  borders  watertight. 

Only  after  we  included  in  the  procedure  a  routine  check 
by  applying  colored  water  to  the  border  strips,  could 
impermeability  be  verified  and  the  result  considered 
trustworthy.  Sometimes,  though,  plots  that  seemed  wa¬ 
tertight  prior  to  the  start  of  rainmaking,  developed 
leaks  during  the  test  or  towards  its  end,  probably  be¬ 
cause  of  a  difference  between  dry  and  wet  hydraulic 
conductivity  of  joints. 

Although  it  is  possible  to  obtain  an  adequate  descrip¬ 
tion  of  runoff  development  and  infiltration  rate  for 
certain  of  the  plots  which  developed  border  leaks,  we 
have  tabulated  and  used  only  the  data  for  those  tests 
which  proceeded  faultlessly  and  are  thus  completely 


-57- 


trustworthy.  The  only  item  of  information  utilized  from 
all  tests  was  the  total  amount  of  water  needed  for  the 
initiation  of  runoff,  for  rocks  other  than  amphibolite. 
This  variable  is  not  affected  by  the  faults  mentioned. 

8.2  The  infiltration  tests 

iD_  ill™  1 3.£l2Q 

The  rainmaking  device  (photographs  1-6)  consists  of  a 

o 

drop  generator,  area  0.8m  ,  which  applies  simulated 

2 

rain  to  a  bordered  micro-plot,  area  0.25m  .  The  drop 
generator  consists  of  a  rectangular  aluminum  frame 
along  whose  width  a  rubber  tube  winds  back  and  forth. 
Syringe  needles  (type  20x  1  |  G)  are  stuck  into  the 
tube  so  that  their  sharp  end  is  inside  and  the  wide 
end  points  downwards.  These  needles  form  the  drops, 
with  the  intensity  controlled  by  water  pressure  in  the 
tube.  A  constantly  moving  netting  installed  below  the 
needles  distributes  the  drops  over  some  area.  Constant 
pressure  is  ensured  by  a  constant  head  tank  (with  a 
float)  which  supplies  the  drop  generator,  and  is  sup¬ 
plied  in  turn  from  a  tank  located  at  a  higher  elevation 

The  device  is  levelled  by  means  of  four  adjustable  legs 
The  plot  is  delimited  by  low  fences  made  of  putty  (on 
bedrock)  or  sheet  metal  (on  colluvium),  and  has  a  fun¬ 
nel  at  its  lower  end  which  facilitates  the  collection 
of  runoff  into  containers.  Starting  with  the  time  of 
runoff  initiation,  the  containers  are  exchanged  every 
minute,  so  that  mean  infiltration  rates  for  60-second 
time  intervals  can  be  obtained. 

As  already  mentioned,  border  effectiveness  and  possible 
leaks  around  the  outlet  funnel  are  checked  by  the  appli 
cation  of  colored  water. 


Some  technical  data; 


Size  of  drops:  about  5mm 

Spacing  between  needles:  7x7  cm 

Size  of  device:  100x93  cm 

Size  of  plot:  50x50  cm 

The  last  two  measures  are  horizontal. 

The  rain  distribution  was  tested  in  a  closed  space 

2 

over  an  area  of  10  cm  ,  with  the  following  results: 


Test  No.  Rainage  No.  Cumulative  rainfall ,  mm 


1  1  3 

2  2  1/2 

2  1  2 

2  1  1/2 

3  1  3 

2  3  1/2 

3  4 


The  distribution  in  the  open  is  assumed  to  be  better 
because  of  wind  effect. 

Rainfall  intensity  was  determined  by  measuring  the  wa 

2 

ter  collected  in  a  measuring  tray,  size  0.5m  ,  i.e. 
identical  to  that  of  the  plot.  Three  measurements, 
one  minute  long  each,  were  carried  out;  two  were  made 
prior  to  the  exposure  of  the  plot  to  rain,  and  the 
third  was  made  at  the  end  of  the  infiltration  trial. 

The  advantages  of  the  rainfall  simulator  are: 

I)  mobility:  the  device  enables  tests  in  areas  unsuita 
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ble  for  vehicular  transport.  It  can  be  operated  by  two 
persons  (though  three  are  better)  and  is  carried  from 
site  to  site  by  the  operators. 

II)  water  economy:  water  is  the  heaviest  component  of 
the  simulator.  At  distances  of  2km  from  the  vehicle, 
the  motto  "conserve  every  drop"  acquires  a  very  special 
significance!  The  device  uses  10-15  liters  to  prepare 
it  for  a  set  of  tests  (this  can  be  done  near  the  ve¬ 
hicle)  and  then  20-30  liters  for  each  infiltration  test; 
the  exact  amount  depends  on  the  time  needed  to  achieve 
the  terminal  infiltration  rate  and  on  possible  faults 

in  the  device  or  in  the  plot  during  the  test. 

III)  selectivity  of  rain  intensity:  the  device  performs 
accurately  in  the  intensity  range  of  15-120mm/hr,  and 
maintains  a  uniform  intensity  (the  devi'Hons  are  1-3%) 
with  a  reasonable  drop  distribution.  Switching  to  dif¬ 
ferent  size  needles  should  enable  the  simulation  of 
rain  of  any  desired  intensity. 

IV)  versatility:  the  device  can  be  used  in  the  laborato¬ 
ry  as  well  as  in  the  field.  It  is  adaptable  to  a  slope 
range  of  0-100%. 

V)  The  device  is  simple  and  inexpensive  to  construct. 

Some  disadvantages  of  the  rainfall  simulator  are  con¬ 
nected  with  the  drop  generator:  the  drops  are  not  ide¬ 
ally  distributed,  they  are  larger  than  those  of  natural 
rain,  and  their  fall  velocity  is  smaller.  The  first  two 
factors  can  be  disregarded  for  rainmaking  on  rocky  sur¬ 
faces,  but  the  third  one  may  indeed  be  a  problem. 

Other  disadvantages  concern  the  plots:  the  waterproo¬ 
fing  in  jointed  rocks  and  in  coarse  colluvium  poses  con¬ 
siderable  problems,  especially  on  steep  slopes.  In  the 
putty-fenced  plots  the  border  strip  is  appreciably  dis¬ 
turbed:  the  inner  part  of  the  fence,  mean  width  0.75  cm, 
produces  an  area  of 
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0.75  x  200  =  150cm2 

2 

i.e.  approximately  6%  of  the  plot  area  of  2500cm  . 

For  reasons  we  cannot  explain,  it  is  not  possible  to 
use  the  device,  over  two  days,  for  more  than  ten  tests. 
Thereafter  the  simulator  begins  to  perform  erratically, 
and  the  intensity  varies  within  a  single  test. 

Another  limitation  of  the  device  is  its  requirement  of 
distilled  or  filtered  water,  to  avoid  the  clogging  of 
the  syringe  needles. 

§;?i?__I§§5.!2§sylts 

The  results  of  13  infiltration  tests  conducted  on  eight 
plots  are  listed  in  Tables  8.1  to  8.11.  The  location  of 
these  plots  is  shown  on  fig.  8.1. 

§i?;.§__Ih§_2orrectionof_the_measured_infiltration_rates 

In  the  tests  the  plot  infiltration  rate  was  calculated 
by  subtracting  the  runoff  collected  in  the  outlet  fun¬ 
nel  from  the  plot  rainfall, over  successive  periods  of 
one  minute  duration.  The  value  of  this  infiltration  rate 
was  corrected  according  to  the  following  criteria: 

I)  For  sloping  plots  the  splash  balance  is  not  zero  be¬ 
cause  the  drop  generator  is  horizontal,  i.e.  at  the 
lower  end  the  drops  fall  from  a  higher  elevation.  Con¬ 
sequently  the  loss  of  water  to  the  plot  by  the  down- 
splash  component  of  the  bottom  inner-border  strip  is 
larger  than  the  gain  to  the  plot  by  the  downsplash  com¬ 
ponent  of  a  top  outer  border  strip  of  identical  width. 
The  width  of  the  rained-on  outer  margins  of  the  plot  is 
about  20cm. 

An  estimate  of  the  loss  of  water  to  the  plot  caused  by 
the  above  factor,  based  on  data  on  the  number  of  drops 
and  their  fall  velocity,  yielded  the  following  correc- 


rainfall  recorder 
stream  gauging  station 
[7J  infiltration  test  site 


x&s, 


05V140I  / 


Fig.  8.1:  Nahal  Yael  infiltration  test  sites  and  relevant  instru 
mentation. 


Table  8.1.  :  Infiltration  test  at  site  No.  7 

Photograph  No. 11  Figure  No.  8.2 


Photograph  No.  11 

L i thology:  granite,  medium  grain, 
un.iointed  _  _ 


Area  of  plot,  horizontal:  2480  crn 
Runoff  begins  after:  0  min  58  sec 


Date  of  test:10.  9  .78 
Slope:  2° 

Ha  in  int.:  36  mm/hr 
41 


T  ime 
min 


Runoff 


55.1 

93.1 


93.4 


Inf.  rate 

mm/hr 

Comments 

18.2 

rain  36  mm/hr 

7.6 

II 

5.0 

II 

2.8 

rain  41  mm/hr 

3.5 

0.62 

0.87 

0.38 

0.62 

Terminal  conditions 

Stable  infiltration  rate:  Rai nfal 1 -runoff  ratio 

.  62  mm/hr  98.4  per  cent 


-64- 


Table  8.2.  :  Infiltration  test  at  site  Mo.  15 


Photograph  No .  9 


Figure  No.  8.3 


Lithology:  granite,  coarse-grained.  Date  of  test:  8.  1.79 

sparsely  jointed _  Slope:  34  0 

Rain  int. :  31. 4  mm/hr 
32. 2  „ 


Area  of  plot,  horizontal:  2355  cm^  Temp. water. 
Runoff  begins  after:  .  1  min  12  sec 


20  °C 


T  ime 

mi  n 

Runoff 

Inf.  rate 

mm/hr 

Comments 

cm^ 

% 

1 

.  2 

90 

73.1 

10.5 

rain  31.4  mm/hr 

3 

95 

7.1 

It 

4 

110 

3.6 

It 

5 

114 

89.7 

3.2 

rain  32.2  mm/hr 

6 

116 

2.8 

II 

7 

II 

8 

II 

9 

II 

10 

II 

11 

116 

91.3 

2.8 

II 

12 

116 

2.8 

II 

13 

14 

15 

16 

17 

18 

19 

20 

Terminal  conditions 

Stable  inf iltra ti on  rate :  Ralnf al 1 -runoff .ratio: 

91.3  per  cent 


2  . 8  mm/hr 


Figs.  8.2  -  8.5:  Nahal  Yael  infiltration  tests  -  granite,  schist, 
colluvium.  P  =  precipitation;  RITS  ■  Rainfall 
amount  Infiltrated  Till  Stability  achieved. 


r  1 

Li _  i 
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Tab  1  e  8.  3 


Infiltration  test  at  site  No. 3 


Photograph  No . 10 

Lithology:  schist,  densely,  jointed; 
joints  perpendicular  to  slope _ 


Figure  No.  8.4 

Date  of  test:  8 .  1 .79 

Slope:  33° 

Rain  int.:  25.2  mm/hr 
35.5 


Area  of  plot,  horizontal:  2316  cm 
Runoff  begins  after:  3  min  05  sec 


Temp:water:  •  20°C 


T  ime 

min 

Runoff 

Inf.  rate 

mm/hr 

Comments 

3 

cm 

% 

1 

2 

3 

4 

5 

38 

30.0 

rain  25.2  mm/hr 

6 

34.0 

16.6 

7 

14.2 

8 

intensity  increased 

9 

71 

17.1 

rain  35.5  mm/hr 

10 

95 

69.3 

10.9 

II 

11 

108 

7.5 

II 

12 

117 

5.2 

II 

13 

116 

84.6 

5.4 

II 

14 

124 

3.4 

II 

15 

no  data 

16 

no  data 

17 

104 

8.6 

rain  35.5  mm/hr 

18 

120 

4.4 

II 

19 

no  data 

20 

no  data 

21 

no  data 

22 

no  dat* 

23 

no  data 

24 

126 

91.9 

2.8 

rain  35.5  mm/hr 

25 

119 

86.9 

4.6 

M 

Terminal  conditions 

ll5H25i2D_!I3te :  Rai  nf  f -E§£i2 : 


3.8  mm/hr  89  per  cent 
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Table  8.4  :  Infiltration  test  at  site  Ho.  1 

Photograph  No.  7 

Lithology:  Colluvium  in  an  area  of  Date  of  testill.  8-78 

amphibolite  and  dykes;  stone  cover  Slope:  11° 

85%;  angular  stones,  size  5-10  mm.  Rain  int.:  51.6  mm/hr 
some  up  to  40  mm _ 

2 

Area  of  plot,  horizontal:  cm 

Runoff  begins  after:  3  min  20  sec* 


Terminal  conditions 

13-6  mm/hr  73.6  per  cent 


*  plot  was  dry  at  onset  of  rain 


-68- 


Table  8.5  :  Infiltration  test  at  site  No.  5 


Photograph  No.  8 


Figure  No.  8.5 


L i thology:  colluvium  in  an  area  of 
schist  and  dykes;  stone  cover  95%; 
angular  stones,  size  10-40  mm _ 


Date  of  test:  9.  1.79 
Slope:  9° 
Rain  int.  :  30.7  mm/hr 


Area  of  plot,  horizontal:  2354  cm*"  Temp. water 
Runoff  begins  after:  5  min  40  sec 


T  i  me 

m  i  n 

Runoff 

Inf.  rate 

mm/  hr 

Comments 

3 

cm 

c' 

/: 

1 

2 

3 

' 

4 

. 

5 

6 

7 

6 

8 

16 

13.2 

26.7 

9 

30 

23.2 

10 

34 

22.2 

11 

40 

20.7 

12 

no  data 

13 

54 

42.6 

17.1 

14 

58 

16.1 

15 

67 

13.7 

16 

71 

12.7 

17 

79 

65.2 

10.7 

18 

78 

11.0 

19 

79 

10.7 

20 

78 

11.0 

21 

82 

67.7 

10.0 

Terminal  conditions 

Stable  infiltration  rate:  Ra i nfal 1 -runoff  ratio: 


10.  5  mm/ hr 


65  per  cent 
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Table  8.6  :  Infiltration  test  at  site  No.  8 


Photograph  No .  12 

L  i  1. 1 1 o  1  o g y  :  amphi  bol  i  te  ;  dense 
joints _ _ _ 


Area  of  plot,  horizontal  :  2022  c in ^ 

Runoff  begins  after:  0  min  50  sec 


Figure  No.  .8.  6 

Date  of  test:  12.  8.78 

Slope:  20° 

Rain  int.:  56.4  mm/hr 
55.2 
54.0 


Time 

m  i  n 

R  u  n  o  f  f 

Inf.  rate 

mm/  hr 

Comments 

3 

cni 

c- 

ft 

1 

O 

L. 

107 

24.6 

rain  56.4  mm/hr 

3 

139 

68.4 

17.6 

II 

4 

145 

13.3 

II 

i 

147 

11.6 

rain  55.2  mm/hr 

(i 

152 

10.1 

II 

7 

152 

SI  .7 

10.1 

II 

O 

155 

8.1 

rain  54.0  mm/ hr 

9 

131 

15.1 

II 

10 

159 

S7.3 

6.8 

II 

I  1 

12 

1  3 

14 

15 

16 

17 

1  hi 

1  9 

|  20 

Terminal  conditions 


It  i  doubtful  whether  stability  has  indeed  been  reached 
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Table  8.  7  :  Infiltration  test  at  site  No. 20 


Photograph  No.  13 

Lithology :  amphibolite;  medium 
dense  jointing _ 


Figure  No.  8.7 

Date  of  test:  9  .  4.79 
Slope:  18° 
Rain  int.:  40.8  mm/hr 


Area  of  plot,  horizontal:  2665  cm^  Temp. water 

Runoff  begins  after:  1  min  45  sec  soil 


26  °QC 
23°C 
31°C 


Time 

Runoff 

Inf.  rate 

Corrected  # 
inf.  rate* 
mm/hr 

m  i  n 

i 

3 

cm 

% 

mm/  hr 

1 

2  1/2 

44 

24.3 

30.84 

33.4 

3 

90 

20.48 

22.1 

4 

97 

18.91 

20.1 

5 

102 

56.3 

17.78 

18.6 

6 

102 

17.78 

18.6 

7 

113 

15.30 

15.7 

8 

125 

12.60 

12.5 

9 

114 

62.9 

15.0 

15.2 

10 

119 

13.95 

14.0 

11 

112 

15.53 

15.7 

12 

122 

67.4 

13.28 

13.2 

13 

118 

14.18 

14.2 

14 

121 

66.8 

13.5 

13.5 

15 

16 

17 

18 

19 

20 

.  _ 

Terminal 

Stable  infiltrationrate: 

13.5  mm/hr 


conditions 

B§infall-runoff_ratio: 

66.8  per  cent 


"see  text 


Table  8.  8 


Infiltration  test"at  site  No. 20 


Photograph  No .  13 

Lithology:  as  in  table  8.7 _  Date  of  test:  9.  4.79 

_  Slope:  18° 

Rain  int. :  19.9  mm/hr 


Area  of  plot,  horizontal:  2665  cm^  Temp. water. 

Oil  y,  • 

Runoff  begins  after:  2  min  30  sec  soil': 


29  °C 
28  °C 
24  °C 


Time 

Runoff 

Inf.  rate 

Corrected 
inf.  rate 
mm/hr 

m  i  n 

cm^ 

a 

/? 

mm/hr 

1 

2 

3 

25 

28.4 

14.2 

15.5 

4 

44 

9.9 

10.8 

5 

58 

6.7 

7.3 

& 

64 

5.4 

5.9 

7 

64 

5.4 

5.9 

8 

64 

72.7 

5.4 

5.9 

9 

62 

5.8 

6.4 

10 

66 

4.9 

5.4 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

- 

Terminal  conditions 

5.  9  mm/hr  72.7  per  cent 

this  is  a  repeat  test  conducted  3  hours  after'  the  first 
test  (table  8.7).  Note  the  difference  in  temperature  bet 
ween  wet  and  dry  soil.  The  reason  for  the  low  terminal 
rate  is  not  clear. 
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Table  8.9  :  Infiltration  test^at  site  No.  20 


Photograph  No . 13 


Figure  No.  8.8 


Lithology :  as  in  table  8.7 


Date  of  test:  10.  4.79 
Slope:  18° 
Rain  int. :  18 .7  mm/hr 


2 

Area  of  plot,  horizontal:  2665  cm 
Runoff  begins  after:  2  min  00  sec 


Temp :wa ter 
air 
soil* 


21°C 
20  °C 
17  °C 


T  ime 

Runoff 

Inf.  rate 

Corrected 
inf.  rate 
mm/hr 

mi  n 

3 

cm 

% 

mm/hr 

1 

2 

3 

25 

30  • 

1 

14.1 

4 

34 

■RSv^; 

11.9 

5 

38 

10.1 

10.9 

6 

43 

9.0 

9.7 

7 

48 

57 

9.9 

10.5 

8 

50 

7.4 

8.0 

9 

52 

7.0 

7.5 

10 

52 

62 

7.0 

7.5 

11 

52 

7.0 

7.5 

12 

52 

7.0 

7.5 

13 

14 

15 

16 

17 

18 

19 

20 

54 

65 

. 

6.5 

7.0 

Terminal  conditions 


Stable  i nf i 1 tra ti onra te :  Ra  inf  5  Hi™  D2f  • 

7.  5  mm/hr  62  per  cent 

this  is  a  repeat  test  conducted  22  hours  after  the  first 
test  (table  8.7)  and  19  hours  after  the  second  test  (ta¬ 
ble  8.8).  Note  the  difference  in  temperature  between  dry 
and  wet  soil. 
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Table  8.10  :  Infiltration  test“at  site  Ho. 20 

Photograph  No. 13 


Lithology:  as  in  table  8.7 

Date  of  test: 

Slope : 

Rain  int.:  42 

10.  4.7  9 
18  0 
mm/hr 

2 

Area  of  plot,  horizontal:  2665  cm 

Temp:water: 
air  : 

23°C 

22°C 

Runoff  begins  after:  0  min  25  sec 


T  iir.e 

Runoff 

Inf.  rate 

Corrected 
inf.  rate 
mm/hr 

mi  n 

cni^ 

% 

mm/hr 

1 

mm 

70.9 

12.2 

13.3 

2 

■fta 

9  .9 

10.7 

3 

144 

9.5 

4 

145 

77.9 

9.2 

9.6 

5 

148 

8.6 

8.6 

6 

148 

79.5 

8.6 

8.6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

144 

9.5 

10.1 

Terminal  conditions 

5table_infiltration_rate:  B?l0f2lllCyD2ff_!C5^l2: 

8  .  6  mm/hr  79.5  per  cent 

this  is  another  repeat  test,  performed  one  hour  afterthe 
previous  one  (table  8.9). 
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Table  8.11  :  Infiltration  test  at  site  No. 40 


Photograph  No. 14 


Figure  No.  8.9 


Lithology :  amphibolite,  medium  .join¬ 
ting.  Rain  deliberately  stopped 
twice  during  test. _ 


2 

Area  of  plot,  horizontal:  2550  cm 
Runoff  begins  after:  0  min  55  sec 


Date  of  test:  19.  9 .79 


Rain  int. :  21 .6  mm/hr 

19.4 

26.4 


Temp: water 
air 
soil 


23.5 

23 

25 


T  ime 

Runoff 

Infiltration  rate 

Rai  n 
int. 

mi  n 

3 

cm 

% 

direct 

mm/hr 

corrected 

mm/hr 

mm/hr 

1 

2 

16 

17.5 

■9 

— 

21.6 

3 

24 

21.6 

4 

26 

16.1 

21.6 

5 

26 

28.5 

15.5 

16.1 

21.6 

6 

28 

15.0 

15.5 

21.6 

7 

26 

15.5 

16.1 

21.6 

8 

31 

14.3 

14.8 

21.6 

9 

28 

15.0 

15.5 

21.6 

10 

33 

36.2 

13.8 

14.2 

21.6 

11 

33 

13.8 

14.2 

21.6 

12 

39 

12.4 

12.7 

21.6 

13 

36 

39.5 

13.1 

13.5 

21.6 

14 

15 

16 

36 

13.1 

13.5 

21.6 

17 

37 

40.5 

13.3 

13.7 

21.6 

TB 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

DTD 

continued  on  following  page 
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Table  8.  11 


Infiltration  test  at  site  No. 40  (contl) 


T  ime 

Runoff 

Infiltration  rate 

min 

cm^ 

* 

direct 

corrected 

mm/hr 

mm/hr 

31 

0.0 

32 

II 

33 

34 

82 

9.9 

10.5 

19.4 

35 

47 

56 

8.5 

8.9 

II 

36 

51 

7.5 

7.9 

II 

37 

50 

7.8 

8.1 

II 

38 

39 

50 

60 

7.8 

8.1 

II 

4  C 

0.0 

41 

42 

43 

44 

45 

46 

47 

48 

76 

68 

8.4 

9.1 

26.4 

■ 

72 

9.4 

10.2 

II 

76 

8.4 

9.1 

II 

»|M|i 

76 

8.4 

9.1 

N 

79 

70.5 

7.7 

8.4 

n 

Terminal  conditions 

Stable  infiltration  rate!  Ra i nf al 1 -runoff  ratio: 


Stage  1:  13.5  mm/hr  39.5  per  cent 

Stage  2:  8-  1  n  60  *»  » 

.Stage  3:  9.  1  »  68  n  n 
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tion  coefficients:  For  plots  sloping  30-35°,  deduct  5% 
from  incoming  rainfall;  for  plots  sloping  15-20°,  de¬ 
duct  2%  from  incoming  rainfall;  for  plots  of  lesser  slo¬ 
pe,  no  correction  is  warranted. 

II)  Disturbance  of  margins:  great  care  was  exercised  to 
construct  the  micro-fence  bordering  the  plots  in  such 

a  way  that  the  horizontal  area  of  its  inward  sloping 

side  becomes  as  small  as  possible.  The  fence  walls  were 

slightly  bent  inward  so  as  to  "cover"  the  putty  strip 

at  its  base.  However,  in  certain  plots  the  putty  strip 

reached  the  substantial  width  of  up  to  1  cm  mean  value, 

2 

totalling  200  cm  per  plot  or  roughly  8  per  cent  of  its 
area . 

Since  the  putty  strip  produces  100%  runoff,  the  water 
amount  contributed  by  it  was  deducted  from  the  total 
collected  at  the  outlet  funnel,  and  the  plot  area  was 
also  reduced  by  subtracting  the  area  of  the  putty  strip. 
This  correction  was  applied  to  the  two  amphibolite  plots 
(No. 20  and  No. 40),  in  which  it  was  expected  to  exercise 
a  considerable  influence  on  the  result.  However,  the 
difference  proved  to  be  slight.  Consequently,  this  fac¬ 
tor  was  neglected  for  all  other  plots. 

III)  The  infiltration  tests  were  conducted  over  diffe¬ 
rent  seasons,  consequently  different  temperature  values 
for  the  simulated  rainwater  as  well  as  for  the  air  and 
the  soil  should  be  considered.  Although  the  effect  of 
temperature  on  infiltration  has  been  well  documented 
for  some  cases  (Shanan  1975;  and  others),  the  complete 
relationship  is  far  from  clear  and  may  be  variable  from 
one  soil  type  to  another.  We  decided,  arbitrarily,  to 
apply  a  conservative  correction:  in  order  to  standar¬ 
dize  the  test  results  for  infiltration  of  cold  rain 
water  in  winter,  we  have  reduced  the  terminal  infiltra¬ 
tion  rates  of  tests  carried  out  in  winter  with  the  artifi- 


r 
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cially  warmer  water  by  10  per  cent;  for  all  other 
tests  the  reduction  was  fixed  at  20  per  cent. 

Table  8.12,  the  data  summary  sheet  of  all  infiltration 
tests,  lists  also  the  values  as  corrected  by  the  above 
procedures . 

§^2 . 4 _ iQf  Utrati  on  _into_  joints 

In  highly  jointed  areas  the  tests  described  above  did 
not  yield  accurate  results  because  it  was  impossible  to 
adequately  waterproof  the  joints  located  under  the 
plot  boundaries. 

Two  sets  of  experiments  were  carried  out  in  order  to  evaluate 
the  infiltration  rate  into  joints: 

I)  The  application  of  simulated  rainfall  on  an  unboun¬ 
ded,  jointed  surface  and  the  subsequent  determination 
of  the  time  of  runoff  initiation  and  of  the  infiltra¬ 
tion  rate  into  the  joints.  The  latter  figure  was  appro¬ 
ximated  by  dividing  the  depth  of  wetting  by  the  dura¬ 
tion. 

II)  The  application  of  water  to  an  intersection  of 
joints  at  a  rate  exceeding  the  infiltration  rate,  for 
a  relatively  long  duration,  and  the  approximation  of 
the  infiltration  rate  as  above. 

InfiItration_test_on_an_unbounded_surface:  The  surface 
selected  for  this  test  is  composed  of  schists  with  wide 
joints  filled  with  fine-grained  material.  The  total  rain¬ 
fall  amounted  to  3.3  mm  and  was  applied  by  the  same  de¬ 
vice  used  in  the  tests  described  earlier.  The  results 


were : 


Nahal  Yael  -  data  sunnary  sheet 


of  20  m/hr  of  rainfall  **  based  on  an  earlier  test 
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Time 

(min) 


Rain  inten-  Observation  on  Comments 

sity  runoff 

(mm/hr) 


0  -  4 


31  Runoff  starts  af¬ 

ter  a  few  seconds 


4-7 


7  -  14 


0  Water  continues  to  rain  stopped 

stand  in  depressi¬ 
ons  formed  by  the 
joints 

12  Runoff 


Ten  minutes  after  the  termination  of  rain,  the  depth  of 
wetting  in  most  points  was  observed  to  be  1-2  cm;  in 
the  wide  joints  -  4.5  cm.  Horizontal  joints,  i.e.  those 
partitioning  pseudostrata,  at  a  depth  of  2-3  cm  were 
all  dry.  This  means  that  no  lateral  flow  in  depth  occur' 
red . 

The  mean  depth  of  wetting  was  determined  to  be  roughly 
3  cm  per  25  min  (till  the  time  of  exposure),  i.e. 

1.2  mm/min  (72  mm/hr).  But  since  the  water  occupies  on¬ 
ly  50%  of  the  joint  volume*  ,  the  actual  rate  is  only 
0.6  mm/min  (36  mm/hr).  Areally  the  joints  occupy  about 
20  per  cent;  therefore  the  rate  of  infiltration  into  the 
joints  is  only  one  fifth.  Accepting  that  the  rock  it¬ 
self  absorbs  only  negligible  amounts  of  water  (Lavee 
1973),  we  arrive  at  an  infiltration  rate  of  7  mm/hr. 


During  vertical  movement  within  joints  under  such  wet 
conditions,  the  "soil"  within  the  crack  can  be  assumed 
to  be  supersaturated;  the  "soil"  texture  is  a  silty 
sand,  with  a  one  cm  thick  vesicular  layer  at  the  top, 
and  a  compact  or  platy  (parallel  to  the  crack  surfaces) 
structure . 
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^2iDt§_ill?®!I§§9^l2Q_£§§^_^2;I  (same  lithology  -  schists 
with  wide  joints):  these  rain  simulation  tests  were  de¬ 
signed  solely  to  observe  the  rate  of  advance  of  the 
wetting  front  in  the  joints.  The  rain  intensity  excee¬ 
ded  by  far  the  infiltration  rate  into  the  joints,  and 
considerable  runoff  occured.  The  surface  depressions 
were  filled  to  capacity  and  created  a  certain  hydraulic 
head,  probably  no  different  from  the  conditions  during 
natural  rain,  as  observed  during  the  test  on  an  unboun¬ 
ded  surface  described  earlier. 

Rain  equivalent  to  0.140  liters  was  applied  over  a  pe¬ 
riod  of  3  minutes  to  an  area  size  8x20  cm.  The  depth  of 
the  wetting  front  was  determined  when  the  water  dis¬ 
appeared  from  the  mi cro-depressions ,  i.e.  5  minutes  af¬ 
ter  the  termination  of  the  rainfall,  to  be  2  cm  (ave¬ 
rage  value).  Considering  that  the  reading  was  made  8 
minutes  after  the  beginning  of  the  rainfall,  this  indi¬ 
cates  a  mean  rate  of  advance  for  the  wetting  front  of 
2.5  mm/min.  Based  on  the  consideration  detailed  earlier, 
the  infiltration  rate  indicated  for  the  joints  surface 
area  is  1.25  mm/min,  or  75  mm/hr. 

1?2  2Dt§_l2t§E§§2t  l2D_£§st_No;2  (a  densely  jointed  dyke 
rock  surface):  As  a  result  of  a  5  minute  long  applica¬ 
tion  of  0.400  liters  of  simulated  rain,  water  flowed  in 
two  microchannels  conditioned  by  joints  at  their  bottom. 
Ten  minutes  after  the  onset  of  rainfall  the  wetting 
front  was  found  to  have  advanced  to  a  depth  of  2  cm, 
on  average,  i.e.  a  rate  of  120  mm/hr.  This  means,  as 
above,  a  water  equivalent  of  60  mm/hr  for  the  joint  sur¬ 
face  area. 

d2i25s_i ntersecti on_ tes t_No . 3  (same  site  as  No. 2):  This 
test  was  carried  out  at  a  T-intersection  of  3  wide 
joints  (width  1-3.5  mm).  Water  was  applied  at  a  relative¬ 
ly  very  high  discharge  -  7.3  liters  over  20  minutes  - 
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from  a  tap.  Three  minutes  after  the  termination  of  the 
water  input,  i.e.  23  minutes  after  its  beginning,  the 
following  data  were  observed: 

In  the  joints  above  which  the  water  flowed,  the  depth 
of  the  wetting  front  ranged  from  3-4  cm  (minimum)  to 
12  cm  (maximum);  mean  value  8  cm.  Therefore,  the  water 
equivalent  infiltration  rate  is,  roughly,  45  mm/hr 
(minimum),  150  mm/hr  (maximum),  and  100  mm/hr  (mean 
value)  . 

Summary_of _tests : 

The  mean  water  equivalent  infiltration  rate  of  all  the 
tests  described  above  is  67  mm/hr.  If  the  last  test,  in 
which  exceptionally  large  amounts  of  water  were  applied, 
is  ignored,  the  value  becomes  55  mm/hr.  We  conclude 
that  the  infiltration  rate  of  60  mm/hr  represents  well 
the  behavior  of  joints. 

Remembering  the  negligible  absorption  of  water  during 
rainstorms  by  rocks  such  as  granite,  schist,  and  quartz 
porphyry  (Lavee  1973),  the  dominant  factor  which  con¬ 
trols  infiltration  in  a  rocky  area  is  the  number  of 
joints  per  unit  area  and  their  width;  or,  the  propor¬ 
tion  of  joint  surface  area  in  the  entire  area  conside¬ 
red.  The  schist  site  tested  is  very  highly  jointed 
(photograph  15),  and  its  joints  occupy  20  per  cent  of 
its  area.  However,  at  a  depth  of  1  cm  and  beyond  the 
aggregate  joint  surface  area  diminishes  to  between 

5  and  10  per  cent.  Thus  the  infiltration  rate  of  such 
an  area  in  its  entirety  is  only  3-6  mm/hr.  In  normally 
jointed  terrain  (table  8.13),  the  cracks  at  depths  of 
1  cm  and  lower  occupy  only  2-3 %  of  the  total  surface 
area,  and  the  areal  infiltration  rate  is  only  1-1.5 
mm/hr.  However,  some  of  the  slope  terrain,  in  particu¬ 
lar  of  granitic  lithology,  is  crossed  by  very  wide 
joints  which  produce  a  veritable  microtopography  of 
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Table  8.13:  Summary  of  infiltration  tests  at  joints 


Rock  type 

Number  of 

jointing 

directions 

Approx, 
width  of 
joints 

(mm) 

Proportior 
of  joint 
area 

(*) 

Comments 

Dyke  rock 

2 

0.5 

2.25 

8-10  joints 
over  40  cm 

Schist 

1  * 

0.5 

3 

'cleavage;  12 
joints  over 

20  cm 

A 

Schist 

1  " 

0.5 

2 

perpendicular 
to  cleavage;  8 
joints  over  20 
cm 

Granite 

1 

1 

1  5 

6  joints  over 

40  cm 

Granite 


2 


2 


2 


3  and  5  joints, 
respectively, 
over  40  cm  each 


-86- 


rills  and  depressions;  these  must  cause  a  considerable 
increase  in  the  infiltration  capacity  of  the  entire 
slope . 

Valuation:  In  the  tests,  the  relation  between  infil¬ 
tration  into  the  joint  and  the  width  of  the  joint  is 
not  treated,  though  some  influence  can  be  expected. 
Likewise,  the  effect  of  rain  intensity  is  disregarded, 
but  this  factor  appears  to  be  of  no  consequence,  as  any 
rainfall  amount  which  generates  runoff  produces  also 
storage  depression  which  covers  the  joints  with  water. 

Observations  during  the  tests  suggested  that  infiltra¬ 
tion  proceeds  very  rapidly  into  the  top  centimeter 
characterised  by  the  vesicular  structure,  and  slows 
down  considerably  below  that  depth,  when  the  wetting 
front  reaches  denser  (platy  or  massive)  material.  There 
for  the  infiltration  rate,  after  a  few  millimeters  of 
initial  rainfall,  is  probably  even  lower  than  the  mean 
values  as  stated  above.  Observations  in  cuts  in  magma¬ 
tic,  volcanic,  and  metamorphic  rocks  disclose  that  most 
joints  virtually  disappear  at  a  depth  of  several  deci¬ 
meters;  at  most  they  become  so  thin  that  no  fill  mate¬ 
rial  is  present  and  hence  no  water  can  pass  through 
them . 

In  conclusion,  after  being  wetted  by  a  continuous  rain 
of  several  hours'  duration,  or  by  a  series  of  several 
shorter  rainstorms,  the  joints  become  saturated  over 
their  entire  depth,  at  which  time  the  infiltration  rate 
of  rocky  jointed  terrain  becomes  nearly  equal  to  that 
of  rocky  unjointed  terrain. 


The  infiltration  rates,  as  determined  by  the  above  tests 
are  lower  by  a  factor  of  2-3  than  the  infiltration  rate 
obtained  in  a  plot  test  of  a  similar  terrain  (plot  No. 3 
table  8.3.).  The  discrepancy  can  be  explained  by  the  re 
latively  low  level  of  accuracy  inherent  in  the  type  of 
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tests  used  for  the  joints.  However,  as  indicated  by  the 
plot  tests  and  independently  substantiated  by  the  "joint" 
tests,  the  infiltration  rate  in  rocky  jointed  terrain 
is  much  lower  than  might  be  expected  from  the  visual 
surface  impression  of  widely  gaping  cracks  with  but 
1 i ttle  fil 1  material . 

8.3  Infiltration  rates  in  various  lithologies 

Granitic  slopes  constitute  only  the  lower  part 
of  Nahal  Yael  ,  and  are,  therefore,  of  little  importance 
as  runoff  contributors  to  the  watershed.  Granite  is, 
however,  very  widespread  in  the  Sinai  as  well  as  in  other 
deserts  elsewhere. 

The  granite  at  Nahal  Yael  is  composed  of  feldspars, 
quartz,  and  biotite,  and  is  medium  to  coarse  (0.5  - 
3  mm).  Some  slopes  are  smooth,  others  are  jointed  with 
a  fairly  sharp  microrelief. 

Out  of  four  tests  carried  out  on  granite,  two  (plot 
No. 7,  table  8.1  and  plot  No. 15,  table  8.2)  proved  de¬ 
pendable.  Plot  7  represents  medium,  unjointed,  resis¬ 
tant  and  smooth  granitic  terrain.  Plot  15  represents 
coarse-grained  granite,  somewhat  jointed  and  slightly 
crumbly  with  a  very  rough  surface.  Averaging  the  re¬ 
sults  of  the  two  plots,  we  conclude  that  0.6  mm  of  rain¬ 
fall  are  necessary  for  the  infiltration  of  runoff;  the 
stable  infiltration  rate,  standardized  to  a  20  mm/hr 
winter  rainfall,  is  1.3  mm/hr. 

It  seems  that  the  infiltration  rate  obtained  for  plot  15 
is  not  the  terminal  one,  the  reason  being  that  in  depth 
the  rock  becomes  harder  and  less  weathered. 

Data  for  four  natural  rainfall  events,  measured  on  a 
steep,  coarse  grained,  thinly  jointed  granitic  site 
near  the  Santa  Katherina  monastery  in  the  Sinai  (A.  Yair, 
personal  communication)  suggest  that,  for  a  rainfall 
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intensity  of  8  mm/hr,  the  mean  infiltration  rate  is  on 
the  order  of  2  mm/hr.  For  the  event  of  11  December  1978 
in  Nahal  Yael ,  a  little  runoff  began  to  flow  on  the 
smooth  granite  slopes  already  at  the  intensity  of  2.2 
mm/hr  (after  an  initial  wetting  of  3  mm);  however,  no 
flow  was  observed  on  jointed  slopes  at  the  same  time. 

Only  when  the  rainfall  intensity  increased  to  4.35  mm/hr, 
considerable  runoff  resulted  on  both  types  of  slopes. 

In  conclusion,  granitic  slopes  generate  runoff  quickly; 
the  infiltration  rate  is  low  and  the  runoff  high.  For 
a  shower  whose  intensity  exceeds  the  infiltration  rate, 
the  granite  absorbs  1-1.25  mm  of  water  over  5-6  minutes 
prior  to  reaching  the  stable  infiltration  rate;  this 
value  is  suggested  for  use  in  short  rainstorms  (lasting 
from  some  minutes  to  a  few  tens  of  minutes).  For  larger 
storms  the  rate  probably  decreases. 

The  terminal  infiltration  rate,  for  a  medi um-i n tens i ty 
rain  (10-20  mm/hr)  is,  approximately : 

Unjointed  granite  1  mm/hr 

Thinly  jointed  granite  1.5-2  mm/hr 

Medium  to  densely- 

jointed  granite  2. 5-3. 5  mm/hr 

Schist^  Most  of  the  Nahal  Yael  catchment  is  schist.  The 
schist  is  pelitic,  highly  foliated,  and  is  composed 
mainly  of  biotite,  muscovite,  plagioclase,  and  quartz 
(Shimron,  1974).  "Dip"  slope  schists  produce  smooth  sur¬ 
faces,  sometimes  jointed  downslope,  and  generate  runoff 
rapidly;  their  terminal  infiltration  rate  is  low.  Slopes 
with  schist  cleavage  planes  which  are  inclined  opposite 
to  the  slope  produce  very  rough  and  densely  jointed  ter¬ 
rain  surfaces,  characterised  by  a  very  high  depression 
storage  as  well  as  a  high  infiltration  capacity. 

The  second  type  is  represented  by  plot  3  (table  8.3). 
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The  trials  indicate  that  runoff  begins  only  after  the  fall 
of  1.26  mm  of  rain  ,  and  that  the  terminal  infiltration  ra¬ 
te  (standardized  as  above)  is  2.9  mm/hr. 

Tests  conducted  on  a  "dip"-slope  schist  site  failed  due 
to  our  inability  to  waterproof  the  wide  joints  under 
the  plot  boundaries.  Our  qualitative  impression  was 
that  the  terminal  infiltration  rate  of  these  slopes  was 
much  lower,  perhaps  1-1.5  mm/hr. 

Amphibolite:  In  the  simulated  rainfall  experiments  spe¬ 
cial  attention  was  devoted  to  the  infiltration  capacity 
of  the  amphibolite,  because  most  of  watershed  05  is  un¬ 
derlain  by  it,  as  well  as  because  it  is  so  different 
from  the  other  rocks  investigated.  The  amphibolite  of 
Nahal  Yael  is  a  massive  rock  of  a  fairly  low  metamor¬ 
phosis  stage.  It  is  very  similar  to  gcbbro  and  some¬ 
times  difficult  to  distinguish  from  it.  Its  mineralogic 
composition  is  actinolite,  and  hornbl ende-pl agiocl ase- 
biotite  (Shimron,  1974).  Texture  is  coarse, 

with  crystals  from  0.5  to  3  mm  in  size.  The  amphibolite 
weathers  more  readily  than  granite,  schist,  or  quartz- 
porphyry;  it  occurs  in  two  terrain  types: 

I)  continuous,  relatively  soft  rock  surfaces,  with  a 
high  degree  of  roughness  and  with  joints  which  cause 
depressed  strips.  This  facies  is  found  on  the  slopes. 

II)  resistant  and  topographically  conspicuous,  jointed, 
rocky  knolls,  located  on  numerous  summit  areas. 

The  "soft"  facies  is  dominant  and,  by  percentage  of 
area  cover,  the  "resistant"  facies  is  nearly  negligible. 
The  massive  appearance  of  the  rock  surfaces  conveys  the 
impression  that  the  amphibolite  is  similar  to  the  granite 
in  infiltration  rate  and  runoff  generation;  however, 
this  impression  is  wrong. 

It  is  clear  from  the  substantial  number  of  tests  carried 
out  in  amphibolite  terrain  (plot  8,  Table  8.6;  plot  20, 
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table  8.7-  8.10;  plot  40,  table  8.11;  and  addi tional  ,  undepen 
dable  plots)  that  this  rock  type  reacts  in  a  manner  in¬ 
termediate  between  rock  and  colluvium.  Its  infiltration 
capacity  is  much  higher  than  that  of  granite  or  dyke 
rocks,  and  is  also  considerably  higher  than  that  of  ve¬ 
ry  jointed  schist. 

The  three  plots  indicate,  as  an  average,  that  runoff  on 
amphibolite  develops  only  after  0.7-0. 8  mm  of  rainfall, 
and  that  the  terminal  infiltration  rate  (standardized 
for  a  20  mm/hr  winter  rainfall)  is  between  5  and  6 
mm/hr . 

After  each  test  on  amphibolite,  a  25  cm  deep  pit  was 
dug  just  outside  each  plot  (photograph  16).  In  all  three 
pits  the  material  was  weathered  rock  in  the  form  of  co¬ 
arse  sand  mixed  with  some  finer  material;  the  largest 
particles,  size  1-4  mm,  were  composed  of  individual 
crystals  or  small  rock  chips.  This"sand"  was  found,  in 
all  3  cases,  to  be  completely  saturated  right  down  to 
the  bottom  of  the  pit.  This  "sand"  possesses  a  very 
high  hydraulic  conductivity,  and  it  seems  that  the  in¬ 
filtration  rate  is  controlled  by  the  less  permeable  top 
layer  and  by  the  unweathered  rock  in  depth. 

Dykes^  Two  perpendicular  systems  of  dykes  criss-cross 
the  entire  area  of  Nahal  Yael  in  a  variable  density. 

They  are  composed  of  different  rock  types,  the  most  im¬ 
portant  being  quartz-porphyry. 

The  dyke  rock  is  very  densely  jointed  and  sometimes 
even  completely  shattered.  For  this  reason  no  dependa¬ 
ble  results  could  be  obtained  in  all  four  tests  conduc¬ 
ted,  as  water  leaked  under  the  plot  boundaries. 

On  the  basis  of  these  rainfall  simulation  test  and  the 
"point"  tests  described  above,  we  assume  that  a  relatively 
small  initial  amount  of  water  (roughly  0.5  mm)  is  nee¬ 
ded  to  satisfy  depression  and  joint  storage;  thereafter 


a  low  infiltration  of  1-2  mm/h.r  applies. 

Colluvi_unK  Extensive  areas  in  the  Nahal  Yael  catchment, 
in  particular  in  its  upper  part,  are  covered  by  a  col¬ 
luvial  veneer.  This  colluvium  is  composed  of  fine  sand 
with  a  large  amount  of  stones  size  20  to  80  mm.  Between 
80  and  100  per  cent  of  its  surface  is  covered  by  angu¬ 
lar  stones  size  10-50  mm.  Observations  following  the 
runoff  event  of  9  February  1979  disclosed  that  during 
natural  rain  a  thin  crust  forms  on  that  part  of  the  sur¬ 
face  not  covered  by  stones.  Hence  it  is  possible  that 
the  simulated  rain,  which  obviously  does  not  permit 
such  a  crust  to  form,  exaggerates  the  natural  infiltra¬ 
tion  rates. 

In  the  tests, runoff  was  initiated  after  3  mm  infiltrated 
The standardi zed  mean  terminal  infiltration  rate  was 
about  9  mm/hr. 

Klein  (1972)  collected  runoff  from  natural  rain  over 
colluvial  plots  in  Nahal  Yael,  but  his  series,  which 
consisted  of  very  small,  1 ow- i ntens i ty  events  with  con¬ 
siderable  interference  from  contributing  dykes  as  well 
as  instrumental  surfaces,  cannot  be  used  for  compari¬ 
son  with  the  present  data. 

The  present  results,  however, agree  well  with  those  of 
Lavee  (1973),  who  applied  simulated  rainfall  to  6  talus 
slopes  in  eastern  Sinai,  three  of  which  (Nos.  1,2,3) 
were  schist  and  granite.  Each  slope  was  subjected  to 
three  spells  of  simulated  rain  totalling  20  mm  over  40 
minutes.  The  mean  terminal  infiltration  rate  for  the 
last  two  spells  on  the  three  above  slopes  was  12.8  mm/hr 
If  standardized  in  the  sense  of  the  present  report,  this 
figure  becomes  10.1  mm/hr,  and  is  comparable  to  the  re¬ 
sults  obtained  for  plots  1  ( tab! e  8.4)and  5  (table  8.5). 
The  difference  of  1.5  mm  might  be  explained  by  the  coar¬ 
ser  material  in  the  Sinai  taluses. 
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For  Lavee’s  third  rain  spell  alone,  the  rate  is  9.4  mm/hr 
(  =  9.1  mm/hr  standardized). 

The  agreement  between  the  two  sets  of  data  is  surpri¬ 
sing,  because  Lavee  simulated  his  rainfall  over  large 
plots  (area  about  100m  ) ,  and  he  stresses  large  diffe¬ 
rence  in  infiltration  rate  within  the  plot  (rills  vs. 
ridges),  while  in  the  present  work  very  small  (0.25m  ) 
homogeneous  micro-plots  were  used. 

8.4  Factors  which  affect  the  infiltration  rate 

all_and_t  ime  _  n  e  §d  ed  _  to  _a  c  h  1  eve  _  term  2.naX_lDZ 
filS^at;ion_rates 

The  following  discussion  and  conclusions  are  based  on 
Table  8.12. 

Rainf §ll_5?!ount^  As  Shown  by  the  two  tests  conducted  on 
granite,  about  1  mm  of  rainfall  is  needed  to  reach  the 
terminal  infiltration  rate.  However,  in  evaluating  the 
effect  of  joints,  a  general  value  of  1.5-2  mm  seems  mo¬ 
re  representative,  because  the  two  plots  tested  had  be¬ 
low  average  jointing. 

The  measured  value  for  schist  -  2.9  mm  -  seems,  on  the 
other  hand,  to  be  on  the  high  side,  as  the  tested  plot 
was  excessively  jointed.  A  reasonable  representati ve  va¬ 
lue  for  schist  slopes  with  outcropping  cleavage  planes 
would  be  2-2.5  mm,  and  less  for  slopes  with  concordant 
schistosity. 

The  mean  value  for  the  amphibolite  plots,  2.3  mm,  being 
the  average  of  3  plots,  seemed  to  be  representative. 

Dependable  trials  for  colluvial  slopes  are  available  at 
only  one  plot,  but  evaluation  of  trials  on  plot  that 
leaked  suggest  that  a  value  of  6  -  1.5  mm  is  a  reasona¬ 
ble  estimate. 


These  results  are  of  considerable  importance  for  the 
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prediction  of  runoff.  Whenever  the  pre-rain  is  fairly 
intense,  i.e.  approaching  the  infiltration  rate,  then 
the  downward  movement  of  water  proceeds  almost  under 
saturated  conditions.  If,  during  such  a  pre-rain,  the 
infiltration  rainfall  amount  suffices  for  achieving  a 
stable  infiltration  rate,  then  only  very  small  addi¬ 
tional  amounts  (tenths  of  mm?)  during  a  very  short  ti¬ 
me  interval  (tens  of  seconds?)  are  needed  to  achieve 
the  terminal  infiltration  rate. 

Time^  The  time  needed  to  achieve  a  stable  infiltration 
rate,  for  various  rainfall  intensities,  is  7  minutes  in 
granite,  10-12  minutes  in  schist  and  amphibolite,  and 
15-20  minutes  in  colluvium.  However,  the  tests  on  two 
of  the  amphibolite  plots  (plot  20,  table  8.10,  and  plot 
40,  table  8.11)  indicate  that,  during  a  rain  duration 
of  15  minutes,  a  stable  infiltration  rate  is  indeed 
achieved,  but  it  :s  not  the  terminal  rate.  Although 
this  could  be  caused  by  inadequacies  in  the  graphic  ap¬ 
proximation  by  the  asymptote  in  relation  to  the  sensivi- 
ty  of  the  measured  data,  a  physical  reason  should  also 
be  considered.  It  is  possible  that  the  wetting  front 
reaches,  as  the  rain  proceeds,  some  different  horizon 
of  lower  hydraulic  conductivity  or  somewhat  different 
moisture  content.  Similar  changes  in  the  "termi nal "  in¬ 
filtration  rate  have  also  been  observed  by  Morin  et  al. 
(1976)  and  Yair  and  Lavee  (1976). 

We  will  assume  that  the  infiltration  rate  of  amphibolite 
becomes  terminal  if  the  pre-rain  lasted  at  least  30  mi¬ 
nutes  and,  of  course,  fell  in  a  sufficient  amount.  For 
all  other  tarrain  types  the  rates  and  times  will  be 
taken  as  stated  above. 

§^4^2 ll£!23ti.gn_ra  te  _  a  s  _  a  unct  i.on_of  _rai  nf  a]_l 

intensity 


Contrary  to  a  decrease  in  terminal  infiltration  rate 
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with  increased  rainfall  intensity,  presumably  because 
of  crusting  on  exposed  soils  (Mol denhauer ,  1960;  Morin 
et  al . ,  1976;  Vaston,  1965),  some  of  the  data  of  this 
study  suggest  a  positive  linear  relationship  between 
these  two  variables  (plots  20  and  40;  tables  8.10, 

8.11). 

The  magnitude  of  this  factor  can  be  illustrated  by  the 
following  example:  on  amphibolite,  if  a  terminal  infil¬ 
tration  rate  of  7.25  mm/hr  applies  for  a  rainfall  in¬ 
tensity  of  12  mm/hr,  then  doubling  the  intensity  to 
24  mm/hr  would  increase  the  terminal  infiltration  rate 
to  8.2  mm/hr.  As  a  first  approximation,  a  linear  rela¬ 
tionship  might  be  developed. 

A  possible  explanation  hypothesizes  lateral  changes 
within  one  plot  of  infiltrating  properties,  similar  to 
those  assumed  by  Shanan  (1975)  to  be  areally  normally 
distri buted . 

As  the  data  base  available  from  this  study  for  the  eva¬ 
luation  of  this  factor  is  too  small,  it  is  not  further 
dicussed  here.  It  was,  however,  taken  into  account  when 
deriving  the  relevant  terminal  infiltration  rates  for 
several  examples  (see  section  8.5  below). 

8.5  Runoff  computed  on  the  basis  of  this  study 
§i:Ll__Methods 

The  factors  which  affect  runoff  are  the  infiltration  rate 
during  the  runoff-producing  rain  spell  and  the  intensi¬ 
ty  of  this  rain.  Also  important  are  rainfall  amount  (and 
intensity)  prior  to  the  effective  rain  spell,  and  ini¬ 
tial  moisture. 

Each  rainstorm  (=  event)  analysed  was  timewise  divided 
as  follows: 

I)  Initial  wetting  rain:  those  parts  of  the  event  pre- 
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ceding  rains  II)  and  III); 

II)  Pre-rain;  a  relatively  intense  rain  which  immedi¬ 
ately  precedes  rain  III); 

III)  Effective  rain:  the  part  of  the  event  which  gene¬ 
rates  runoff; 

IV)  Post-rain:  that  part  of  the  event  which  immediately 
fol 1 ows  rain  III). 

Whenever  the  effective  rain  fell  on  dry  terrain,  the  time 
ofrunoff  initiation  was  calculated  by  dividing  the  amount 
of  water  needed  to  achieve  a  stable  infiltration  rate 
by  the  rain  intensity,  using  the  data  obtained  from  the 
infiltration  tests.  Effective  rain  preceded  by  low- 
intensity  rainfall  presents  some  difficulty,  as  in  all 
our  experiments  the  simulated  rain  applied  was  much  more 
intense  than  the  terminal  infiltration  rate.  The  proce¬ 
dure  adopted  was,  in  cases  where  stages  I)  and  II)  were 
both  present, to  calculate  the  time  of  runoff  initiation 
from  rain  II)  as  above,  and  to  "postpone"  its  time  of 
initiation  by  10  to  30  seconds.  In  cases  of  an  inter¬ 
mission  in  the  rain  between  stages  I)  and  II),  the  cal¬ 
culation  was  done  as  if  the  ground  were  dry,  and  the 
time  of  runoff  "initiation"  was  "postponed"  by  an  arbi¬ 
trary  length  of  time. 

The  terminal  infiltration  rate  was  calculated  by  using 
the  results  of  the  simulated  rainfall  tests,  corrected 
in  relation  to  effective  rain  intensity  as  illustrated 
i n  section  8.4.2. 

The  relatively  short  interval  between  the  time  of  run¬ 
off  initiation  and  the  achievement  of  stable  infiltra¬ 
tion  -  about  8  to  15  minutes  -  exerts  a  confining  influ¬ 
ence  on  the  form  of  the  infiltration  curve.  Although  it 
is  impossible  to  compute  the  infiltration  curves,  using 
the  Horton  formula,  on  the  basis  of  the  three  variables 
alone  (time  of  runoff  initiation;  infiltration  rate  at 
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that  time;  terminal  infiltration  rate),  the  curves  as 
drawn  after  these  data  and  in  accordance  with  the  gene¬ 
ral  form  obtained  during  the  tests,  actual ly  leave  a 
margin  of  error  of  no  more  than  -  0.1  mm.  The  infiltra¬ 
tion  rates  as  used  for  these  evaluations  are  summarized 
in  table  8.14.  The  influence  of  the  post-rain  was  neg¬ 
lected  at  this  stage,  and  will  be  discussed  later. 

Whenever  the  pre-rain  or  the  post-rain  were  intense 
enough  to  yield  runoff  on  non-amphibolite  rocks,  we  in¬ 
cluded  in  the  calculations  one  half  of  the  excess  vo¬ 
lume  of  these  spells,  on  the  assumption  that  such  a 
condition  decreases  the  losses  during  the  effective 
rain  spell. 

The  discharge  data  for  watershed  05  have  been  measured 
over  the  last  13  years.  Because  of  minimal  routing  los¬ 
ses  due  its  small  size  and  rocky  channels,  the  data  for 
this  watershed  provide  an  opportunity  of  comparison 
with  runoff  as  computed  on  the  basis  of  the  infiltra¬ 
tion  data  and  detai 1 ed  in  the  previous  sections. 

A  special  geomorphic  mapping  project  was  undertaken  to 
provide  the  detailed  terrain  information  required  for 
the  development  of  a  deterministic,  hydro-geomorphi c 
dynamic  model  (see  chapter  5  of  this  report).  This  infor¬ 
mation  can  be  adapted  for  use  in  our  comparative  cal¬ 
culation,  as  follows: 

A  R  E  A 


hectares 

per  cent 

Watershed  05,  total 

5.0 

100 

Colluvium,  contiguous  areas 
Colluvium  +  rocky  surfaces 

1.0 

20 

(areal  ratio  2  to  1) 

3.0 

60 

Rocky  surfaces,  contiguous 

1.0 

20 

Table  8.14:  Infiltration  values  for  various  slope  ter¬ 
rains  used  for  calculating  excess  rainfall 


Amphibo¬ 

lite 

Other . 
rocks 

Collu¬ 

vium 

Stable  infiltration  rate 
standardized  to  a  rain¬ 
fall  of  20  mm/hr 

5.66 

1.83 

8.75 

mm/hr 

Duration  of  rainfall  un¬ 
til  achievement  of  sta¬ 
ble  infiltration  rate 

30* 

9 

16 

mi  n 

water  amount  infiltrated 
until  infiltration  rate 
stabilized,  standardized 
to  winter  conditions 

3.1 

1.6 

5.9 

mm 

Water  amount  infiltrated 
until  initiation  of  run¬ 
off 

0.76 

0.6** 

3.05 

mm 

Evaluated  as  an  average 
of  plots  no. 

8,20,40 

7,15,35 

“  1,5 

estimated  on  the  basis  of  the  follow-up  rains  in 
plots  20  and  40. 

average  of  eight  plots,  of  which  only  three  have 
been  reported  here. 

a  finer  distinction  could  be  obtained  by  using  the 
available  data  for  individual  rock-terrain  types 
(see  section  8.2  above);  however,  a  similar  value 
would  have  been  obtained  for  most  rock-terrain  asso¬ 
ciations  which  include  components  of  some  moderately 
cleaved  schist,  some  gneiss-granite  and  some  jointed 
quartz  -porphyry  dyke  material. 


Disaggregating  the  "mixed"  area  and  differentiating  the 
rock  types,  we  get: 

AREA 


hectares 

per  cent 

Colluvium 

2.5 

50 

Amphibolite 

1.75 

35 

Other  rocks 

0.75 

15 

This  is  only  a  gross  estimate  for  a  number  of  reasons: 
the  areal  measurements  were  not  very  accurate  and  the  rock 
characteristics  vary  also  within  the  lithologic  units. 

In  particular,  the  depth  and  composition  of  the  collu¬ 
vium  vary.  Also,  the  spatial  distribution  of  the  litho¬ 
logic  units  is  disregarded,  though  thi'  factor  probab¬ 
ly  exercises  a  decisive  influence  on  runoff  generation. 

§i§.2__Ihe_event_of_25_March_1971 

The  event  of  25  March  1971  is  taken  below  as  a  compu¬ 
tational  example.  This  runoff  event  was  generated  by  a 
uniform,  hi gh- i ntens i ty  rainfall  which  followed  imme¬ 
diately  a  fairly  intense  pre-rain.  It  produced  a  good 
approx ima t i on  of  the  unit  hydrograph  of  the  catchment. 

Rainfall  data: 


Time  from  begin  of 
event  (min  ) 

Rain  stage 

Amount 

(mm) 

Duration 

(min  ) 

Intensity 

(mm/hr) 

0  -  15 

pre-rain 

1 

15 

0-30 

15  -  18 

effective  rain 

3.25 

3 

66 

18  -  31 

post  rain 

0.5 

17 

1.8 

The  intensity  of  the  pre-rain  rose  gradually. 


Runof According  to  our  results,  run¬ 
off  from  amphibolite  is  initiated  after  0.6  mm  infiltra- 
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te  into  the  ground  (this  amount  includes  depression 
storage).  The  duration  until  a  terminal  infiltration 
rate  is  achieved:  10  minutes.  In  the  event  discussed 
the  runof f-produci ng  rain  had  an  intensity  of  66  mm/hr. 
Prior  to  it  already  1  mm  of  rainfall  had  infiltrated, 
and  the  last  part  of  this  pre-rain  was  already  intense. 
Thus  it  is  realistic  to  assume  immediate  runoff  as  soon 
as  the  effective  rain  sets  in,  at  a  rate  controlled  by 
the  stable  infiltration  rate  for  a  15  minute  rain  du¬ 
ration,  corrected  for  the  intensity  of  66  mm/hr.  This 
rate  comes  out  to  be  10.9  mm/hr. 

From  these  three  data  -  initial  losses  0.6  mm;  time  to 
stability  10  min;  stable  infiltration  rate  10.9  mm/hr  - 
a  continuous  infiltration  curve  for  the  event  is  con¬ 
structed  (fig.  8.10).  Although  the  curvature  of  the 
graph  could  be  evaluated  from  the  Horton  equation,  we 
preferred  to  draw  the  curve  in  accordance  with  the 
mean  curvature  obtained  from  the  field  data. 

Measurement  of  the  shaded  area  on  fig.  8.10  results,  by 
the  usual  procedure, in  the  runoff  volume.  In  the  pre¬ 
sent  case  this  value  is  found  to  be  1.25  mm;  multiplied 
by  the  amphibolite  area  of  1.75  ha,  we  get  an  excess 

3 

rainfall  volume  of  22m  . 

Note  that  in  fig.  8.10  the  infiltration  curve  has  been 
shifted  40  seconds  to  the  left,  to  account  for  the  effect 
of  the  pre-rain  which  must  have  caused  the  immediate  ge¬ 
neration  of  runoff  with  the  onset  of  the  effective  rain. 

By D2 f f _f C2'D_2i!]|!I_r29b§i  T*ie  infiltration  curve  used  is 
based  on  the  results  from  three  plots  -  two  on  granite 
and  one  on  schist.  Although  the  non-amphibolite  rock  as¬ 
semblage  in  watershed  05  is  not  identical  to  this  compo¬ 
sition,  it  is  sufficiently  close  to  warrant  this  gene¬ 
ralization. 

Initiation  of  runoff  after  the  infiltration  of  0.6  mm; 


IJ 


Figs.  8.10  -  8.14:  Nahal  Yael  watershed  05:  infiltration  and  excess  rainfall  R, 

events  of  25.3.71,  12.12.77,  and  11.12.78. 


Figs.  8.15 


a  ,,  Nahal  Yael  watershed  05:  hyetograph  (rain  gauges  25,  26,  31 

•  and/or  33,  and  scaled  hydrograph,  events  of  11.12.78  &  9.2.79. 


Time  to  infiltration  stability  -  9  min; 

Stable  infiltration  rate  -  3.25  mm/hr. 

These  data  yield  an  excess  rainfall  of  3.0  mm  over  the 

3 

partial  area  of  0.75  ha,  i.e.  22.5m  excess  rainfall 
volume . 

The  total  rainfall  excess  from  the  rocky  surfaces  of 

3 

watershed  05  was,  for  this  event,  44.5m  . 

By Q9f f _f £°2_£2lly ^ ly The  basic  data  which  govern  run¬ 
off  from  colluvium  were  determined  to  be: 

I)  At  rain  intensity  of  60  mm/hr,  runoff  will  start 
after  3  minutes; 

II)  The  terminal  infiltration  rate  is  9  mm/hr.  Figure 
8.11  shows  that  the  excess  rainfall  is  only  0.1  mm  or, 

3 

over  2.5  ha,  2.5m  . 

In  comparison,  Klein's  (1972)  colluvial  plot  102  yiel¬ 
ded,  for  the  same  event,  at  least  0.75  mm.  Extended  to 
the  2.5  ha  of  colluvium  in  watershed  05,  this  is  equi- 

3 

valent  to  18m  ,  at  least.  The  reason  for  the  discre¬ 
pancy  cannot  be  determined  with  certainty. 

ByD2ff.volume_from_entire_watershed:  Totalling  the  con¬ 
tributions  of  the  lithologic  units,  we  get  a  watershed 
.  3 

rainfall  excess  of  47m  ,  or  -  using  Klein’s  (1972)  data 

3 

-  of  62.5m  .  The  hydrograph  at  station  05  yielded  a 

3 

runoff  volume  of  37.5m  .  Thus,  the  transmission  losses 

3 

were  about  10m  ,  or,  if  we  accept  Klein’s  (1972)  data, 
25m^  or  more. 

.  3__The_event_of  _1 1  _December_197§ 

This  event  consisted  of  a  low  flow  at  stations  05  and 
04.  O.S.  was  present  in  Nahal  Yael  at  the  time  of  the 
event. 

Rainfall^  The  rainfall  data  for  the  event  were  are  de¬ 
tailed  in  table  8.15.  The  intensity  of  the  so-called 
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Table  8. 

15:  Rainfall 

data , 

event  of  11 

December 

1978 

Time  from 
beginning 
of  event 

Type 
(phase) 
of  rain 

Amount 

Duration 

Intensity 

Cumul . 
amount 

(hr:min) 

(mm) 

(min) 

(mm/hr) 

(mm) 

0 

0:00 

pre-rain  A 

4.5 

2.2 

4.5 

2:20 

pre-rain  B 

1.8 

19.5 

5.5 

6.3 

2:40 

medium 

fee- 

0.6 

2.25 

15.9 

6.9 

2:42 

Intense  t1ve 

2.3 

4.5 

30.6 

7.6 

2:47 

post-rain 

0.7 

35 

1.3 

8.3 
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wetting  rain,  which  lasted  over  two  hours,  is  only  an 
average  value,  representing  an  actually  variable  inten¬ 
sity.  Likewise,  the  intensity  of  the  post-rain  is  also 
an  average  value.  . 

§§2idGlHS_of_!iy Q2tf  1  Because  of  the  relatively  Intense  pre¬ 
rain  (B),  slope  runoff  was  initiated  immediately  with 

i 

the  onset  of  the  effective  rainfall,  and  it  stopped  right 
at  the  end  of  the  effective  rainfall  period. 

The  terminal  infiltration  rates  for  a  rain  lasting  se¬ 
veral  hours  and  corrected  to  an  intensity  of  30  mm/hr 
are: 

for  amphibolite  6.5  mm/hr; 

for  other  rocks  2.1  mm/hr. 

As  pre-rain  B  was  fairly  intense,  and  the  entire  pre¬ 
rainfall  totalled  6.3  mm,  the  terminal  rate  applies, 
approximately,  to  all  rock  types  from  the  onset  of  the 
medium-effective  rain. 

Excess  rainfall,  amphibolite:  2.15  mm  x  1.75ha  *  37. 6m^; 

3 

Excess  rainfall,  other  rocks:  2.6  mm  x  0.75  ha  =  16m  . 

3 

Total  excess  rainfall  on  non-colluvial  surfaces:  54m  . 

The  estimation  of  the  infiltration  rate  into  the  collu¬ 
vium  is  especially  difficult,  because  it  cannot  be  as¬ 
signed  the  terminal  rate  from  the  onset  of  the  effective 
rain  stage.  No  data  are  available  to  enable  the  assess¬ 
ment  of  the  effect  of  the  pre-rains  on  the  infiltration 
rate.  In  practice,  however,  several  reasonable  assump¬ 
tions  which  bridge  this  gap  of  information  yield  quan¬ 
titatively  similar  results.  We  assumed,  arbitrarily 
that  the  medium-effective  rain  spell  is  needed  to  initi¬ 
ate  runoff,  and  thereafter,  only  one  half  of  the  dura¬ 
tion  needed  in  fig.  8.11  is  necessary  to  achieve  sta¬ 
bility;  i.e.  5  minutes.  The  resulting  curve  (fig.  8.12) 
yields  an  excess  rainfall  of  1.2  mm,  equivalent  for  the 

3 

colluvial  area  of  2.5  ha  to  30m  . 


In  case  a  crust  on  the  colluvium  begins  to  develop,  due 
to  the  relatively  long  rainfall  duration,  this  excess 
rainfall  might  become  larger. 

The  non-amphi bol i tic  rocks  yielded  some  runoff  during 

the  pre-rain  stages.  The  added  excess  rainfall  from  this 

3 

sources  is  estimated  to  be  3.5m  . 

3 

Totalling  the  excess  rainfall,  we  get  88m  ,  in  relation 

3 

to  the  runoff  volume  of  51m  measured  at  the  gauging 
station  of  watershed  05;  i.e.  the  transmission  losses 

3 

were  33m  . 

Expressed  in  a  different  way,  the  excess  rainfall  for 
this  event  constitutes  61%  of  the  runof f-produci ng  rain, 
but  the  watershed  yield  is  only  35%. 

§i:L4__The_event_gf_9_Februar¥_i979 

Rainfalls  The  rainfall  data  listed  in  table  8.16  are 
mean  values  for  rain  recorders  26  and  33.  A  rainless  in¬ 
terval  of  one  hour  followed  the  initial  wetting  phase. 
Evaluation^  The  pre-rain  did  not  yield  any  appreciable 
runoff.  In  order  to  reach  stable  infiltration  rate, 

2.3  mm  of  water  must  infiltrate  into  the  amphibolite  as 
saturated  flow  (Table  8.14).  The  intensity  of  the  pre¬ 
rain  exceeded  slightly  that  stable  rate.  This  means  that 
infiltration  proceeded  as  unsaturated  flow  at  the  begin¬ 
ning  of  the  pre-rain  phase,  and  approached  saturation 
conditions  towards  its  end.  Therefore  it  is  permissible 
to  assign  terminal  infiltration  rates  to  the  amphibolite 
terrains  (as  well  as  to  the  other  rocky  surfaces)  as  of 
the  onset  of  the  effective  rain,  or  very  shortly  there¬ 
after. 

The  stable  infiltration  rate  for  amphibolite  was  deter¬ 
mined  to  have  been  5.2  mm/hr,  and  for  the  other  rocky 
surfaces  -  1.6  mm/hr.  A  correction  to  an  intensity  of 
16  mm/hr  was  taken  into  account  in  the  deviation  of  these 
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Table  8.16:  Rainfall  data,  event  of  9  February  1979 


Time  from 
beginning 
of  event 

Type 
(phase) 
of  rain 

Amount 

Duration 

Intensity 

Cumul . 
amount 

(hrrmin) 

(mm) 

(min) 

(mm/hr) 

(mm) 

0 

0:00  -  0:40 

initial 

wetting 

1.6 

40 

2.4 

1.6 

1:40  -  2:05 

pre-rain 

2.9 

23 

7.8 

4.5 

2:05  -  2:11 

effective 

1.6 

6 

16 

6.1 

2:11  -  2:25 

post-rain 

0.4 

15 

1.8 

6.5 

2:25  -  2:33 

post-rain 

0.5 

8 

3.5 

7.0 
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values.  The  resulting  excess  rainfall  for  the  non- 

3 

colluvial  parts  is  29.7m  . 

The  colluvial  parts  present  some  complications.  Under 
saturation,  5.9  mm  needed  to  reach  a. stable  infiltration 
rate  (Table  8.14).  The  pre-rain  supplied  one  half  of 
this  amount,  and  even  if  the  entire  wetting  rain  is 
added,  we  are  still  1.4  mm  short.  Therefore  it  is  impos¬ 
sible  to  determine  neither  an  infiltration  rate  nor  an 
infiltration  curve  for  the  effective  rain  spell,  be¬ 
cause  the  range  of  the  reasonable  assumptions  is  too 
wide  and  may  cause  large  errors  . 

This  event  is  quite  similar  in  character  to  the  one  des¬ 
cribed  in  the  previous  section  (II  December  1978). 

Hence  we  assume  also  a  similar  response  of  the  colluvi¬ 
um,  and  estimate  that  it  yielded  roughly  one  half,  or 
slightly  less,  of  the  water  yield  of  the  rocky  terrains, 

3 

i.e.  about  10-15m  . 

Since  the  effective  rain  of  this  event  fell  only  at  one 
half  of  the  intensity  of  the  previous  event  (11  Decem¬ 
ber  1978),  it  is  possible  that  the  water  yield  off  the 
colluvium  is  even  less.  However,  this  effect  may  be 
compensated  by  some  crusting  on  the  small,  inter-stone 
surfaces  which  had  been  preserved  from  the  event  which 
occured  four  weeks  earlier. 

All  these  factors  considered,  we  arrive  at  a  total  ex- 

3 

cess  rainfall  for  watershed  05  of  50-55m  .  Station  05 

3  3 

recorded  31m  ,  which  leaves  20-25m  for  transmission 
losses. 

The  rainfall  excess,  including  parts  of  the  pre-rain 
and  the  post-rain,  constitute  62-6925  of  the  effective 
rainfall;  the  watershed  runoff  volume  amounts  to  39% 
thereof. 
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inf  a.ii_§v  §nt_and_  indications  _of_runoff_thre; 
§112.1^1  _I?_December^  197  7 

During  this  event  no  flow  reached  gauging  station  05, 
but  the  channel  immediately  upstream  of  the  weir  had 
low  flow  marks  along  a  reach  of  approximately  50  meters. 

The  rainfall  data  for  the  event,  as  evaluated  from  rain- 
gauges  26  and  33,  were: 


Type 

Amount,  mm 

Duration,  min 

Mean  int.,  mm/hr 

Effective 

2 

2 

60 

Medium 

1 

4.5 

15 

Weak 

0.5 

7.5 

4 

Figures  8.13  and  8.14  relate  the  hyetog^aph  to  the  mean 
infiltration  curves  of  amphibolite  and  of  other  rocks, 
respectively.  Note  that  the  "medium"  rain  phase  actu¬ 
ally  had  a  gradually  decreasing  intensity. 

On  colluvial  surfaces  no  runoff  was  generated.  Excess 

3 

rainfall  from  amphibolite  amounted  to  4.9m  ,  and  from 

3 

other  rocks  -  6.2m  .  Most  of  the  total  rain  excess* 

3 

which  amounted  to  10m  ,  was  absorbed  by  the  colluvium. 
Small  residues  were  absorbed  by  several  minute  allu¬ 
vial  pockets  in  the  channel,  in  conjunction  with  the 
very  low  flow  mentioned  above. 

§i§i§__6§l!2f5ll_§vent,_6_February_1978 

On  6  February  1978  a  rain  event  of  eight  hours' inter¬ 
mittent  duration  occured  in  Nahal  Yael .  The  total  rain¬ 
fall  of  8  mm  was  insufficient  to  generate  runoff,  nor 
were  low  flow  marks  evident.  Thus  it  was  below  the  thre¬ 
shold  of  runoff. 

The  rainfall  data  for  the  early  and  relevant  part  of  the 
event,  based  as  usual  on  raingauges  26  and  33,  were: 
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Time  from 
start  of 
event 
hrsrmin 

Type  of 
rain 

Amount 

mm 

Duration 

min 

Intensity 

mm/hr 

Cutnul. 

Amount 

mm 

0:00 

0 

wetting 

2.6 

1:40 

2.6 

intermission 

0 

2:10 

2.6 

intense 

1.4 

2 

42 

2:12 

4.0 

post-rain 

0.3 

5 

3.6 

2:17 

4.3 

The  wetting  phase  actually  had  quite  variable  intensi¬ 
ties  with  some  intervening  rainless  spells.  The  post¬ 
rain  actually  had  a  gradually  declining  intensity.  The 
rainfall  amount  was  also  spatially  fairly  variable,  with 
4.3  mm  being  a  mean  value. 

Assuming  that  the  rain  fell  on  a  completely  dry  surface, 
the  following  figures  result: 

Amphibolite:  Runoff  initiated  65  seconds  after  onset  of 

3 

"strong"  rain.  Excess  rainfall  0.12  mm  or  2.1m  . 

Other  rocks:  Runoff  initiated  50  seconds  after  onset  of 

3 

"strong"  rain.  Excess  rainfall  0.17  mm  or  1.2m  . 
Colluvium:  No  runoff. 

3 

Total  excess  rainfall:  3.3m  . 

Some  small  amount  must  be  added  to  this  figure,  to  ac¬ 
count  for  the  possible  contribution  of  the  post-rain, 
facilitated,  perhaps,  by  the  cumulative  effect  of  the 
long  wetting  phase.  However,  even  the  doubling  of  the 
above  excess  rainfall  will  not  generate  any  channel 
f  1  ow. 


8.  Ek7 _ Summary 

Table  8.17  lists,  for  the  Nahal  Yael  watershed  05,  the 
volumes  of  excess  rainfall  for  five  events,  as  computed 
on  the  basis  of  the  present  study,  and  relates  them  to 
the  measured  runoff  volume  at  the  gauging  station.  For 
the  three  runoff  events,  the  transmission  losses  amount 
to  between  10  and  37m  per  event.  The  recipients  of  this 
water  are:  I)  storage  in  the  few  small  patches  of  allu¬ 
vium  in  the  channel;  II)  filling  of  tiny  rock  pools  in 
the  main  channel  bed  and  in  some  tributaries;  and  III) 
infiltration  into  those  parts  of  slope  terrain  charac¬ 
terized  by  excessive  infiltration  capacity  (rill  zones 
and  negative  dykes  filled  with  coarse  colluvium). 

Theoretically,  losses  of  types  I)  and  II)  should  be  con¬ 
stant,  while  those  of  type  III)  should  relate  in  some 
consistent  manner  to  the  runoff  which  actually  occurs 
on  the  slopes.  From  our  data,  the  upper  limit  of  losses 
is  37m3,  and  the  lower  limit,  about  10m3.  Further  stu¬ 
dies,  including  perhaps  some  additional  measurements,  are 
needed  to  verify  the  above  reasoning  and  to  elucidate 
its  quantitative  relationships. 

As  a  summary  estimate,  the  mean  transmission  losses  for 

3 

a  one  hour  flood,  volume  30-60m  ,  in  watershed  05,  ge¬ 
nerated  by  an  effective  rainfall  duration  of  5-10  minu- 

3 

tes,  are  approximately  25m  . 
§.::L§__Il3§-§ff®£t_of_the_post-rain 

The  post-rain  was,  generally,  disregarded  in  our  compu¬ 
tations  of  excess  rainfall,  because  its  intensity  is 
lower  than  the  terminal  infiltration  rate;  hence  it  can¬ 
not  contribute  runoff.  However,  in  a  catchment  of  the 
size  of  Nahal  Yael  05,  the  concentration  time  of  excess  rain¬ 
fall  from  the  slopes  to  the  mouth  amounts  to  several  tens  of 
minutes.  Therefore,  in  certain  rainstorm  patterns,  some 
post-rain  may  fall  directly  on  the  remnants  of  the  water 
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Table  8.17:  Watershed  05  -  excess  rainfall  and  runoff 
volume 


Date 

Duration  of 

effective 

rainfall 

(min) 

Computed 

rainfall 

excess 

(m3) 

Measured 

runoff 

volume 

(m3) 

Difference: 

transmission 

losses 

(m3) 

25.3.71 

3 

47 

37.5 

9.5 

(62.5) 

(25) 

11.2.78 

7 

88 

51 

37 

9.2.79 

6 

50  -  55* 

31 

19  -  24* 

12.12.77 

- 

10 

0** 

10 

6.2.78 

- 

5-10 

0*** 

5-10 

Mean 

value  per  runoff  event 

23 

(  )  incorporating  the  rainfall  excess  from  colluvium 

terrains  as  determined  by  Klein  (1972); 

ft 

ftft 


A  A 


low  accuracy,  inherent  in  method  used; 
marks  of  a  very  low  flow  found  in  channel  bed; 
no  marks  of  flow  found. 
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"sheets"  which  flow  down  the  slopes.  In  this  situation 
these  raindrops  of  the  post-rain  join  the  body  of  flow¬ 
ing  water  and  thus  are  directly  contributed  to  the  run¬ 
off,  in  analogy  to  direct  rain  which  falls  on  a  flow¬ 
ing  river.  A  quantitative  estimate  of  this  contribution 
is  attempted  below. 

Figure  8.16  shows  that  the  post-rain  of  the  event  of 
11  December  1978  cannot  affect  the  channel  flow.  The 
hyetographs  reveal  that  the  part  of  the  post-rain  con¬ 
tiguous  with  the  effective  rain  had  an  intensity  of 
0.66  mm/hr.  During  this  rain,  a  flow  started  at  the  05 

3 

gauging  station  at  a  discharge  of  0.038m  /sec;  then  it 
gradually  receded.  At  17  minutes  after  the  beginning  of 
flow  at  the  station,  the  rainfall  intensity  rose  to 
4.3  mm/hr,  but  without  any  appreciable  effect  on  the  hy¬ 
drograph,  which  continued  its  recession  much  like  in 
the  event  of  25  March  1971  (the  "unit  hydrograph"  event) 
until  its  extinction  approximately  one  hour  after  its 
begi nni ng . 

In  the  event  of  9  February  1979  the  intensity  of  the 
post-rain  was  higher  -  1.8  mm/hr  over  15  minutes,  later 
increasing  to  3.5  mm/hr.  Here  too  a  possible  effect  on 
the  hydrograph  is  not  discernible  (fig.  8.17). 

A  post-rain,  mean  intensity  1.8  mm/hr,  followed  also  the 
effective  spell  of  the  storm  of  25  March  1971.  Its  quan¬ 
titative  effect  will  be  evaluated  below. 

Rainfall  data,  post-rain:  amount  -  0.5  mm;  duration  - 
17  min;  mean  intensity  -  1.8  .„m/hr. 

Assumptions:  I)  for  the  first  7  minutes  -  5%  of  the 
catchment  are  covered  by  flowing  water;  II)  for  the 
subsequent  10  minutes  -  3 %  of  the  catchment  are  cove¬ 
red  by  flowing  water. 

Result:  The  total  volume  contributed  by  the  post-rain 
to  these  proportional  areas  is  1.7m  -  a  negligible 

amount . 
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Similary  computed,  the  post-rain  contribution  to  the 

3 

event  of  11  December  1978  was  1.0m  and  to  that  of 

3 

9  February  1979  -  approximately  1.5m  (over  the  first 
15  mi nutes ) . 

The  addition  of  the  potential  excess  rainfall  from  those 
small  areas  whose  terminal  infiltration  rate  is  less 
than  the  relatively  low  intensities  of  the  post-rain 
likewise  adds  only  a  negligible  amount. 

Therefore  in  calculating  excess  rainfalls  for  this  study, 
the  contribution  of  the  post-rain  was  disregarded. 

§^5^9 _ §®^lment 

The  amount  of  sediment  in  the  runoff  collected  at  the 
outlet  funnel  of  the  plots  during  the  infiltration 
tests  described  above  defies  simple  analysis,  because 
of  its  dependence  on  drop  size  and  fall  velocity  through 
splash.  It  is  possiole,  though,  to  meaningfully  evaluate 
the  variation  in  sediment  concentration  over  time  in  a  single 
plot  during  a  simulated  rainstorm.  A  preliminary  evalu¬ 
ation  of  this  kind  is  attempted  on  the  basis  of  the  da¬ 
ta  for  plot  20,  on  which  four  successive  storms  had 
been  rained  over  a  time  span  of  24  hours  (Table  8.18). 

Within  each  of  the  sub-storms  .produced  by  simulated  rain 
lasting  10-15  minutes  each,  no  systematic  trend  in  sedi¬ 
ment  concentration  is  evident.  The  somewhat  lower  con¬ 
centration  values  in  the  samples  of  the  two  last  sub¬ 
storm  in  relation  to  the  two  first  ones  is  probably  cau¬ 
sed  by  lessened  availability  of  suitable  particles  due 
to  a  "cleaning"  effect  of  the  first  flows. 

Contrary  to  expectations,  the  concentration  values  for 
the  first,  high-intensity  sub-storm  are  not  significant¬ 
ly  higher  than  those  for  the  second,  low-intensity  storm. 
No  explanation  for  this  surprising  phenomenon  is  offerred 
at  this  stage.  The  sediment  aspect  of  this  study  will  be 


Tine,  in  minutes 
after  begin  of 
simulated  rain 

1 

1521 

2 

875 

479 

3 

1648 

916 

626 

4 
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512 
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5 

• 
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1312 

6 

1213 

1193 
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7 

1120 

451 
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521 

\ 

8 

1078 
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1040 

9 

1150 

740 

10 

734 

1 

11 

j! 
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4 

12 

i 

1 

13 
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-1 14- 


treat  ed  in  detail  and  in  a  wider  context  in  future. 
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9.  THE  AREAL  FRAMEWORK  OF  THE 
NAHAL  YAEL  MODEL 


by  Gabriel  Sharon  and  Asher  P.  Schick 


9.1  Introduction 

This  chapter  summarizes  a  preliminary  effort  to  de¬ 
velop  an  areal  framework  for  a  distributed,  geosystem- 
oriented  model  for  a  watershed,  based  on  the  example  of 
Nahal  Yael .  In  this  paragraph  the  aim  and  general  struc¬ 
ture  of  the  model  are  described,  as  well  as  the  initial 
stages  of  its  development.  Subsequent  sections  deal 
with  descriptive  aspects  of  the  model  and  explain  how  it 
is  set  to  accept  information  (9.2.),  and  with  a  theore¬ 
tical-mathematical  description  of  the  method  adopted 
for  the  simulation  of  the  terrain  topography  and  of  its 
drainage  network  on  the  basis  of  elevation  data  for  grid 
points  (9.3.  and  9.4.).  Section  9.5.  discusses  the  so¬ 
lutions  described  in  9.3.  and  9.4.  as  treated  by  compu¬ 
ter  programs. 

For  the  sample  input  to  these  computer  programs  we  have 
selected  the  uppermost  part  of  the  Nahal  Yael  watershed, 
with  the  incorporation  of  the  entire  watershed  remai¬ 
ning  to  be  done  during  a  future  stage  of  the  work.  The 
methods  of  computation  and  of  data  input  into  the  com¬ 
puter  programs  as  well  as  an  analysis  of  the  results  are 
presented  in  section  9.5. 

In  addition,  all  along  this  chapter  we  have  dwelt  on  ge¬ 
neral  trends  and  evaluations  as  they  gradually  evolved 
during  the  development  of  the  model  framework.  Clearly, 
this  is  only  an  initial  step  in  the  development  of  the 
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model .  Therefore  we  have  freely  included  ideas  and  sug¬ 
gestions,  to  be  worked  out  in  detail  during  the  future 
stages  of  the  project. 

Ib§_2^i§S5lY§_2f  _  t  h  e  _n?o  d  e  1 

The  Nahal  Yael  model  is  designed  as  a  distributed  model, 
which  should  ultimately  develop  into  a  vehicle  which 
will  evaluate  not  only  runoff  but  also  erosion-sedimen¬ 
tation  processes  as  well  as  effects  thereon  of  short  and 
long  term  pulses  such  as  high-magnitude  events,  landuse 
changes,  climatic  shifts,  and  man-made  intervention.  The 
model  is  constructed  in  such  a  way  as  to  accommodate  the 
fact  that  variations  both  in  time  and  in  space  of  the 
input  data  substantially  affect  runoff  and  associated 
output  data. 

The  areal  framework  chosen  for  the  model  is  one  of  rela¬ 
tively  small  squares  equal  in  area.  It  is  assumed  that 
within  each  grid  element  (=  "square")  the  geophysical  vari 
ables  are  homogeneous.  The  parameters  vary  from  one  square 
to  the  next.  The  dominant  processes  in  the  watershed 
occur  in  two  phases:  the  diffuse  slopeflow  phase,  regar¬ 
ded  as  taking  places  over  the  entire  surface  of  the  squa¬ 
res,  and  the  concentrated  channelflow  phase,  simulated 
by  a  trajectory  composed  of  square  sides. 

The  objective  is  to  construct  a  general  framework  for 
the  quantitative  description  of  attributes  and  processes 
within  a  watershed  in  such  a  way,  that  it  will  be  adaptable 
with  little  change  and  only  nr'nor  preparatory  effort  to 
other  watersheds  and  regions.  Thus  it  is  necessary  to 
divide  the  area  very  mechanically  into  geometric  elements 
(such  as  squares);  a  subjective  areal  division,  such  as 
one  based  on  geomorphic  or  other  terrain  criteria,  is  to¬ 
tally  unsuitable  for  tie  universality  requirement. 

It  is  very  important  to  use  in  the  model  easily  obtai¬ 
nable  input  data,  of  such  kinds  not  requiring  elaborate 


measurements  or  special  research.  The  information  avai¬ 
lable  is  stored  in  relation  to  the  individual  grid  ele¬ 
ments.  Likewise,  the  output  -  water  depth  and  yield,  se¬ 
diment  depth  and  yield  -  will  be  obtained  discretely 
for  each  individual  grid  element  or  channel  element  and 
for  every  time  unit,  as  required.  Thus,  the  model  we  aim 
at  is  a  distributed  deterministic  model. 

Upon  the  completion  of  the  model,  it  should  be  possible, 
in  addition  to  the  uses  already  mentioned,  to  follow  in 
detail  the  evolution  of  landforms  in  diverse  regions. 

The  evaluation  of  the  effects  of  changes  which  occur  in 
one  area  within  the  watershed  on  another  area  in  it  will 
also  be  much  simpler. 

The  numerous  iterative,  though  mostly  simple,  computa¬ 
tions  which  this  model  presupposes  are  possible  only 
with  the  electronic  computer.  Such  computations  will  be 
done  separately  for  each  grid  element,  and  will  be  exe¬ 
cuted  in  addition  to  other  computational  stages  common 
to  models  such  as  optimization  of  parameters. 

This  section  describes  only  the  first  stage  accomplished 
-  the  adoption  of  the  grid  network  for  the  areal  frame¬ 
work  of  the  model  and  the  construction  thereon  of  the 
simulated  topography  and  drainage  system.  To  this  frame¬ 
work  it  is  possible  to  attach  the  input  data. 

9.2  The  Rectangular  Grid 

As  it  is  impossible  to  construct  a  model  which  will  han¬ 
dle  differently  every  single  point  within  an  area  under 
study,  we  have  to  assume  homogeneity  of  geophysical  attri¬ 
butes  within  small  areal  elements.  In  particular,  we  as¬ 
sume  homogeneity  in  the  slope  angle  and  slope  orientation 
of  each  areal  element,  i.e.  we  are  dealing  with  an  assem- 
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blage  of  planar  surfaces. 

Because  of  the  requisite  of  an  arbitrary  division  Into 
unit  elements  mentioned  above,  we  need  some  perfect 
polygon  for  the  unit  area.  We  selected  the  square ,  with 
the  grid  predetermined  by  the  map  coordinate  system,  and 
the  cardinal  N-S  and  E-W  coordinate  lines  coinciding  with 
grid  lines.  This  facilitated  easy  reading  of  data  off 
the  map  for  assigning  to  the  squares.  The  parallel  align¬ 
ment  of  all  the  sides  of  all  the  squares  constitutes  a 
decided  advantage  in  comparison  with  other  perfect  poly¬ 
gons  such  as  the  triangle  or  the  hexagon;  this  is  espe¬ 
cially  felt  in  the  preparation  of  the  data  for  the  compu¬ 
ter  and  in  the  ease  of  reading  the  evaluations  off  the  map 
at  grid  intersections. 

It  might  well  be  that  the  equilateral  triangle  _  and  its 
analog,  the  equilateral  hexagon,  both  constitute  a  more 
"natural"  unit  areal  element,  and  perhaps  they  provide 
more  directional  flexibility  in  areal  transfers  of  mass 
from  one  grid  element  to  the  next.  But  the  i nconveni ence 
of  handling  and  managing  non-cartesian  elements,  coupled 
with  the  more  involved  mathematics  necessary  for  simulating 
mass  transfers,  by  far  outweigh  the  advantages  to  be 
gained.  In  any  case,  if  the  unit  area  is  small  enough,  and 
the  transfer  equations  are  handled  properly,  the  exact 
form  of  the  unit  area  is  of  no  consequence,  particulary 
if  we  recall  that  every  element,  however  small,  is  only 
an  approximation  of  the  real-world  terrain  surface. 

?i?il__Q2D§lderati onsf or_the_sel ecti gn_of _the_sguare_s1ze 

In  principle,  the  smaller  the  unit  area,  the  larger  will 
be  the  degree  of  homogeneity  within  each  square.  Therefore 
the  tendency  in  distributed,  deterministic  modelling  is 
to  have  unit  areas  as  small  as  possible.  However,  there 
are  a  number  of  constraints:  size  of  the  area  (watershed) 
under  study,  degree  of  detail  of  the  information 


-120- 


(geophysical  data),  and  computer  capacity  and  time  avai¬ 
lable. 

For  very  large  areas  the  degree  of  detail  will  be  small, 
and  the  division  into  tiny  units  cannot  be  justified,  be¬ 
cause  the  rate  of  change  of  the  geophysical  parameters 
over  adjoining  unit  areas  will  be  very  small.  On  the 
other  hand,  the  square  side  will  always  have  to  be  at 
least  one  order  of  magnitude  smaller  than  Schumm's  con¬ 
stant  of  channel  maintenance  (=  the  inverse  of  drainage 
density),  in  order  to  ensure  that  the  mean  slopeflow 
trajectory  from  divide  to  channel  is  covered  by  at  least 
five  unit  elements. 

In  our  sample  area  -  described  in  more  detail  in  section 

9.5  below  -  we  have  elected  to  use  unit  elements  of 

12.5  meters  square,  delineated  on  a  detailed,  photogram- 
metry-based  topographic  map.  The  scale  of  the  map  is 
1:1250  and  the  contour  interval  is  one  meter.  At  the  pre¬ 
sent  stage,  the  information  input  for  each  square  con¬ 
sisted  only  of  the  elevation  of  each  corner  of  each  square 
(=  every  intersection  for  80  grid  lines  per  kilometer  of 
N-S  and  E-W  length).  Thus,  neither  the  density  of  geophy¬ 
sical  information  nor  considerations  of  computer  time 

and  storage  had  any  influence  on  the  size  selection. 
However,  if  both  the  size  of  the  study  area  and  the  num¬ 
ber  of  geophysical  parameters  increase  substantially 
during  future  stages  of  the  development  of  the  model  - 
the  size  selected  for  the  unit  area  will  have  to  be  re¬ 
cons  i dered . 

Although  we  have  divided  the  terrain  continuum  in  the 
present  work  into  equal-area  elements,  it  is  quite  possible 
that  in  the  course  of  further  development  the  need  for 
unit  area  flexibility  will  arise.  This  could  be  attained 
by  starting  with  an  equal-area  basic  grid,  but  certain 
squares  or  group  of  squares  could  be  further  subdivided 
on  the  basis  of  an  objective  criterion  reflecting  the 


rate  of  change  over  adjacent  squares  of  a  given  geophy¬ 
sical  property.  For  example,  let  us  look  at  at  topogra¬ 
phic  relief,  i.e.  rate  of  change  of  elevation  over  dis¬ 
tance.  We  might  wish  to  subdivide  the  standard  square 
over  areas  of  cliffs,  canyons  or  badlands. 

Assume  hj,  h^  and  h^  to  be  the  elevations  of  the  cor¬ 
ners  of  a  standard  square,  and  hm  to  be  the  elevation  of 
the  square  center  (fig.  9.1).  Read  from  the  map,  or 
from  the  air  photograph,  those  five  values  and  compute 

4  2 

a  =  E  (h.  -hm)*  (9.2.1) 

i  =  1 

for  all  the  standard  squares.  Squares  with  a  large  a 
will  be  those  with  a  considerable  inner  relief,  and  those 
may  be  divided  into  four  small  squares  each  (fig.  9.1). 

In  practice  a  threshold  value  a  -  dependent  on  the  fre¬ 
quency  distribution  of  the  computed  a  values  -  could  be 
determined,  and  a  further  subdivision  made  wherever 
a  <  a . 

The  application  of  a  criterion  of  this  type  enables  us 
to  choose  a  relatively  large  standard  square,  and  at  the 
same  time  to  augment  the  sensitivity  of  the  model  to 
changes  in  the  geophysical  input.  However,  the  resulting 
complications  in  the  mathematical  representation  of  the 
mass  transfer  from  square  to  square  (not  considered  in 
the  context  of  this  chapter)  may  offset  the  advantages. 

9^2^2 _ Mode_of_simulati^on_with_regard_to_the_grid 

The  simulation  of  process  in  the  studied  area  will  be 
tied  to  the  grid  as  follows:  the  square  areas  will  accom¬ 
modate  the  assumedly  unconcentrated  slopeflow,  and  some 
of  the  square  sides  will  carry  along  them  the  disconti¬ 
nuities  of  concentrated  channel  flow,  as  well  as  those 
of  the  divides. 
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All  parameters,  such  as  elevation,  terrain  composition, 
rainfall  input,  and  others,  will  be  referenced  to  the 
grid  elements.  For  static  parameters  this  will  be  con¬ 
stant  for  all  time  units,  and  for  dynamic  parameters 
they  will  change  with  time.  The  same  will  be  done  for 
those  square  sides  which  simulate  channels. 

Flow  processes  will  occur  diffusely  from  square  to  squa¬ 
re,  until  they  hit  a  channel  side.  From  then  on  flow 
will  proceed  along  the  channel  network  simulated  by  a 
connective  assemblage  of  square  sides. 

In  this  manner  another  constraint  for  the  size  of  the 
square  would  be  the  actual  width  of  the  channel.  It  is 
undesirable  in  this  model  for  a  real-world  channel  to  be 
wider  than  a  square  side. 

9.3  The  description  of  the  topography  by  means  of  the 
grid  data 

From  our  basic  premise  according  to  which  every  square 
unit  element  approximates  a  plane  surface,  there  follows 
immediately  the  first  step  of  computing  the  orientation 
and  slope  angle  of  that  surface  on  the  basis  of  the  ele¬ 
vation  of  its  corners.  The  problem  is  analogous  to  the 
computation  of  dip  and  strike  values  for  a  geological 
stratum  on  the  basis  of  four  outcrops  whose  location  in 
space  is  known. 

The  orientation,  which  will  be  also  named  below  the  azi  - 
mu  th ,  is  the  direction,  in  spa^e,  of  the  maximum  slope  of 
the  plane  surface.  Every  square  can  thus  receive  an  azi¬ 
muth  value  between  0° ( nor thward )  through  90°(eastward) 
and  further  clockwise  all  the  way  back  to  360°(again 
northward,  i.e.  =  0°).  Clearly  the  azimuths  are  orthogo¬ 
nal  to  the  terrain  contours.  The  slope  angle,  subsequent¬ 
ly  also  called  the  dip,  is  the  rate  of  elevation  loss  of  the 
plane  surface  with  regard  to  the  horizontal  xy  plane  of 
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the  grid. 

Three  points  suffice  to  unilaterally  define  a  plane  sur¬ 
face  in  space.  Since  we  have  four  points,  the  computa¬ 
tion  of  both  the  azimuth  and  the  slope  angle  was  made 
separately  for  every  combination  of  three  out  of  these 
four  points.  The  azimuth  and  dip  values  assigned  to  the 
square  are  the  mean  of  these  four  combinations  (details 
below). 

Hence  the  computational  procedure  consists  of  two  sta¬ 
ges  : 

I)  Computation  of  the  azimuth  and  the  dip  for  three 
given  points;  and 

II)  Estimation  of  the  azimuth  and  the  dip  for  the  grid 
square . 

9i3il__Computation_of_a2imuth_and_dip_for_three_giv|nI 

P22!]£§_i!)_§!?§£§ 

The  computation  is  based  on  considerations  of  analytical 
geometry  in  three  dimensions  as  follows:  there  are  given 
three  different  points  in  the  XYZ  space.  XY  is  the  map 
plane,  with  the  x-axis  directed  east-west  and  the  y-axis 
north-south;  z  is  the  axis  of  topographic  elevation. 

Let  us  assign  the  following  coordinates  to  the  three 
points: 

A  =  (Xj,  yj,  Zj) 

B  -  (^2*  y2,  z 2 ) 

C  =  ( x3  ’  ^3  »  z3 )  (fi9-  9-2) 

The  equation  of  the  unilaterally  defined  plane  which  in¬ 
cludes  A,  B  and  C  is,  in  determinant  notation 
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x  y  z 


X1  yl  Z1 


x2 

x3 


y2 

y3 


z2 

z3 


1 

1 

1 

1 


=  0 


(9.3.1) 


All  the  points  (x,y,z)  for  which  equation  (9.3.1)  holds 
are  located  on  the  plane  defined  by  ABC.  Evaluating  the 
determinant  after  the  first  line: 


*1  Z1  1 

xi  zi 1 

X1  yl  1 

X1  yl  Z1 

y2  z?  l 

-  y- 

X2  Z2  1 

+  z- 

xj  y2  l 

- 

x  2  y  2  z  2 

*3  Z3  1 

x3  z3  1 

x3  ^3  1 

x3  y3  z3 

(9.3.2) 

For  convenience  let  us  substitute  for  the  coefficients 
of  equation  (9.3.2)  as  follows: 


y,  z,  i 

xi  zi  1 

a  = 

y2  1 

b  =  - 

x2  z2  1 

y3  7 3  1 

x3  z3  1 

x,  y,  i 

X1  yl  Z1 

c  = 

x  2  y  2  ^ 

d  =  - 

CVJ 

fst 

CVJ 

>> 

CVJ 

X 

x3  y 3  1 

x3  y3  z3 

the  result  is  the  equation  of  a  plane  surface: 

ax  +  by  +  cz  +  d  =  0  (9.3.3) 

with  a,b,c  and  d  derived  directly  from  the  above  by  using 
the  usual  determinant  development  of  a  3x3  matrice  as 
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fol 1 ows : 

a  b  c 

d  e  f  =  a(e*k  -  f-h)  -  b(d*k  -  f*g)  +  c(e*g  -  d-h) 
g  h  k 

9J_3^2_Comgutation_of_the_azimuth 

We  are  interested  in  the  orientation  downslope;  hence 
we  have  to  distinguish  in  the  computation  of  the  azimuth 
among  several  cases. 

First  assume: 

a  ^  0  ,  b  M  ,  c  /  0  ,  d  ^  0 

This  situation  is  shown  in  fig.  9.3,  in  which  A,  B,  and 
C  are  the  three  given  points;  P,  Q,  and  R  are  the  inter¬ 
section  points  of  the  plane  with  the  axes;  the  origin 
0  =  (0,0,0);  and  the  vector  US'  signifies  the  direction 
of  the  azimuth. 

According  to  equation  9.3.3: 

P  =  (-d/a,  0,  0) 

Q  =  (0,  -d/b,  0) 

R  =  (0,  0,  -d/c) 

Now  consider  each  case  for  triangle  POO  in  relation  to 
the  sign  of  -d/a,  -d/b,  and  -d/c. 

Case  1 :  As  in  figure  9.4: 

Assume 

-d/a  >  0 
-d/b  <  0 
-c/c  >  0 

The  projection  on  plane  xy  (z  =  0)  is  shown  in  fig.  9.4. 
a  is  the  required  azimuth. 


The  downs! ope  direction  of  the  plane  is  determined  accor- 


ding  to 

the  sign 

of  -d/c 

,  because 

;  it  signifies  whether 

point  R 

on  the  z- 

axis  is 

above  or 

below  the  origin. 

In  this 

case  the 

sol uti on 

is,  cons 

iidering  the  constraints 

mentioned  above: 

tan  a 

=  -a/b 

90  <  a 

<  180 

Case  2:  Assume 

-d/a  >  0 
-d/b  >  0 
-d/c  >  0 


The  projection  on  the  xy-plane  is  shown  on  fig. 
The  solution: 


Case  3:  Assume 


9.5. 


-d/a  >  0 
-d/b  >  0 
-d/c  <  0 


It  can  be  seen  from  the  projection  on  the  xy-plane 
(fig. 9. 6)  that  this  resembles  case  2,  but  the  sign  of 
-d/c  is  reversed.  Therefore  the  direction  changes  180°. 


The  solution: 


It  is  possible  to  analyse  in  a  similar  way  all  the  cases. 
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3 

2  =  8  in  total:  four  for  d/c<0  and  four  for  d/c>0. 

Up  to  this  point  -  for  all  cases  discussed  -  a,b,c,or  d 
were  assumed  never  to  be  zero.  Let  us  now  analyse  one 
special  case,  say: 

a  =  0 ,  b  =  0 

c  t  0,  d  i  0 

The  planar  equation  reduces  to 

cz  +  d  =  0 

or  z  =  -d/c 

This  is  a  plane  which  parallels  the  xy-plane,  that  is  to 
say  that  the  elevation  of  all  three  points  is  identical. 
They  are  situated  on  a  horizontal  plane  for  which  the 
azimuth  is  not  defined,  and  the  dip  is  0. 

Although  c  can  be  zero  theoreti cally ,  it  would  mean  that 
the  plane  defined  by  our  three  points  parallels  the  z- 
axis  and  is  perpendicular  to  the  xy-plane.  In  practice 
this  can  only  happen  if  the  three  points  fall  on  a  straight 
line,  which  never  happens  in  our  procedure.  Hence  al¬ 
ways  c  /  0.  Therefore  a  Iso  a  =  b  =  d  =  0  is  impossible, 
because  then  also  c  would  be  0. 

The  solutions  for  the  various  applicable  cases  are  sum¬ 
marized  in  Table  9.1.  Parts  A  and  B  list  conditions 
whose  solutions  are  shown  in  part  C,  in  the  same  column. 

For  example,  condition  (d  =  0,  c/a<0,  c  /  b  <  0 )  is  found 
in  part  A,  column  2.  Hence  the  solution  is  tan  a  =  a/b 
and  180<a<  270.  Special  cases  are  treated  in  Table  9.2. 

9i3i3__Comgutation_gf_the_sloge_angle 


Referring  to  fig.  9.3  we  get  the  dip,  or  maximum  slope 
angle,  from  triangle  ROS  as  in  the  two-dimensional 


Table  9.1:  Solutions  for  the  orientation  in  space  of  a 

plane  defined  by  three  points  -  general  cases 

a  is  the  required  azimuth; 

-  means  that  the  value  of  the  expression  in 
the  left  column  is  less  than  zero; 

+  means  that  the  value  as  above  is  positive. 


Conditions:  a/0  b/0  c/0  d=0 


and 

c/a 

(-) 

(-) 

(  +  ) 

<  +  ) 

and 

c/b 

(  +  ) 

(-) 

(-) 

(  +  ) 

Conditions:  a/0  b/0  c/0  d/0 


and  d/a 

( +) or ( - ) 

(+)or( -) 

(  - ) or (  +  ) 

(-)or(+) 

and  d/b 

(  - ) or  (  +  ) 

(+)or(-) 

(+)or(-) 

(-)or(+) 

and  d/c 

(-)or(+) 

( -)or(+) 

( -)or (+) 

Solution 


tan  a  = 

-a/b 

a/b 

-a/b 

a/b 

Azimuth  range 

270<a<360 

180<a<270 

90<a<180 

0<o<90 

Solution 


Solution 
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figure  9.7,  wherein  the  dip  is  denoted  by  <t>. 
From  geometric  considerations 


cos  4> 


(a2  +  b2  + 


C2)0 


and  therefore  the  dip  is 


4> 


arc  cos 


,  2  .  .2  .  2.0*5 

(a  +  b  +  c  ) 


(9.3.4) 


?.:3.:4__Eya]uation_of_the_slope_angle_and_of_its_orien; 

5§5l22_f2C.5_§9y2!22 

It  has  been  shown  above  how  it  is  possib’:  to  obtain  the 
coefficients  of  an  equation  of  a  plane  defined  by  three 
points  in  space,  and  how  these  coefficients  enable  us  to 
establish  the  values  of  the  azimuth  and  the  slope  angle. 
On  the  basis  of  this  procedure,  the  problem  of  assigning 
slope  angle  and  azimuth  values  to  a  grid  square  can  be 
treated  as  follows. 

In  square  ABCD,  compute  azimuth  and  slope  angle  as  ex¬ 
plained  above  for  each  of  the  triangles  ABC,  ABD,  BCD, 

ACD  (fig.  9.8).  Now  let  us  define  the  slope  angle  of  the 
square  as  the  mean  of  the  four  slope  angles  obtained  for 
each  of  the  above  triangles.  Similary,  the  azimuth  of 
the  square  will  be  defined  as  the  mean  of  the  four  azi¬ 
muths  obtained  for  each  of  the  triangles.  Although  we 
get  in  this  way  only  an  approximation  of  a  segment  of  the 
actual  surface  by  a  plane,  it  may  actually  be  more  re¬ 
presentative  of  the  general  trend  of  the  surface,  be¬ 
cause  it  is  not  so  dependent  on  a  single  elevation  value 
as  would  have  been  the  case  with  a  triangle. 

If  all  four  points  happen,  by  chance,  to  be  exactly  on 


the  plane,  we  would  get  the  real  values  for  the  plane. 

But  in  the  general  case,  the  points  will  diverge  from 
the  plane  fitted  to  them.  The  fit  however,  tends  to  be 
good.  An  application  of  the  method  to  the  Nahal  Yael  wa¬ 
tershed  has  shown,  in  the  majority  of  cases,  a  good  re¬ 
production  of  the  slope  directions  as  expected  on  the 
basis  of  the  topographic  map.  In  certain  grid  squares 
internal  relief  causes  i naccuraci es .  For  these  the  inhe¬ 
rent  assumption  of  a  pseudo-planar  surface  is  wrong,  and 
the  only  available  solution  is  a  further  subdivision  of 
the  unit  area  (section  9.2.  above).  Such  a  subdivision 
could  be  made  to  depend  on  an  objective  criterion,  such 
as,  for  example,  a  standard  deviation  of  the  component 
triangle  angles  and  azimuths  above  a  given  threshold  va- 
1  ue . 

In  summary,  the  above  procedure  yields  a  simulation  of 
the  terrain  topography  based  on  a  square  grid  system. 

The  next  step  in  the  development  of  the  model  -  the  simu¬ 
lation  of  the  divides  and  drainage  lines  on  the  basis  of 
the  above  plane  assemblage  -  will  be  discussed  in  the 
next  section. 

9.4  Determination  of  divides  and  channels 

The  simulated  divides  and  channels,  to  be  generated 
from  the  square  grid  data,  are  designed  to  enable  the  dy¬ 
namic  functions  of  the  model.  The  grid  square  network 
will  merge  into  a  stream  channel  network,  and  through 
these  mass  transfers  (water  and  sediment)  will  take  place. 

In  our  model,  the  drainage  network  as  well  as  the  divides 
are  represented  by  sides  of  grid  squares.  A  square  side 
which  forms  a  segment  of  a  stream  channel  is  named  a 
stream  barrier,  and  a  square  side  which  forms  a  segment 
of  a  divide  is  called  a  divide  barrier.  Both  types  con- 
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stitute  barriers  in  the  sense  that  they  do  not  transmit 
mass  (water  and  sediment)  through  them.  At  most,  in  the 
case  of  stream  barrier,  mass  is  conveyed  along  the 
barri er . 

The  process  of  fluvial  mass  transfer  in  a  watershed  mo¬ 
delled  along  these  lines  is  thus  divided  into  two  phases: 
I)  the  unconcentrated  phase  (i.e.  slopeflow)  occurs  over 
the  square  surfaces;  and  II)  the  unconcentrated  phase 
(i.e.  channel  flow)  occurs  along  the  square  sides. 

The  criterion  for  determing  whether  a  given  square  side 
constitutes  a  barrier  depends  on  the  directional  trends 
(azimuths)  of  the  two  adjoining  squares.  Whenever  those 
azimuths  trend,  generally,  in  opposing  directions,  a 
barrier  is  indicated.  In  cases  where  the  azimuthal  trends 
are,  more  or  less,  parallel,  then  the  squares  they  re¬ 
present  are  part  of  the  same  slope. 

An  extreme  case  of  non-opposing  directional  trends  oc¬ 
curs  when  the  trends  of  two  neighboring  squares  are  iden¬ 
tical  (fig.  9.9).  Here  the  azimuth  of  ABCD  is,  just  like 
that  of  CDEF ,  270°.  CD  cannot  constitute  a  barrier  be¬ 
cause  the  flow  within  CDEF  is  totally  toward  CD,  and 
within  ABCD  the  flow  is  directed  from  CD  toward  AB.  A 
continuous  slopeflow  over  both  squares  is  thus  indicated. 

If  the  trends  are  generally  opposing,  two  cases  are  dis¬ 
tinguished:  I)  convergence;  and  II)  divergence. 

In  converging  trends  a  stream  barrier  is  formed,  while 
for  diverging  trends  we  get  a  divide  barrier.  Extreme 
cases  of  converging  and  diverging  trends  are  illustrated 
in  figures  9.10  and  9.11,  respectively.  In  both,  CD  is 
the  barrier. 

For  the  general  case  a  threshold  angle  can  be  used  to  de¬ 
fine  the  relationship  between  two  adjacent  opposing 
squares.  If  the  angle  between  the  azimuths  of  the  two 


i 
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squares  is  greater  than  the  threshold  -  a  barrier  exists. 
For  practical  purposes,  it  is  better  to  establish  two 
different  threshold  values  -  one  for  channel  barriers 
(called  STBAR)  and  one  for  divide  barriers  (called  DVDBAR). 

The  determination  of  the  inter-azimuthal  angle  is  shown 
in  figure  9.12.  The  trend  of  ABCD  is  designated  o.  The 
trend  of  CDEF  is  B.  The  angle  between  the  azimuths  is  de¬ 
signated  4>. 

According  to  the  various  directions  of  a,  B  we  get  by 
analytical  substitution  the  following  equation  for  the 
inter-azimuthal  angle: 

j  a  -  B  for  |  a  -  B  £  180 

<D  =  (9.4.1) 

360  -|  a  -  B  for  «  -  B  >  180 

The  conditions  for  the  existence  of  a  barrier  are: 

Channel  barrier  -  converging  azimuths  and  <*>  >  STBAR 
Divide  barrier  -  diverging  azimuths  and  $  >  DVDBAR. 

For  the  threshold  angles  there  exists  the  following  theo¬ 
retical  constraint: 

0  £  STBAR,  DVDBAR  6  180 

Clearly,  the  smaller  the  threshold  angles,  the  larger 
the  number  of  barrier  sides.  In  the  extreme  cases  of 

STBAR  =  DVDBAR  =  0 

a  channel  barrier  is  simulated  at  any  side  toward  which 
there  is  convergence,  however  slight,  and  a  divide  barrier 
will  result  at  any  side  from  which  there  is  even  a  minute 
degree  of  divergence. 
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Alternately  for 


STBAR  =  DVDBAR  =  180 
no  barriers  are  formed  at  all. 

After  finding  the  angle  between  two  azimuths,  it  is  ne¬ 
cessary  to  ascertain  whether  it  is  converging  or  diver¬ 
ging,  because  the  definition  of  <t>  (equation  9.4.1)  does 
not  enable  the  distinction,  in  itself,  between  a  thalweg 
and  a  crest  line.  As  an  example,  in  figure  9.10  <J>  =  180, 

and  in  figure  9.11  0=180  as  well;  however,  the  first 

is  diverging  and  the  second  is  converging. 

The  distinction  between  a  diverging  and  converging  barrier 
can  be  made  either  by  the  method  of  segments  or  by  com¬ 
paring  elevations. 

I )  Ih§_"§egment"_method:  Referring  to  figure  9.13,  let 
us  define 


azimuth  A5CD  =  a 

"  CDEF  =  & 

"  ADGH  -  y 

Further,  let  us  define  side  TTH,  if  indeed  it  qualifies 

as  a  barrier,  a  longitudinal  barrier.  Similary  let  us  de¬ 
fine  an  east-west  directed  side,  if  it  qualifies  as  a  bar¬ 
rier,  a  latitudinal  barrier ,  like  A U  in  fig.  9.13. 

The  divergence/convergence  test  is  done  separately  for 
longitudinal  and  latitudinal  barriers. 

!:209ity(!i[!c;!l_barr i§rs :  The  azimuths  of  two  east-west 
neighbor  squares  (such  as  ABCD  and  CDEF  in  fig.  9.13) 
converge  if  the  azimuth  of  the  eastern  square  (CDEF)  is 
between  180°  and  360°  and  the  azimuth  of  the  western 
square  (ABCD)  is  between  0°  and  180°.  If  both  conditions 
are  not  met,  there  is  no  barrier.  If  both  conditions  are 
satisfied  and  $  (the  inter-azimuthal  angle)  exceeds  STBAR, 
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a  longitudinal  channel  barrier  is  formed. 

The  azimuths  diverge  only  when  conditions  opposite  to 
convergence  prevail,  i.e. 

azimuth  ABCD  is  between  180°  and  360°; 
azimuth  CDEF  is  between  0°  and  180°! 

If,  in  addition,  the  inter-azimuthal  angle  exceeds  the 
threshold  angle  for  divides,  a  longitudinal  divide  barrier 
is  formed,  such  as  along  side  ITU  (fig.  9.13). 

^5L'rl§!2§:  For  a  latitudinal  barrier  such  as 
7HJ  in  fig.  9.13,  convergence  or  divergence  are  determina¬ 
ted  in  an  analogous  manner: 

Convergence  exists  if  90  4  a  4  270 

and  y  ^  90  or y  £  270;  and 
Divergence  exists  if  a  ^  270  or  a  £  90 

and  90  £  7^270. 

For  any  other  set  of  azimuthal  relationships  no  barriers 
between  adjacent  squares  are  formed  and  the  flow  over 
them  is,  accordingly,  slopeflow  only. 

The  application  of  the  “segment"  method  is  summarized  in 
table  4.1. 

The  solutions  in  table  9.3  refer  to  the  usual  case  of 
the  simulation  network  paralleling  the  map  coordinate  net¬ 
work.  If  another  direction  is  chosen  for  the  grid  of  the 
simulation,  the  system  must  be  rotated  by  an  angle  whose 
magnitude  is  the  difference  between  the  map  worth  and 
the  longitudinal  axis  of  the  model  grid. 

II)  Comgar i ng_el evati ons :  The  second  method  usable  for 
distinguishing  between  convergence  and  divergence  relies 
on  the  data  bank  of  topographic  elevations  as  well  as  on 
the  azimuth  values. 

Two  adjacent  squares  converge  toward  the  side  common  to 
both  if  the  mean  elevations  of  both  pairs  of  outer  cor¬ 
ners  are  higher  than  the  mean  elevation  of  the  two  ends 
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Table  9.3:  Conditions  for  the  existence  of  barriers  - 
the  "segment"  method. 

See  figure  9.13. 
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of  the  common  side.  Figure  9.14  shows  two  adjacent  squa¬ 
res  ABCD  and  CDFF,  with  TTU  as  their  common  side.  Let 
Ha,  Hg,...,  Hp  be  the  elevation  of  each  corner.  Conver¬ 
gence  toward  ITU  prevails  if 

?(HC  +  HD)  <  2(HA  +  HB> 

and 

^(Hj.  +  Hg)  <  ^(Hp  +  Hp ) 

Actually,  the  factor  ^  may  be  omitted.  Thus,  divergence 
away  from  ITU  will  prevail  if 

Hc  +  HD  >  ha  +  Hp 

and 

HC  +  HD  >  HE  +  HF 

If,  in  addition  to  the  two  convergence  criteria  detailed 
above,  we  have 


4>  >  STBAR 

then  ITU  is  a  channel  barrier. 

Likewise,  if,  in  addition  to  the  two  divergence  criteria 
listed  above,  we  have 

0  >  DVDBAR 

then  ITU  i  s  a  divide  barri  er . 

Both,  the  segment  method  and  the  elevation  method  make 
sure  by  their  structure  that  channel  barriers  and  divide 
barriers  are  mutually  exclusive:  it  is  impossible  for  a 
given  side  to  be  both  at  one  and  the  same  time.  The  seg¬ 
ment  method  has  the  advantage  of  not  requiring  another 
step  of  computation  beyond  the  determination  of  the  azi¬ 
muths  and  their  evaluation  by  segments.  In  handling  a 
large  network,  this  may  result  in  some  economizing  of 
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computer  time.  Its  main  disadvantage  is  the  fact  that 
many  borderline  cases  do  not  qualify,  in  this  method,  as 
barriers,  when  barriers  actually  do  exist  in  nature. 

This  point  will  be  elaborated  and  illustrated  in  the  next 
section. 

The  elevation  method,  though  requiring  another  computatio 
nal  step,  is  less  prone  to  border  cases  and,  furthermore, 
is  not  dependent  on  network  orientation.  The  determina¬ 
tion  if  a  given  barrier  as  a  channel  or  a  divide  is  imme¬ 
diate.  However,  inaccuracies  result  from  the  use  of  an 
elevation  average  value  as  a  criterion  for  comparison. 

The  threshold  angles  STBAR  and  DVDBAR 

Throughout  the  discussion  above  we  have  taken  for  gran¬ 
ted  the  applicability  of  the  threshold  ("critical")  an¬ 
gles  STBAR  and  DVDBAR,  which  were  supposed  to  consistent¬ 
ly  determine  whether  a  barrier  between  two  adjacent  squa¬ 
res  exists  solely  by  comparing  their  inter-azimuthal 
angle  with  these  thresholds.  This  assumption  will  be  dis¬ 
cussed  below. 

From  the  evaluation  of  maps  and  of  trial  runs,  it  has 
been  found  that  the  angles  which  yield  the  best,  optimal 
results  for  the  Nahal  Yael  watershed  are: 

For  STBAR  40°  -  45°; 
for  DVDBAR  50°  -  60°. 

The  procedure  for  the  "manual"  network  involved  the  tra¬ 
cing  of  the  drainage  lines  and  the  divides  from  a  topo¬ 
graphic  map,  and  approximating  them  by  the  nearest  grid 
link  (=  square  side).  By  combining  various  values  of 
STBAR  with  various  values  of  DVDBAR,  and  making  a  link 
by  link  comparison  of  the  resulting  computer-simulated 
network  with  the  manually  approximated  map-based  network, 
the  optimal  values  listed  above  have  been  determined. 
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Using  the  optimal  threshold  angles  does  not  guarantee 
complete  connectivity  of  the  simulated  network.  Even  for 
the  combination  of  angles  which  produces  the  largest 
number  of  barriers,  a  certain  number  of  links  are  not 
reproduced,  while  elsewhere  unjustified  links  appear  in 
the  simulated  network.  An  extreme  example  of  the  latter 
occurs  for  the  theoretical  case  of 

STBAR  =  DVDBAR  =  0 

As  a  means  of  "filling  the  holes"  in  the  network,  it  is 
necessary  in  future  to  develop  a  technique  which  will 
connect  disjointed  links  of  the  same  type,  i.e.  stream 
barriers  with  stream  barriers  and  divide  barriers  with 
divide  barriers.  Such  a  procedure  might  be  based  on  the 
principle  of  searching  out  the  nearest  barrier  of  the  sa¬ 
me  type  and  connecting  to  it  without  crossing  another 
barrier . 

Whenever  the  network  as  simulated  by  the  optimal  angles 
method  has  only  few  network  discontinuities,  it  may  be 
more  practical,  for  a  local  project,  to  achieve  complete 
connectivity  manually. 

Evaluation  of  the  results  obtained  for  various  combina¬ 
tions  of  threshold  angles  revealed  that  very  often  a 
barrier  is  not  generated  at  links  which,  in  fact,  are  bar 
riers,  when  the  adjoining  square(s)  have  a  considerable 
internal  relief.  This  means  that,  for  those  squares,  the 
corner  elevations  do  not  adequately  represent  the  terrain 
therefore,  the  azimuth  -  a  sinle  parameter  based  solely 
on  these  elevations  -  is  likewise  not  adequately  repre¬ 
sentative,  and  a  further  subdivision  is  thus  called  for. 

In  summary,  although  fully  connected  networks  of  stream¬ 
lines  and  divides  are  not  generated  even  by  the  optimal 
STBAR  -  DVDBAR  angle  combination,  a  large  part  of  the 
network  is  indeed  produced.  As  a  future  development,  the 
threshold  angles  could  be  made  to  vary  over  area  as  a 
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function  of  the  slope  angle  relationships  between  adja¬ 
cent  squares,  or  even  in  conjunction  with  some  non-geo- 
metric  terrain  characteristic. 

9.5,  Substitution  of  solutions  in  the  computer 
9^5^!. _ iQ^roductign 

Following  the  theoretical  solutions  arrived  at  for  a 
single  square  (section  9.3)  and  for  barriers  between 
squares  (section  9.4),  a  computer  program  in  Fortran  IV 
has  been  written  for  the  application  of  these  solutions 
to  a  field  of  squares.  The  main  input  of  the  program  are 
the  elevation  data  of  the  square  corners  as  read  from 
the  topographic  map.  The  main  output  of  the  program  is: 

I)  The  mean  elevation  for  each  square; 

II)  The  azimuth  of  each  square; 

III)  The  slope  angle  of  each  square; 

IV)  The  barrier  network  resulting  from  a  given  combina¬ 
tion  of  inter-azimuthal  threshold  angles. 

This  output  was  also  converted  to  a  graphic  display  by 
means  of  a  Gerber  plotter. 

In  this  section  the  computer  programs  are  explained  and 
evaluated,  the  methods  of  computation  are  detailed,  and 
the  results  are  discussed. 

_area  _?!3?L§ome_technical  _remarks 

The  area  whose  map  served  as  tne  base  for  trying  the 
computer  programs  is  located  in  the  upper  parts  of  the 
Naha!  Yael  Research  Watershed,  and  comprises  the  whole 
of  watershed  05  and  a  part  of  watershed  04.  A  square 
grid  drawn  on  a  map,  scale  1 : 1 ,250, yielded  about  1,000 
squares  12-5  meters  square.  Each  unit  element  had  an  or- 
thogonally  projected  surface  area  of  156-25  m  . 


The  corner  elevation  of  all  squares  was  read  manually 
off  the  map  and  punched  on  a  data  card.  (Automatic  rea¬ 
ding  devices  might  perform  the  same  function  either  off  a  map  or 
off  an  air  photograph  whenever  the  model  is  applied  to 
large  areas).  The  axes  of  the  grid  used  run  parallel  to 
the  map  coordinate  system.  Point  (250,800)  is  the  south-  , 
western  corner  of  the  source  map  and  was  designated  (1,1) 
in  the  model.  From  this  followed  the  model  coordinates 
in  such  a  way  that  the  index  of  the  easternmost  corner 
was  41  and  that  of  the  northernmost  corner  was  26.  Based 
on  this  grid,  an  individual  square  is  identified  by  the 
coordinates  of  its  southwestern  corner.  Numerical  results 
for  each  square  appear  in  the  output  in  tabular  form, 
with  rows  and  columns  parallel  to  the  coordinate  network. 

Index  I  denotes  the  latitudinal  ( east-wes+ 1  axis,  with 
the  highest  value  designated  II. 

Index  J  denotes  the  longitudinal  ( north-south)  axis, with 
the  highest  value  designated  OJ. 

The  grid  point  network  must  clearly  exceed  the  boundaries 
of  the  studied  area,  so  that  the  entire  divide  zone  is 
covered  by  full  squares,  as  in  figure  9.15.  All  corners 
of  all  shaded  squares  must  be  included  in  the  data  input, 
even  though  most  of  them  are  located  outside  the  divide. 

In  the  area  used  here  as  an  example,  considerable  areas, 
mainly  in  the  SW  part  of  the  grid  actually  lie  beyond 
the  divide  of  the  Nahal  Yael  watershed. 

In  the  following  pages  the  computer  program  and  its  tabu¬ 
lar  and  graphic  outputs  are  described  in  more  detail. 

The  program  is  available  upon  request. 

•i?.::!i3__Comgt2ter_printout_Noi_l_{SJFCOAY) 

This  program  includes  the  main  inputs  and  outputs.  It  is 
based  on  the  scanning  by  loops  of  the  corners  which  be¬ 
long  to  given  squares  in  the  main  program.  For  every 
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square  reached  by  the  scanning  process,  the  required  out¬ 
put  is  arrived  at  by  a  suitable  sub-program  of  a  function 
or  subroutine  type.  All  the  computations  relating  to  an 
individual  square  (e.g.  the  slope  angle)  are  performed 
within  one  major  loop  (from  statement  13  to  statement  65, 
inclusive).  Another  loop  performs  the  interrelationships 
between  two  adjacent  squares  -  the  barrier  decisions 
(from  statement  66  to  statement  72).  After  these  two 
main  parts  there  follow  the  printout  statements:  first 
for  the  graphic  output  (statement  73:  CALL  6ERBRS)  and 
then  the  numerical  output  which  appears  in  the  printout. 

PROGRAM  STAGES 

I^  Readlng  the  data ;  The  data  -  elevations  of  the  square 
corners  -  are  read  into  an  E(II,JJ)  matrice  (in  our  exam¬ 
ple  E(41,26))  and  these  data  appear  as  the  first  table 
of  the  output.  The  reading  is  done  line  by  line.  The 
length  of  each  line  is  II  and  the  number  of  the  lines  is 
JJ.  This  rectangular  field  may  include  areas  for  which 
there  are  no  data,  as  they  are  outside  the  watershed  and 
its  divide  zones.  For  these  points  some  impossible  eleva¬ 
tion  value  is  punched-,  in  our  case,  where  all  real  ele¬ 
vations  exceed  200m,  we  have  punched  0  (zero).  For  the 
indices  which  refer  to  these  points  no  computations  were 
performed  and  all  the  variables  computed  for  them  were 
put  as  "blank".  The  computer  cell  which  represents  any 
one  of  these  points  does  not  contain  any  value.  The  check 
whether 

E(I,J)  =  0 

is  performed  before  the  beginning  of  the  computations 
(statement  15). 

In  the  output,  every  location  for  which  such  a  blank  is 
indicated,  there  appears  a  character  different  from  all 
others:  N  or  -0. 
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_II_)_Com2ut.ing_ithe _n}e a n _elevat-ion_for_each  square:  As 

each  square  is  identified  by  the  coordinates  of  its  lower 
left  corner,  the  mean  elevation  of  square  (I,J)  is  (see 
figure  9.16): 

0 • 25 ( E ( I ,  J  )  +  E ( 1  +  1 ,  J  )  +  E ( I +1 , J  +  l )  +  E ( I ,  J  +  l ) ) 

The  computation  is  performed  by  FUNCTION  AVG  which  cal¬ 
culates  the  arithmetic  mean  of  N  arrayed  values.  The  ma- 
trice  of  the  mean  elevations  is  designated  AH.  Its  size 
is  (II-l.JJ-l)  or,  in  our  example,  AH(40,25).  It  appears 
as  the  second  table  of  the  output. 

II  Il_Calculat2ng_azimut.h  and_s2ope  angle  f  or  _each_sguare : 

First  a  test  is  made  for  a  special  case,  in  which  all 
corners  of  the  given  square  are  of  equal  plevation  (sta¬ 
tement  21).  In  such  a  case  the  azimuth  is  undefined,  be- 
cause  we  have  a  horizontal  plane.  The  azimuth  is  then 
designated  blank  and  the  slope  angle  is  zero.  In  our  set 
of  1,000  squares  there  was  no  such  case.  Evidently  it  is 
unlikely  in  fluvially  conditioned  topography  in  general. 

In  all  other  cases  it  is  necessary  to  calculate  azimuth 
and  slope  angle  separately  for  each  combination  of  three 
out  of  the  four  corners  (section  9.3  above).  Each  com¬ 
bination  was  designated  M,(M  =  1,2, 3, 4). 

Assume  we  have  a  given  M  combination  of  3  points.  The 
values  for  the  coefficients  of  the  plane  equation  AT(M), 
BT(M),  CT(M),  DT(M)  are  obtainable  through  FUNCTION  DET, 
which  calculates  the  determinant  of  a  mat)  ice  of  the  or¬ 
der  3x3  (see  paragraph  9.3.1  above).  The  resulting  co¬ 
efficient  values  have  a  scale  distortion,  because  the 
horizontal  distance  between  adjacent  corners  is  one  unit 
(fig.  9.16),  and  the  vertical  difference  between  corner 
points  is  in  meters;  however,  this  distortion  does  not 
affect  our  procedure,  because  we  are  interested  only  in 
the  relative  values  of  the  coefficients  (except  for  one 
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case  to  be  discussed  below). 

The  four  coefficients  AT(M),  BT(M),  CT(M),  DT(M)  consti¬ 
tute  the  input  to  SUBROUTINE  AZIDIP,  which  uses  them  to 
calculate  azimuth  and  slope  angle  (see  paragraph  9.3.2 
above).  In  the  process,  the  variables  TAZT(M)  -  the  "tem¬ 
porary"  azimuth  for  the  particular  M  combination  and 
TDIP(M)  -  the  "temporary"  dip  for  same,  are  calculated. 

By  repeating  this  four  times,  i.e.  for  each  M  combina¬ 
tion,  we  get: 

TAZT(l), . TAZT  ( 4 ) 

and 

TDIP(l), . TDI P  (4) 

The  final  mean  values  are  obtained  through  AVG,  and  the 
resulting  azimuth  and  slope  angle  for  square  ( I , J)  are 
designated,  respectively 

AZIMUTH  (I,J),  DI P ( I  ,  J  )  . 

A  special  case  occurs  whenever  a  certain  square  has  three 
corners  of  equal  elevation.  Then  only  the  slope  angle  is 
calculated  by  AZIDIP,  while  the  azimuth  depends  entirely 
on  the  fourth  point,  which  will  be  either  higher  or  lower 
than  the  three  others  (statements  42  to  56  inclusive). 
Such  a  case  indeed  occurred  in  our  example,  square  (24, 
14). 

By  comparing  the  results  as  computed  by  the  model  with 
the  source  map,  a  generally  dependable  reproduction  of 
the  terrain  is  indicated,  both  for  orientation  and  for 
slope  angle.  Clearly,  the  more  terrain  homogeneity  with¬ 
in  a  square  (extreme  homogeneity:  parallel  and  equidis¬ 
tant  contours),  the  higher  the  representativeness  of  the 
computed  azimuth  and  dip.  /> 

In  practice,  the  azimuth  and  dip  values  represent,  for 
most  squares,  the  characteristics  of  the  topography.  In¬ 
ternal  relief  within  individual  squares  -  such  as  caused 
by  very  densely  dissected  terrain,  junction  areas  of 
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stream  channels,  or  very  localized  cliffs  -  detracts  from 
the  representati vi ty ,  since  such  varied  squares  require 
more  than  one  set  of  two  numbers  to  adequately  portray 
their  topographic  composition. 

The  computer-generated  topography  is  generally  more  mode¬ 
rate  than  the  real-world  prototype,  i.e.  its  slope  an¬ 
gles  are  somewhat  lower.  This  is  so  because  the  randomly 
placed  grid  intersections  often  miss  the  very  highest, 
as  well  as  the  very  lowest,  points  in  the  terrain.  Usu¬ 
ally  these  special  points  cover  only  a  small  fraction  of 
the  terrain,  but  their  role  in  landscape  conditioning 
and  evaluating  is  cons i derabl e .  The  upper  part  of  the 
Mahal  Yael  watershed  has  a  dissected  and  rapidly  changing 
topography.  It  is  possible  that  even  better  representa- 
tivity  would  have  been  obtained  for  an  example  area 
with  a  more  rounded  and  orderly  fluvial  topography. 

In  summary,  we  conclude  that  the  simulated,  computer-ge¬ 
nerated  terrain  is  sufficiently  representati ve  of  the  ac¬ 
tual  topography  and  is  this  suitable  for  the  modelling 
of  unconcentra ted  slope  runoff.  The  terrain  simulation 
would  be  further  improved  by  the  introduction  of  a  square 
size  which  will  vary  as  a  function  of  the  internal  relief 
(see  section  9.2  above). 

SUBPROGRAM  AZIDIP:  The  input  to  AZIDIP  are  the  coeffi¬ 
cients  of  the  planar  equation  A,B,C,D  and  the  output  are 
the  azimuth  and  the  dip  (=  slope  angle).  The  program  is 
executed  after  Tables  9.1  and  9.2  and  after  paragraph 
9.3.3  above.  The  sub-program  consists  essentially  of  com¬ 
paring  coefficients  and  evaluating  their  signs.  The  re¬ 
sult  is  multiplied,  for  the  azimuth  as  well  as  for  the 
dip,  by  180/r. ,  since  the  computer  calculates  angles  in 
radians. 

In  calculating  the  dip,  statements  28  and  29  read,  res¬ 
pectively 
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A  =  A/12-5 
B  =  B/12-5 

These  statements  correct  the  scale  distortion  brought 
in  during  the  calculation  of  A,B  in  the  determinant,  by 
multiplying  first  by  1/1250  (the  map  scale)  and  further 
by  100  (the  number  of  centimeters  in  one  meter).  Remem¬ 
ber  the  map  distance  between  adjacent  corners  in  1  cm. 

For  C  this  correction  is  unnecessary  because  it  is  cal¬ 
culated  solely  in  units  of  cm,  while  D  is  not  involved 
at  all  in  the  calculation  of  the  dip.  In  evaluating  the 
azimuth  there  is  no  scale  distortion;  it  is  calculated 
from  coefficient  relationships  so  that  scale  is  irrele- 
v  a  n  t . 

IV^_Calculationof_barriers ;  This  program  utilizes  the 
method  of  segments  (section  9.4).  Its  input  are  the  azi¬ 
muths  of  the  squares  (from  paragraph  9.3.3).  The  main 
program  scans  all  sides  common  to  any  two  squares  in  a 
network.  The  decisions  whether  a  given  side  is  a  barrier 
or  not,  and  if  it  is  -  of  which  type,  is  done  through 
the  BAR  function. 

The  scanning  procedure  verifies  the  barrier  property, 
for  square  (I,J),  by  testing  the  i nter-azimuthal  angle 
in  relation  to  the  threshold  angle(s)  (section  9.4)  imme¬ 
diately  above  (north)  and  immediately  to  the  right  (east) 
of  point  ( I ,  J ) .  Referring  to  figure  9.17,  while  scanning 
(I,J),  i.e.  point  A,  the  BAR  function  determines,  first, 
whether  side  4B  is  a  latitudinal  barrier  and,  if  so,  of 
what  type.  Then  it  determines  whether  side  AD  is  a  lon¬ 
gitudinal  barrier  and,  if  so,  of  what  type.  The  scanning 
process  covers  all  points  of  the  network.  Thus,  every 
side'  of  the  network  is  covered  and  the  pattern  of  barriers 
is  oh  Id  Hitd  , 

In  the  numerical  output,  the  numbers  which  appear  in  the 
barrier  table  indicate  the  various  possible  combinations 
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of  barriers  for  the  pair  of  sides  tested  for  each  grid 
point.  In  the  main  program  the  barriers  were  designated 
by  BARIER(I,J)  according  to  Table  9.4. 

From  the  way  the  barriers  are  listed,  it  follows  that 

BARI ER  (1,1)  =  0 

This  value  is  assigned  a  priori,  because  sq  (1,1)  has  no 
barrier  along  the  two  sides  which  meet  at  point  (1,1), 
which  is  the  southwestern  corner  of  the  whole  network. 

In  evaluating  barriers  for  the  first  row  (J  =  1),  only 
longitudinal  barriers  may  result.  Likewise,  for  the  first 
column  (I  =  1),  only  latitudinal  barriers  are  permitted. 

The  input  for  function  BAR  is  as  follows: 

AZIMUTH  (I  ,J)  labelled  AZIU 
AZIMUTH  (1-1,0)  "  AZIU 

AZIMUTH  (1,0-1)  "  AZIJ1 

There  is  also  an  input  parameter  labelled  IN.  Whenever 
IN  =  11,  then  J  =  1 ,  we  are  in  the  first  row,  and  conse¬ 
quently  latitudinal  barriers  are  out.  Likewise,  whenever 
IN  =  1  we  are  in  the  first  column,  and  consequently  lon¬ 
gitudinal  barriers  cannot  result.  In  normal  cases  IN  =  0. 

Further  input  are  the  values  of  the  threshold  angles 
STBAR,  DVDBAR  which  result  from  the  COMMON  statement.  The 
output  of  the  BAR  function  is  also  named  BAR,  and  the 
stages  of  calculation  follow  Table  9.4,  as  described 
bel ow : 

First,  equation  9.4.1  is  used  to  calculate  the  inter-azi¬ 
muthal  angle.  DELV  is  the  angle  between  AZIJ  and  AZIIJ, 
and  OELH  is  the  angle  between  AZIJ  and  AZI01.  We  set 
BAR  =  0  and  test  all  the  possible  cases  one  by  one.  As  a 
start,  we  ask:  is  there  a  longitudinal  channel  barrier? 
This  question  is  analogous  to: 

I)  are  the  relevant  azimuths  converging?  and 

II)  is  DELV  >  STBAR? 
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Table  9.4:  The  listing  of  barriers  in  the  output 


BARIER(I,J) 

Meaning 

Graphic 

0 

No  barrier  between  sq(I,J)  and  either 
its  western  or  its  southern  neighbour 

n 

1 

Longitudinal  channel  barrier  between 
sq(I,J)  and  its  western  neighbour 

I 

2 

Latitudinal  channel  barrier  between 

sq(I,J)  and  its  southern  neighbour 

■ 

1 

3 

Channel  barrier  of  type  1  and  2  to- 
gether 

1 

-1 

Longitudinal  divide  barrier  between 
sq(I,0)  and  its  western  neighbour 

1 

■ 

1 

-2 

Latitudinal  divide  barrier  between 

sq(I,J)  and  its  southern  neighbour 

I 

-3 

Divide  barrier  of  type  -1  and  -2  to- 
gether 

1 

■ 

1 

4 

Longitudinal  channel  barrier,  type  1 
and  latitudinal  divide  barrier, type  -2 

1 

1 

1HBBM 

Latitudinal  channel  barrier,  type  2 
and  longitudinal  divide  barrier, type  -1 
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If  both  conditions  exist,  then 

BAR  =  1 

Now  we  move  on  to  the  question:  is  there  a  latitudinal 
channel  barrier?  Again,  first  we  check  if  there  is  con¬ 
vergence,  and  then  if 

DELH  >  STBAR 

If  the  answer  to  both  is  positive,  we  deduce 

BAR  =  BAR  2 


Now,  if  previously 

BAR  =  0 

we  get 

BAR  =  2 

But  if  previously 

BAR  =  1 

we  get 

BAR  =  3 

The  procedure  of  testing  for  divide  barriers  is  exactly 
analogous.  After  going  through  the  entire  procedure  of 
testing  for  all  possibilities,  we  finally  determine  BAR 
to  be  one  of  the  nine  values  listed  in  Table  9.4. 

For  computer  printout  No.l  the  following  values  were 
substituted  for  the  threshold  angles: 

STBAR  =  30 
DVDBAR  =  60 


The  best  simulation  of  the  drainage  network  is  obtained 
somewhere  in  the  vicinity  of  these  values  (see  section 
9.4),  and  the  results  are  presented  numerically  in  the 
last  matrice  of  the  output.  This  table,  which  describes 
the  barriers  numerically,  will  serve  in  future  for  the 
transfer  of  flow  along  the  channel  barriers.  It  also  pro¬ 
vides  the  basis  for  the  graphic  display. 
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V_i_Graghic_outgut_Noil^  This  is  the  last  stage  of  the 
present  program.  It  translates  the  numerical  results  in¬ 
to  a  graphic  display,  scale  1:1,250,  code  name  SJFCOAY. 

The  instrument  used  was  a  Gerber  plotter,  and  all  the 
read-in  operations  for  the  Gerber  were'  carried  out  by  a 
special  sub-program  named  GRBER  8. 

The  input  for  GRBER  8  is  the  elevation  of  the  grid  points, 
the  azimuth  of  the  squares,  the  barriers,  and  the  threr 
shold  angles  STEAR  and  DVDBAR  (in  the  present  case, 

ST  =  30  and  DVD  =  60).  The  graphic  output  is  produced 
step  by  step,  in  the  same  scale  as  the  source  map,  and 
consists  of:  the  axes  of  the  rectangular  grid,  its  boun¬ 
daries,  and  a  designation  of  its  cardinal  directions; 
plots  of  arrows  of  a  standard  length  from  the  center  of 
each  square  in  the  direction  of  the  azimuth;  plotting  the 
divide  lines  according  to  the  divide  barriers  with  a 
full  black  line;  plots  of  stream  lines  according  to  the 
channel  barriers  with  a  full  blue  line,  with  an  added 
blue  arrow  pointing  towards  the  lower  end  of  each  barrier 
line. 

As  stated  earlier,  the  result  is  not  perfect.  In  addition 
to  occasional  lack  of  connectivity  even  for  the  "best" 
threshold  angles,  the  channel  segments  composed  of  a  num¬ 
ber  of  grid  links  are  not  always  unidirectional.  Both  in¬ 
adequacies  are  caused,  ultimately,  by  the  fact  that  we 
use  four  points  to  represent  the  entire  square.  As  long 
as  we  insist  doing  that,  we  have  to,  subsequently ,  go 
through  a  process  of  connection  and  redirection  of  links, 
be  it  by  a subjecti ve  or  an  objective  process. 

Although  the  barrier  calculation  does  not  differentiate 
between  large  and  small  channels,  the  square  size  governs 
the  sensitivity  to  the  smaller  channels.  In  our  simulation 
some  very  small  channels  were  obliterated  because  the 
contour  indentation  effect  of  their  valleys  was  totally 
absorbed  within  a  single  square. 
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The  graphic  output  shows  some  divide  barriers  touch  a 
channel  barrier,  i.e.  they  have  a  common  end  point.  These 
cases,  represented  by  BAR  =  4  and  BAR  =  -4  (Table  9.4), 
are  inconsistent  with  normal  fluvially  conditioned  topo¬ 
graphy  and  a  procedure  must  be  devised  to  eleminate  them. 

Figure  9.18  is  a  corrected  and  "smoothed"  reproduction  of 
this  output. 

9^5^4 _ Comguter_gri!2t;out_No_;_2_(_SJFCOEK^ 

In  its  main  aspects  this  program  is  similar  to  program 
No.l.  The  difference  is  that  here  no  graphic  output  was 
produced,  and  the  barriers  were  identified  by  comparing 
elevations.  Using  this  method,  the  barriers  were  evalu¬ 
ated  for  any  combination  of  angle  values  for  STBAR  and 
DVDBAR ,  with  each  being  assigned  successive  values  at 
20°  intervals  (i.e.  0,  20,  40,  60,  ...,  180). 

The  procedure  employs  the  same  notation  and  order  of 
calculation  for  the  two  sides  which  meet  at  the  grid 
point  of  square  (I,J)  as  in  program  No.l.  The  barriers 
are  evaluated  via  the  BRR  function,  supplemented  by  the 
NC  function.  The  barrier's  notation  is  BRR. 

The  input  ingredients  for  BRR  are:  the  elevations  of  the 
intersections  of  the  rectangular  grid;  the  azimuths  re¬ 
quired  for  the  barrier  test  AZIMUTH  (I,J),  AZIMUTH  (1-1,0), 
AZIMUTH  (I,J-1);  the  parameter  IN,  which  functions  in  the 
same  way  as  in  function  BAR  (computer  printout  1);  the 
index  (I,J)  of  the  given  square;  and  the  size  of  the  ele¬ 
vation  matrice  E,  i.e.  E  (II, JO). 

First  the  inter-azimuthal  angles  DELV  and  DELH  are  de¬ 
fined,  as  in  BAR.  Now,  after  figure  9.17,  the  test  whe¬ 
ther  AD  is  a  barrier  is  performed  in  two  steps.  First, 
function  NC  is  called  to  evaluate  the  elevations  of  the 
end  points  of  EF  and  CB  versus  AD  and  to  determine  whe¬ 
ther  there  is  convergence  (NC  =  10),  divergence  (NC  =  -10), 
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or  unconcentrated  si opef 1 ow. (NC  -  0). 


If  NC  =  0,  the  program  passes  on  to  test  a  possible  la¬ 
titudinal  barrier  at  AB. 


If  NC  =  -10,  i.e.  there  is  divergence,  the  existence  of 
a  possible  divide  is  tested,  i.e.  whether 

DELV  >  DVDBAR 

and  BRR  gets  its  value  accordingly. 

If  NC  =  10,  i.e.  there  is  convergence,  and  then,  if  also 

DELV  >  STBAR 


then 


BRR  =  1 


and  then  the  program  moves  on  to  test  AB  in  a  similar 
way  for  the  possible  existence  of  a  latitudinal  barrier. 


Through  NC  the  trend  of  either  convergence,  divergence, 
or  slopeflow  is  determined,  and  the  final  value  of  BRR 
is  assigned  according  to  the  relationships  between  each 
of  the  two  threshold  inter-azimuthal  angles  and  DELH. 


In  the  numerous  barrier  matrices  we  obtained,  it  is  stri¬ 
king  to  observe  the  diminution  in  the  number  of  barriers 
with  the  increase  of  the  threshold  angle  values,  until,  at 

STBAR  =  DVDBAR  =  180 

none  remain  -  a  necessary  corollary  of  the  principle 
used  for  the  establishment  of  barriers.  It  is  also  pos¬ 
sible  to  observe  the  tendency  of  STBAR  to  be  smaller 
than  DVDBAR;  for  higher  angles  an  increase  in  the  thre¬ 
shold  values  is  especially  felt  in  a  drastic  decrease  in 
the  number  of  channel  barriers. 

A  rigorous  objective  criterion  for  the  determination  of 
optimal  values  for  the  inter-azimuthal  threshold  angles 
has  not  been  developed  at  this  stage. 


\ 


I 


MANUAL  SIMUt  AT  ION  OF  STREAM  BARRIERS 
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9 .5^5 __Ihe_graghlc  outgo t_f2r_var iou spangles 

Graphic  output  No. 1:  This  output  has  been  obtained  from 
the  program  of  printout  No.l.  It  determines  barriers  via 
the  method  of  sectors  for  STBAR  =  30,  DVDBAR  =  60,  and 
displays  them  at  the  scale  of  1:1,250  (this  same  scale 
is  used  for  all  graphic  outputs). 

Graphic  output  No. 2  (SJFC09T):  This  is  a  manual  simula¬ 
tion  performed  on  a  part  of  the  example  area  by  tracing 
the  rectangular  grid  onto  the  source  map  and  assigning 
barriers  to  square  sides  subjectively  judged  to  be  most 
closely  representative  of  the  generally  irregular  real- 
world  "barrier"  lines  {fig.  9.  4  ).  Because  of  time  con¬ 
straints  only  the  stream  channel  lines  were  simulated  in 
this  way.  The  channel  barriers  so  determined  were  then 
transferred  in  the  usual  notation  to  the  NREAL  matrice. 

Although  the  resulting  barrier  network  faithfully  repro¬ 
duces  the  salient  features  of  the  real-world  drainage 
network,  here  too  some  directional  distortions  do  occur. 

Graphic  output  No. 3  (SJFC094):  This  is  the  barrier  map 
for  STBAR  =  0,  DVDBAR  =  0,  i.e.  a  barrier  is  indicated 
wherever  the  trends  are  opposing,  however  slightly.  The 
result  shows  many  closed  depressions  as  well  as  numerous 
superfluous  channel  barriers  and  divide  barriers.  The 
barriers  are  of  type  BAR. 

Graphic  output  No. 4  (SJFCODY):  This  also  shows  the  bar¬ 
riers  for  STBAR  =  0,  DVDBAR  =  0,  but  the  barriers  are  of 
type  BRR.  The  output  is  similar  to  No. 3.  Also  here  we 
have  many  depressions  which  do  not  exist  in  nature  as 
well  as  many  superfluous  barriers. 

Graphic  output  No. 5  (S0FC0AW):  Shows  barriers  of  type 
BAR  for  STBAR  =  DVDBAR  =  150.  It  clearly  emerges  that  this 
angle  is  too  large:  barriers  are  produced  only  when  the 
inter-azimuthal  angle  for  two  adjoining  squares  exceeds 
150°,  i.e.  nearly  completely  opposing  trends  (fig.  9*19). 
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SJFCOflW  25/10/79 


Fig.  9.19:  Computer- si mu la ted  drainage  network,  upper  Nahal  Yael 
watershed.  STBAR  *  150°;  DVDBAR  *  150°.  Gerber  output 
SJFCOAW. 
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Function  BRR  yields,  for  the  threshold  angles,  even  less 
barriers. 

Graphic  output  No. 6  (SJFC0D5):  Barrier  computation  based 
on  STBAR  =  DVDBAR  =  45,  type  BRR. 

Graphic  output  'No. 7  (SJFC0D8):  Barriers  based  on  STBAR  =  30 
DVDBAR  =  60  (fig.  9.20). 

In  summary,  low  values  of  threshold  angles  yield  -  for 
either  barrier  type  (BAR  and  BRR)  -  numerous  superfluous 
barriers.  This  is  best  evident  from  a  comparison  of 
graphic  outputs  No. 3  and  4  with  the  manual  simulation 
and  with  the  source  map.  As  we  increase  the  values  of  the 
two  threshold  angles,  the  number  of  barriers,  and  espe¬ 
cially  of  superfluous  barriers,  decreases.  Very  high 
threshold  values,  such  as  used  for  graphic  output  No. 5, 
yield  only  a  relatively  small  number  of  barriers,  but  all 
the  barriers  indicated  do  agree  with  real-world  barriers. 
Conversely,  however,  in  cases  that  produce  many  super¬ 
fluous  barriers,  there  still  remain  certain  grid  links 
which  are  barriers  in  the  real-world,  as  represented  by 
the  numeral  simulation,  but  are  not  indicated  as  such  by 
the  program.  Therefore,  the  most  significant  set  of  thre¬ 
shold  angles  can  be  defined  as  the  one  generating  the 
largest  number  of  barriers  which  include  the  smallest 
number  of  superfluous  barriers. 

The  probable  effect  of  the  sharp  and  dissected  topogra¬ 
phy  of  the  prototype  area  in  relation  to  the  square  size 
has  already  been  mentioned.  To  this  might  be  added  the 
fact  that  much  of  the  prototype  area  is  in  reality  a  di¬ 
vide  zone  (see  graphic  output  No. 2).  So  further  tests 
for  other  topographic  regions  will  be  needed  if  the  model 
is  adopted  for  more  universal  use. 

Comparing  barriers  resulting  from  type  BAR  definitions 
with  those  resulting  from  type  BRR  definitions,  it  emer¬ 
ges  that  the  comparative  elevations  method  requires 
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"optimal"  threshold  angles  which  are  lower  in  value  than 
required  for  the  sector  method.  Compare  graphic  outputs 
No.l  and  No. 7,  where,  for  the  same  limiting  angles,  we 
get  more  barriers  from  the  sector  method  than  from  the 
method  using  comparative  elevations. 

Whenever  real-world  barrier  lines  follow  a  north-south 
or  an  east-west  direction,  even  approximately,  the  re- 
presentati vi ty  by  the  BAR  procedure  (=  the  sector  method) 
is  good.  Conversely,  for  reaches  of  streams  (and  divide 
lines)  oriented  roughly  at  45°  to  one  of  the  cardinal 
directions,  barrier  production  is  only  partial.  This  ef¬ 
fect  is  considerably  increased  whenever  internal  relief 
causes  at  least  one  of  the  two  azimuths  involved  in  a 
barrier  determination  to  be  unrepresentative. 

In  summary,  similar  networks  are  generated  by  the  model 
for  both  methods  for  barrier  determination  (by  sectors 
and  comparing  elevations),  though  the  optimal  threshold 
values  are  different  for  each  method.  In  neither  case 
complete  connectivity  of  drainage  lines  or  divide  lines 
could  be  achieved. 

_ ^BBllS?l5i2!]_2f!_^he_grggrani_to_other_areas 

One  objective  of  our  work  on  this  simulation  model  was 
to  enable  its  use  in  other  areas  with  minimal  changes, 
and  the  computer  programs  have  been  written  with  this  in 
mind.  We  assume  that,  for  other  regions,  the  same  in¬ 
put  -  elevations  of  grid  points  (rectangular)  as  well  as 
the  cardinal  directions  (N-S,  E-W)  will  be  used,  so  that 
no  rotation  will  be  necessary. 

The  following  statements  have  to  be  changed  when  apply¬ 
ing  the  program  to  another  region  (see  printout  No.l): 

I)  Statement  3: 

DATA  II, JJ/ . / 

which  identifies  the  northernmost  index  of  the  grid  JJ 
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and  the  easternmost  index  II  (remember  the  southwestern 
grid  point  is  indexed  .( 1 ,1 )) . 

II)  Statements  1  and  2,  which  signify  the  dimensions  of 
the  variables  for  II,  JJ  should  be,  for  the  grid  point 
matri ce 

E  ( 1 1 ,  J  J  ) 

and  for  the  other  matrices 

AH ( 1 1  —  1 , J J-l ) ,  DIP( 1 1 -1  ,JJ-1) , 
AZIMUTH(II-1,JJ-1) ,  B AR I E R ( 1 1  —  1 ,JJ-1) 

III)  Statements  7  and  8  should  be  changed  according  to 
the  selected  values  for  the  threshold  angles. 

IV)  In  sub-program  AZIDIP,  statements  28  and  29  (A  = 
A/12-5;  B  =  B/12-5)  should  be  changed  according  to  the 
scale  of  the  map  and  the  units  of  grid  point  elevations 
used . 

V)  In  subprogram  GRBER  8  the  statement  FACT0R(  0,394) 
(number  9)  is  also  contingent  upon  the  scale  of  the  map 
used . 

9.6  Summary  and  conclusions 

The  work  reported  here  is  only  the  first  stage  in  the 
development  of  a  comprehensive  deterministic  and  distri¬ 
buted  geo-system  oriented  model.  Most  of  the  work  dealt 
with  the  "infrastructure"  of  such  a  model  -  the  simula¬ 
tion  of  the  topography  and  the  construction  of  a  rectan¬ 
gular  reference  system.  Via  this  system  the  information 
input  will  be  fed  into  the  model,  and  through  it  also 
the  output  will  be  received. 

The  terrain  simulation  was  accomplished  by  calculating 
azimuth  and  slope  angle  for  each  grid  square.  The  drai¬ 
nage  network  and  divide  lines  were  simulated  on  the  ba¬ 
sis  of  assuming  the  existence  of  threshold  inter-azimu¬ 
thal  angles  for  pairs  of  adjacent  squares,  which  define 
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convergence  to  a  channel  or  divergence  from  a  divide. 
Optimal  threshold  angles  (30  -  45°  for  channels,  45  - 
60°  for  divides)  yield  in  general  acceptable,  though  in¬ 
complete  networks,  but  some  inconsistencies  still  exist. 

Two  tasks  left  to  be  done  concern  the  development  of  a 
criterion  for  the  optimization  of  the  threshold  angles 
and  the  development  of  a  "connectivity  function"  for  the 
resulting  drainage  network  and  divide  lines. 

Subjective  manual  filling  in  of  lines  and  occasional  re¬ 
versing  of  computer-indicated  flow  direction  could  sub¬ 
stitute,  for  a  specific  area,  for  the  "connectivity  func¬ 
tion",  but  it  would  detract  from  the  generality  of  the 
procedure. 

A  considerable  improvement  of  the  final  result  -  the 
drainage  and  divide  network  -  can  be  expected  with  the 
introduction  of  a  variable-size  grid,  to  vary  automati¬ 
cally  with  variation  in  relevant  input  data,  especially 
of  the  relief. 
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10.  RESULTS  OF  THE  APPLICATION  OF 
A  SIMPLIFIED  VERSION  OF  THE  STANFORD  WATER  - 
SHED  MODEL  TO  NAHAL  YAEL  WATERSHED  05 


by  Harriet  Feldman 


As  a  background  to  the  main  thrust  of  the  research  in 
the  direction  of  a  fully  distributed  model,  it  was  at¬ 
tempted  to  run  a  few  of  the  Nahal  Yael  data  through  a 
simplified  version  of  the  Stanford  Watershed  Model.  The 
basin  selected  was  once  . again  watershed  05. 

The  main  results  aro  listed  in  Tables  10.1,  and  10.2. 

The  rainfall  data  are  listed  in  Table  10.3.  Note  that 
the  printout  lists  erroneously  the  last  nine  months  of 
a  given  water  year  as  if  they  had  the  same  year  as  the 
first  three  months  (October  -  December),  i.e.  "year  68, 
month  12"  should  be  followed  by  "year  69 ,  month  1",  etc. 

It  can  be  seen  that  the  model  functions  fairly  well  on 
an  annual  basis  (the  error  is  15%).  However,  if  one 
breaks  the  results  down  into  a  monthly  level,  a  diffe¬ 
rent  picture  is  obtained.  During  the  48  months  long  pe¬ 
riod  considered,  events  occurred  only  during  5  months. 

Out  of  these: 

Two  months  showed  a  computed  flow  of  zero  for  actual- 

3 

ly  recorded  flow  volumes  (264  and  135  m  ,  respectively) 

3 

One  month  showed  a  computed  flow  of  366  m  for  an  ac¬ 
tual  flow  volume  of  zero; 

Out  of  the  remaining  two  months,  one  had  a  48  per  cent 

3  3 

too  high  computed  value  (974  m  computed  versus  673  m 
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recorded);  and  in  only  one  month  was  the  computed 

3 

flow  volume  reasonably  accurate  (556  m  computed  ver 

3 

sus  573  m  recorded) . 

Tentatively,  the  conclusion  is  that  the  Stanford  Water¬ 
shed  Model,  and  probably  also  similar  models,  are  not 
suitable  for  successfully  simulating  the  mechanics  of 
runoff  as  needed  for  this  (and  simi 1 ar) projects . 


I 


Table  10.1:  Annual  runoff  in  cubic  meters,  recorded 

and  computed  by  a  simplified  Stanford  Wa¬ 
tershed  Model  procedure,  Nahal  Yael  water¬ 
shed  05,  1968/69  -  1971/72 


ANNUAL  RUNOFF  IN  CUBIC  METERS 
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Table  10.2:  Monthly  runoff  in  cubi c.  meters ,  recorded 
and  computed  by  a  simplified  Stanford  Wa¬ 
tershed  Model  procedure,  Nahal  Yael  water¬ 
shed  05,  1968/69  -  1971/72 


YEAR 

MONTH 

COMPUTED 

RECORDED 

68 

11 

0 

264 

68 

01 

97-4 

.  673 

70 

01 

36  6 

0 

70 

03 

0 

135 

71 

12 

556 

573 

V 


£ 
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Table  10.3:  Recorded  daily  rainfall  and  computed  run¬ 
off,  Naha!  Yael  watershed  05,  1968/69  - 
1971/72 


YEAR 

MONTH 

DAY 

RAINFALL  (MM.)  RUNOFF  (MM.) 

68 

11 

20 

0.9 

24 

6.2 

25 

7.2 

6B 

12 

26 

1.5 

68 

1 

21 
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